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EVALUATION 


Current  DoD  plans  indicate  a continued  transition  from  predominantly 
analog  communications  to  digital  communications.  This  has  mandated  a need  to 
re-examine  all  channel  models  for  accuracy  when  applied  to  digital  communica- 
tions signals. 

In  LOS  microwave,  fading  is  usually  assumed  flat.  Availability  calcula- 
tions imply  that  if  a channel  is  calculated  to  have,  for  example,  a 40  dB 
fade  margin,  that  it  will  be  usable  within  40  dB  of  the  average  received 
signal  level.  This  program  has  demonstrated  for  the  first  time  that  group 
delay  is  its  most  severe  at  deep  fades.  If  the  microwave  radio  is  not 
designed  to  handle  large  excursions  in  group  delay,  a high  error  rate  will 
result  several  dB  earlier  than  anticipated.  There  can  thus  be  a wide 
variation  between  calculated  and  actual  availability. 


Host  empirical  data  on  multipath  fading  has  been  taken  for  link  distances 
on  the  order  of  30  miles  or  less.  Links  in  the  Defense  Communication  System 
at  times  exceed  100  miles.  This  program  has  documented,  mathematically, 
conditions  necessary  for  multipath  fading  to  exist  on  long  links. 

Diversity  switching  on  LOS  microwave  is  usually  accomplished  by  a 
selection  switch.  A question  that  always  arises  is  "How  far  should  the 
normal  channel  fade  before  the  standby  channel  is  switched  in?"  Results  of 
this  program  demonstrated  how  group  delay  can  cause  losses  in  bit  count 
integrity  during  diversity  switching  and  proves  that  the  best  strategy  is  to 
switch  to  the  best  channel  at  all  times. 

Much  work  remains  to  be  done  in  further  examining  the  models  developed 
in  this  program  from  a theoretical  point  of  view  and  in  developing  additional 
experimental  data  to  quantify  system  availability  as  the  calculations  apply 
to  digital  conmiunications. 


BRIAN  HENDRICKSON 
Project  Engineer 
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SECTION  1 


1 


SUMMARY 


Continued  evolution  of  the  Defense  Communication  System 
(DCS)  to  an  all  digital  communications  network  and  previous 
indications  that  fading  conditions  on  line-of-sight  (LOS)  micro- 
wave  relay  links  may  exhibit  a high  degree  of  frequency- 
selectivity  have  motivated  the  investigation  reported  here.  The 
report  describes  the  results  of  an  investigation  of  frequency- 
selective  fading  on  LOS  links  to  determine  the  effect  on  digital 
communications  performance.  In  the  modeling  effort  and  con- 
current computation  particular  emphasis  has  been  placed  on  up- 
dating and  revising  previous  work  to  include  diversity  communi- 
cations . In  the  effort  aimed  at  determining  modem  performance 
the  primary  emphasis  has  been  placed  on  determining  the  extra 
SNR  degradation  due  to  the  frequency-selectivity  of  the  multipath 
fading.  Calculations  have  been  performed  for  two  modems,  a base- 
band modem  operating  with  conventional  FM  radios  (DAU  modem)  and 
a narrow  band  phase-continuous  FSK  modem  (IF  modem).  This  sec- 
tion of  the  report  provides  a list  of  our  conclusions  (Section 
1.1  below),  a list  of  our  recommendations  for  further  work 
(Section  1.2),  and  a summary  of  the  report  contents  (Section  1 . 3) . 

1 . 1 Conclusions 


As  a result  of  this  study  the  following  conclusions  have 
been  reached  relative  to  the  effects  of  frequency-selective 
fading  on  DCS  microwave  LOS  links . 

a)  The  largest  group  delay  excursions  occur  during 
the  deepest  fades  (Section  2.3.3).  Consequently, 
diversity  switching  from  one  fading  diversity 

channel  to  another  channel  with  little  or  no  | 

fading  on  it  can  cause  a sudden  large  discon-  j 

tinuity  in  group  delay  depending,  among  other 
things,  on  the  fade  depth  threshold  for  switch- 
ing. Because  of  the  correspondence  between  the  i 

modem  bit  sync  tracking  behavior  and  the  group  i 

delay  (Section  4.2),  and  the  fact  that  calcu-  ‘ 

lated  discontinuities  can  be  comparable  to  and  I 

exceed  the  s)mibol  duration  (Section  4.1),  loss  , 

of  bit  count  integrity  will  occur  in  some 
circumstances . 
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b)  For  diversity  switching,  use  of  the  rule 
"switch  to  other  channel  if  stronger  than 
current  channel,"  is  recommended  over  use 
of  other  rules  that  postpone  the  switching 
until  the  signal  fades  to  a prespecified 
threshold  (Section  4.1). 

c)  According  to  the  theoretical  calculations 
frequency-selective  fading  on  the  100-mile 
Hohenstadt-Zugspitze  link  is  more  extreme 
than  on  any  other  link  we  have  encountered. 
Nevertheless,  it  is  felt  that  some  of  the 
physical  assumptions  (e.g. , regularity  of  layer 
structure  over  the  path)  may  be  near  their 
limits  of  applicability  for  such  a long  link 
(Section  2.4.1).  The  intuitive  notion  that 
such  an  extreme  geometry  can  lead  to  extreme 
fading  is  supported  by  the  program  of  mea- 
surements described  in  [1.1]  wherein  it  was 
reported  that  fading  over  the  Hohenstadt- 
Zugspitze  link  path  is  greater  than  over  any 
other  tested  in  the  DEB  Stage  I system. 

d)  Relative  to  our  evaluation  of  SNR  degrada- 
tion of  modem  performance  due  to  frequency- 
selective  fading,  performance  losses  vary 
widely  with  change  of  physical  conditions. 

In  all  cases  the  SNR  degradation  increases 
as  the  fading  on  the  link  becomes  more 
frequency-selective  (as  indicated  by 
increased  delay  spreads  of  the  interfering 
rays).  Due  to  the  pronounced  frequency- 
selectivity  of  some  fading  situations 
complete  degradation  can  occur  in  a 20  dB 
fade  (see  Tables  1-1  and  1-2  in  Section  1.3.4). 

e)  SNR  degradation  losses,  due  to  frequency- 
selective  fading,  were  quite  similar  for 
the  two  modems  considered  (DAU  modem  and  IF 
modem) . During  conditions  of  mild  fre- 
quency-selectivity our  calculations  indicate 
that  both  modems  perform  well  and  under 
conditions  of  strong  frequency-selectivity 
performance  for  both  is  severely  degraded. 

In  a narrow  transition  region,  the  IF  modem 
has  the  superior  edge  (Section  3) . 
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f)  Whenever  refractive  multipath  fading  takes 
place  on  a link  the  most  probable  situation 
is  that  the  fading  is  due  to  interference  of 
no  more  than  three  rays  (this  is  demon- 
strated in  Section  2.2.4).  In  some  of  our 
computations,  a fourth  and  fifth  ray  have 
appeared  (see  Section  2.4.2).  These  "extra" 
rays,  which  typically  exist  over  only  small 
intervals  of  layer  height,  are  extremely 
close  to  one  another  in  delay  and  are  close 
also  to  the  layer  path  that  is  greatest  in 
delay.  As  a result,  frequency-selectivity 
of  the  link  still  exhibits  a three-path 
character. 


g)  There  are  no  useful  measurements  of  refrac- 
tive index  gradient  taken  in  conjunction 
with  measurements  of  signal  distortion,  i.e., 
the  P-parameters  (Section  2.3.2).  Thus,  a 
valid  comparison  of  model  calculations  with 
experimental  results  is  not  possible. 

Relative  to  our  propagation  calculations, 
however,  comparing  the  character  of  our 
calculated  results  with  the  character  of 
experimental  records  indicates  some  striking 
similarities  (see  Section  2.5)  and  no  discrepancies. 

1.2  Recommendations  for  Further  Work 


To  comprehensively  determine  the  effects  of  frequency- 
selective  multipath  fading  on  estimates  of  system  availability 
it  is  recommended  that  the  following  studies  and  experiments  be 
carried  out: 

a)  Model  Augmentation 

To  increase  the  usefulness  of  the  model  the 
following  tasks  should  be  undertaken: 

• As  pointed  out  in  Section  2.1.2  below, 
the  effects  of  atmospheric  turbulence, 
normally  not  a matter  of  serious  con- 
cern, may  be  magnified  in  a refractive 
multipath  fading  condition.  A study 
of  this  phenomena,  to  determine  its 
effects  on  system  performance,  is 
recommended.  i q 


• Tilted  atmospheric  layers  alter  the  polar- 
ization components  of  waves  that  interact 
with  them,  and  for  this  reason,  their  effect 
on  the  performance  of  "frequency  reuse" 
systems  (those  which  use  two  orthogonally 
polarized  waves  to  approximately  double  link 
capacity)  should  be  determined. 

• To  increase  the  application  of  our  model  to 
a wider  variety  of  physical  situations,  the 
"thin  layer"  case  should  be  included.  This 

s the  case  wherein  the  top  and  bottom  of  the 
atmospheric  layer  are  both  within  the  ver- 
tical region  between  the  two  terminals . 

• Because  of  the  sensitivity  of  multipath 
fading  to  some  of  the  more  important  phy- 
sical parameters  (e.g.,  refractive  index 
gradient  and  layer  height)  and  the  wide  range 
of  values  these  quantities  can  take  on,  sta- 
tistical modeling  of  these  quantities  should 
be  incorporated  into  the  refractive  multipath 
model . 

b)  Probability  of  Antenna  Decoupling 

Prolonged  signal  attenuation  (antenna  decoupling) 
not  exhibiting  the  usual  more  rapid  character 
of  multipath  fading  has  resulted  in  "blackouts" 
of  some  links  in  Germany.  The  analytical  results 
of  this  report  (and  the  companion  computer  pro- 
grams) are  singularly  appropriate, for  calcula- 
tion of  the  relevant  angles.  For  given  or 
postulated  antenna  patterns  we  recommend  a deter- 
mination of  the  lapse  rates  required  to  cause 
antenna  decoupling  and  hence  the  probability  of 
its  occurrence. 

c)  Diversity  Design  Study 

Since  the  study  results  that  have  been  tradi- 
tionally used  to  determine  diversity  spacings 
are  based  on  empirical  data  for  a limited  number 
of  paths,  the  generalization  to  ether  paths, 
other  frequencies , and  other  areas  involves  some 
degree  of  risk  (as  pointed  out  and  summarized  in 
[1.2]).  A design  study  that  takes  into  account 
all  of  the  wide  variety  of  parameters  in  a reasonable 
and  consistent  way  should  be  undertaken.  The  follow- 
ing tasks  and  questions  should  be  addressed: 
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• Further  analysis  of  space  diversity  oper- 
ation, with  special  emphasis  on  an  analy- 
tical determination  of  optimal  antenna 
spacing. 

• All  candidate  diversity  combining  schemes 
should  be  evaluated— especially  new 
schemes  that,  in  addition  to  the  standard 
fade  level,  use  measurements  of  signal 
distortion  in  order  to  take  into  account 
the  extra  SNR  degradation  due  to 
frequency-selective  fadng  and  (when 
appropriate)  the  possible  discontinuities 
during  diversity  switching. 

• What  effect  does  turbulence  [see  (a) 
above]  have  on  the  probability  of  simul- 
taneous deep  fading  on  diversity  branches? 

d)  Availability  Study 

The  probability  that  a link  will  be  available 
is  a basic  parameter  relating  to  link  design. 

We  recommend  an  extensive  effort  to  calculate 
link  availability  via  the  three-step  program 
outlined  below. 

• First,  determine  the  probability  that 
multipath  fading  will  take  place  on  a 
link.  The  results  would  initially  take 
the  form  of  bounds  of  various  quanti- 
ties, e.g.,  lapse  rate  of  refractive 
index  gradient  and  layer  tilt.  Based 
on  available  statistical  measurements 
for  these  quantities  (quite  extensive 
in  the  case  of  lapse  rate,  the  most 
Important  one)  the  results  would  yield 
the  probability  that  multipath  fading 
occurs  on  the  link. 

• Given  the  condition  that  multipath 
fading  occurs,  determine  the  probability 
that  the  P-parameter  values  (see  Sec- 
tion 3)  are  such  as  to  result  in  unaccept- 
able error  rate  (determined  through  cal- 
culations of  modem  performance  such  as  we 
describe  in  this  report). 
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• In  a diversity  context,  given  the  unaccept- 
able performance  of  one  diversity  channel, 
determine  the  probability  that  performance 
of  the  other  diversity  channel  is  also 
unacceptable. 

The  three  probabilities  above,  appropriately  com- 
bined, yield  the  probability  that  the  link  will 
fail  and,  hence,  the  availability. 

e)  Channel  Measurement 


A carefully  planned  program  of  system  function 
and  distortion  parameter  measurement  coupled  with 
refractive  index  measurements  is  needed  so  that 
the  relationship  of  the  multipath  and  frequency- 
selective  character  of  the  link  to  the  link 
parameters,  geometry,  and  atmospheric  conditions 
can  be  pinned  down. 

1 . 3 Summary  of  Report  Contents 

In  this  section  we  summarize  the  report  contents. 

1.3.1  Channel  Model 

A dependable  model  of  the  communication  channel  is  a base- 
line requirement  in  any  attempt  to  evaluate  relative  modem 
performance.  The  channel  of  interest  here  is  the  DCS  LOS  micro- 
wave  relay  channel.  Because  other  causes  of  fading  can  be 
eliminated  by  simple  design  techniques,  the  modeling  effort  in 
this  investigation  centered  on  a phenomena  known  as  refractive 
multipath.  Loosely,  it  can  be  said  that  this  phenomena  derives 
from  the  fact  that  electromagnetic  energy  does  not  travel  in 
straight  lines  in  a medium  characterized  by  a spatially  varying 
refractive  index.  As  it  happens,  the  atmosphere  in  its  usual 
condition  is  characterized  by  a slight,  approximately  linear 
decrease  (or  lapse)  of  refractive  index  with  height.  This  slight 
lapse  is  not,  in  itself,  enough  to  cause  multipath  fading.  On 
occasion,  however,  there  are  regions  of  vertical  extent  char- 
acterized by  a more  rapid  decline  than  usual,  i.e.,  by  a larger 
lapse  rate.  These  regions  are  referred  to  as  "layers"  through- 
out this  report.  When  there  is  the  possibility  that  electro- 
magnetic energy  traversing  a link  may  intersect  with  a layer, 
the  equations  of  geometrical  optics  indicate  that  multiple  ray 
paths  may  exist.  In  such  a case  the  multiple  paths  interfere  in 
the  classical  manner  and  result  in  multipath  fading.  The 
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existence  of  this  condition  depends  on  the  lapse  rate  of  refrac- 
tive index  gradient,  the  path  length,  the  transmitter  and  receiver 
elevations,  the  layer  height,  the  layer  thickness  and  the  layer 
tilt.  In  terms  of  these  parameters  the  equations  of  geometrical 
optics  can  be  used  to  compute  the  number  of  paths,  their  relative 
delays,  and  their  relative  amplitudes,  thus  providing  a complete 
characterization  of  the  refractive  multipath  channel. 

Calculation  of  the  quantities  required  for  this  character- 
ization are  described  in  Section  2.2.  Also  descri  e 
two  basic  modifications  of  the  model  developed  earlier  at  CNR 
[1  31.  One  of  these  is  the  incorporation  of  layer  tilt  into  the 
model  rsection  2.2.3).  The  other  is  the  inclusion  of  the 
physical  situation  illustrated  in  Fig.  1.1(b)  wherein  the  layer 
is  below  the  highest  terminal,  i.e.,  the  highest  terminal  is 
iiKnersed  in  the  layer  (Section  2.2.2).  ^th  of  these  ™°difica- 
tions  were  deemed  necessary  for  the  modeling  of  links  that  are 
extremely  long  and/or  are  characterized  by  an  extrmely  large 
vStt^al  sepaLtion  of  receiver  and  transmitter.  The  correctness 
of  this  conclusion  is  demonstrated  in  connection  with  the 
Hohenstadt-Zugspitze  link  in  Section  2.4.1. 

Though  there  are  several  layer  parameters  affecting  the 
multipath  character  of  fading  on  an  LOS  link  it  has  long  been 
recognized  that  the  lapse  rate  of  refractive  index  gradient 
plays  a dominant  role;  measurements  of  this  ^antity  have  been 
made  on  a world-wide  basis  for  many  years.  The  primary  data 

sources  for  probability  distributions  of  this  par^eter  are  Bean 

and  Dutton  [1.4],  [1.5],  and  Samson^ [1.6] . The  latter  is  of 
particular  use  here  in  that  it  provides  data  for  Stuttgart, 
Germany  that  can  reasonably  be  applied  to 

Zugspitze  link  discussed  in  Section  2.4.1.  Though  the  data 
referred  to  here  constitutes  the  most  comprehensive  source, 
certain  inadequacies  in  its  collection  could  possibly  lead  tr 
slightly  erroneous  conclusions  (see  Section  2.1.1). 

Though  large  lapse  rates  are  known  to 
fading  there  seems  heretofore  to  have  been  no  definitive  state- 
ment !bout  the  relationship  between  lapse  rate  and  the  existen 
of  multipath.  We  have  determined  (see  Section  2.2.4)  an 
existence  condition  for  multipath  due  to  refraction  in  a hori 
zontal  layer.  This  is  given  in  terms  of  the  inequality 
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where  is  the  lapse  rate  of  refractive  index  gradient,  R^  is  the 
radius  of  curvature  of  ray  paths  in  the  normal  atmosphere 
(approximately  5300  feet),  L is  the  link  length  and  A is  the 
height  of  the  layer  above  the  arithmetic  mean  of  the  transmitter 
and  receiver  heights.  Thus,  we  see  that  the  lapse  rate  must 
exceed  some  easily  calculable  value  if  multipath  fading  is  to 
take  place. 

It  should  be  pointed  out  that  (1.1)  holds  only  for  condition 
"a"  of  Fig.  1.1.  Time  did  not  permit  the  determination  of  a 
similar  inequality  for  condition  "B". 

Some  simple  observations  can  be  made  about  the  existence 
condition  above.  First,  as  the  layer  rises  above  the  arithmetic 
mean  of  the  receiver  and  transmitter  heights  a greater  lapse  rate 
is  required  in  order  to  support  multipath.  Secondly,  as  the  link 
length  increases,  less  of  a lapse  rate  is  required  in  order  to 
support  multipath.  These  trends  are  evident  in  the  graphical 
depiction  provided  in  Fig.  1.2. 

One  of  the  most  powerful  applications  of  the  existence 
condition  in  (1.1)  is  in  the  determination  of  the  probability  of 
multipath  fading.  Consider  any  particular  link  for  which  the 
right-hand  side  of  (1.1)  can  be  calculated.  Designate  this  value 
«cQ.  Then  the  probability  that  multipath  fading  exists  on  the 
link  is  just  the  probability  that  exceeds  a^Q.  But,  this  prob- 
ability is  exactly  that  provided  by  the  world-wide  program  of 
measurements  described  above.  Hence,  the  probability  of  multi- 
path  fading  on  the  link  is  quickly  calculable. 

As  reiterated  in  our  recommendations  for  future  work  the 
derivation  of  the  existence  condition  a|)plying  to  condition  "B", 
so  that  an  integrated  view  can  be  obtained,  is  strongly  recom- 
mended. Such  an  analysis  would  provide  upper  bounds  (i.e., 
worst-case  indications)  of  the  probability  of  outage. 
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From  a physical  point  of  view,  the  creation  of  atmospheric 
layers  is  a highly  variable  and  complicated  process,  the  study 
of  which  falls  more  within  the  bailiwick  of  meteorologists  than 
engineers.  Nevertheless,  an  understanding  of  the  basic  mechan- 
isms resulting  in  layer  formation  is  crucial  to  attainment  of  a 
complete  computer  model.  For  this  reason  a discussion  of  weather 
effects  and  layer  formation  is  provided  in  Section  2.1.2  of  this 
report.  We  have  also  provided  in  that  section  a graphic  des- 
cription of  "heating  from  below",  one  of  the  more  common  causes 
of  layer  formation.  The  overall  result  of  heating  from  below  is 
the  formation  of  three  distinct  atmospheric  regions. 

1)  A narrow  surface  region  near  the  ground  charac- 
terized by  strong  gradients  but  slight  turbulence 
(limited  effect  on  communications,  if  any). 

2)  A thick  central  region  characterized  by  strong 
turbulence  due  to  rapid  upward  convection  of 
air.  This  region  constitutes  a well-mixed 
homogeneous  layer  and  is  characterized  by  the 
standard  lapse  rate  of  refractive  index  gradient 
(about  40  Nu/km) . 

3)  A stable  upper  layer  characterized  by  larger 
than  normal  lapse  rates. 

The  intermediate  central  region  in  this  description  is  the 
region  through  which  the  direct  or  pseudo-direct  path  of  Fig.  1.1 
passes.  The  existence  of  turbulence  in  this  region  motivates  con- 
sideration of  adding  to  the  phasor  normally  representing  this 
path  a small  random  phasor  representing  the  turbulence. 

1.3.2  Fading  Characteristics  Due  to  Refractive  Multipath 

That  refractive  multipath  can  occur  with  enough  regularity 
to  warrant  concern  is  beyond  dispute.  The  fact  that  it  is  inher- 
ently a frequency-selective  phenomena  seems  to  be  less  understood. 
The  in-depth  study  described  in  this  report  indicates  that  the 
frequency-selectivity  of  multipath  fading  can  be  of  profound 
importance.  Some  of  the  highlights  of  this  report  relative  to 
the  character  of  fading  on  microwave  LOS  links  (diversity  dis- 
cussion is  deferred  until  Section  1.3.3)  are  summarized  below. 
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A first  step  on  the  road  to  obtaining  an  overview  of  refrac- 
tive multipath  fading  is  to  examine  its  dependence  on  those 
parameters  that  are  best  defined  probabilistically  and  retain  as 
fixed  parameters  those  quantities  for  which  the  probabilities  are 
best  known.  Thus,  one  can  reasonably  ask  "given  a certain  lapse 
rate  for  which  the  probabilities  are  relatively  well  known,  what 
is  the  dependence  of  fade  level  on  layer  height  and  frequency?" 
The  answer,  to  some  degree,  is  provided  by  the  contour  plot  of 
Fig.  1.3  which  corresponds  to  a 30  mile  link  with  the  specific 
link  geometry  shown  and  an  atmospheric  layer  with  a relatively 
strong  lapse  rate  of  refractive  index  gradient  (420  Nu/km) . This 
figure  was  pieced  together  from  12  computer  plots  and  constitutes 
a dual  presentation  of  frequency-selectivity  and  spatial  selec- 
tivity (operating  frequency  on  the  horizontal  axis  and  layer 
height  on  the  vertical  axis).  Essentially,  it  shows  fade  level 
as  a surface  over  the  height- frequency  plane;  each  contour 
depicts  the  combination  of  layer  height  and  frequency  that  will 
result  in  a fade  to  the  depth  with  which  each  contour  is  designated. 

In  discussing  contour  plots  it  is  helpful  to  refer  to  the 
regions  inside  the  smallest  contours  (the  20  dB  contours  in  Fig. 
1.3)  as  the  deep  fade  loci.  It  is  clear  from  the  contour  diagrams 
that  the  deep  fade  loci  are  few  and  far  between;  in  fact,  only  five 
of  them  appear  in  Fig.  1.3.  Additionally , we  note  that  the  10  dB 
fade  contours  extend  over  a few  hundred  MHz,  whereas  the  20  dB 
fade  contours  are  much  smaller.  This  is  in  keeping  with  experi- 
mental data  pointed  out  by  Kaylor  [1.7]  who  has  noted  that  the 
sharp  deep  fades  occur  during  shallower  fades  that  can  extend 
over  a band  approximately  400  MHz  wide. 

Figure  1.3  was  generated  using  operating  frequencies  in  the 
4 GHz  band.  It  is  of  interest  to  compare  these  results  with  the 
results  of  calculations  for  the  8 GHz  band  which  we  provide  in  Fig. 
1.4.  One  conclusion,  borne  out  by  analysis,  is  that  the  fre- 
quency spacings  between  the  deep  fade  loci  in  the  plots  does  not 
depend  on  the  bands  examined,  i.e.,  given  a layer  height  value  the 
deep  fade  loci  occur  with  the  same  frequency  spacing  in  the  4 GHz 
band  as  they  do  in  the  8 GHz  band.  In  fact,  we  have  derived  an 
expression  which,  given  a constant  layer  height  hL,  yields  the 
frequencies  at  which  the  deep  fades  will  occur  (see  Section  2.3.3) 
and  this  expression  indicates  that  at  most  two  deep  fades  may 
occur  every  \AT(hL)]“^  Hz,  where  At  is  the  delay  difference 
between  the  two  rays  refracted  in  the  layer.* 


*0ur  arguments  here  relate  to  multipath  interference  between  a 
total  of  three  rays.  It  is  demonstrated  in  Section  2 . 2 . 4 that  the 
existence  of  a higher  number  of  paths  is  considerably  less 
probable  an  occurrence. 
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REIOTIVG  LfiTCR  HEIGHT  (ET) 

QO  MO  dO  12  0 l&O  200  2N0  2H  0 320  3bO  MOO  iMO  MdO  S20 


PRtQUENCT  CGH2) 

GROUND  RANGE  = 30  MILES 

REFRACTIVE  INDEX  GRADIENT  = -420  Nu/km 

IJVYER  THICKNESS  >62.6  FEET 

HEIGHT  OF  RECEIVER  #1  = 80  FEET  (ABOVE  TRANSMITTER) 

Figure  1.4  Frequency  Constraints  on  Deep  Fade  Loci  (The  deep 
fade  loci  at  any  layer  height  must  fall  on  one  of 
the  dotted  lines.) 
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The  expression  referred  to  above  [ (2 . 156)  in  the  text]  has 
been  used  to  generate  the  dotted  lines  in  Fig,  1,4,  Note  that 
these  lines  pass  directly  through  the  deep  fade  loci,  and  hence 
result  in  a considerable  reduction  of  effort  in  determining  the 
deep  fade  loci  which  are  the  regions  of  greatest  interest, 

A second  conclusion  that  can  be  made  from  a comparison  of 
the  two  contour  diagrams  is  that  the  relative  occurrence  vs , layer 
height  of  the  deep  fades  does  depend  on  the  frequency  band.  Com- 
paring Figs.  1.3andl.4we  see  that  the  vertical  separation  of  the 
deep  fade  loci  in  the  4 GHz  band  is  approximately  twice  what  it  is 
in  the  8 GHz  band.  This  dependence  can  be  demonstrated  generally 
by  appealing  to  the  approximately  linear  dependence  of  ArCh^)  on 
layer  height,  (see  Section  4.1).  When  this  is  done,  we  see  that 
there  can  be  at  most  two  deep  fades  every  [ f At(0) ] "Ihjnax  feet,  i.e., 
the  layer  height  separation  of  the  deep  fade  loci  is  inversely 
proportional  to  f.  Thus,  spatial  selectivity  becomes  more  impor- 
tant at  the  higher  frequencies.  If  it  can  be  assumed  in  the  rough- 
est sense  that  there  is  some  regularity  to  the  layer  motion  then 
one  can  conclude  that  fading  will  be  more  rapid  in  time  at  the 
higher  frequencies . 

It  should  be  emphasized  that  the  contour  diagrams  of  Figs. 

1.3  and  1.4  are  really  just  "snapshots"  of  one  particular  physical 
situation.  The  actual  situation  is  that  the  refractive  index 
gradient  and  layer  height  have  a complicated  temporal  interdepen- 
dence; e.g.,  as  layer  height  and/or  refractive  index  change  in  time 
the  contours  of  Fig,  1.4  change  shape  and  the  deep  fade  loci  move. 
This  motivates  statistical  modeling  of  the  two  basic  physical 
parameters,  layer  height  and  refractive  index  gradient.  A similar 
consideration  should  also  be  made  of  layer  tilt  (more  important 
on  the  longer  links). 

Possible  loss  of  bit  count  integrity  during  diversity  switch- 
ing (to  be  discussed  in  more  detail  in  the  following  section), 
motivates  consideration  of  the  behavior  of  group  delay  during  a 
fading  condition.  As  it  happens,  IF  modem  bit  s3nich  tracking 
closely  follows  the  center-of-band  group  delay,  although  small 
deviations  occur.  This  point  is  illustrated  in  Fig.  1.5  which 
shows  both  group  delay  and  bit  synch  tracking  as  the  angle  of  the 
transfer  function  at  the  center  of  the  band  is  varied  over  360°  (a 
discussion  of  the  significance  of  this  approach  is  provided  in 
Section  3.2.1).  The  correspondence  between  group  delay  and  modem 
tracking  is  even  more  pronounced  in  the  case  of  the  baseband  modem. 

Because  of  the  correspondence  cited  above  it  is  of  interest 
to  examine  the  behavior  of  the  group  delay  in  fading  situations. 
This  topic  is  discussed  in  detail  in  Section  2.3.3.  Here,  we 
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merely  point  out  that  the  group  delay  must  in  general  satisfy  the 
inequality 


N 

L 

i=l 


“i  '^'i 


(1.2) 


where  L is  the  center-of-band  fade  level,  N is  the  number  of  inter- 
fering paths  and  ai  and  Ti  are  respectively  the  amplitude  and  delay 
on  the  i'th  path.  An  important  point  in  this  connection  is  that 
the  upper  bound  in  (1.2)  cannot  be  reached  unless  the  frequency  is 
such  that  the  fade  level  is  at  a minimum.  This  indicates  a close 
correspondence  between  the  deep  fade  loci  and  the  most  extreme 
values  of  group  delay— a correspondence  that  has  been  validated 
without  exception  in  our  plots  of  these  quantities.  Another  point 
is  that  the  group  delay  is  unsubstantial  when  signal  level  is  high, 
L^l.  In  such  a situation  most  of  the  links  we  have  studied 
(extremely  long  links  excluded)  indicate  that  the  group  delay 
magnitude  remains  below  10  ns.  In  a deep  fade,  however,  the  group 
delay  can  be  extremely  large  (see  Section  2.4  and  Section  1.3.4 
below)  and  often  comes  quite  close  to  achieving  the  upper  bound 
in  (1.2). 


1.3.3  Diversity  Operation 

The  fading  characteristics  discussed  above  impact  diversity 
operation  simply  by  virtue  of  the  fact  that  they  differ  from 
diversity  channel  to  diversity  channel.  To  get  the  full  benefits 
of  diversity  operation  one  ideally  would  like  to  have  the  situ- 
ation where  there  is  always  a strong,  clear  signal  on  the  standby 
channel  whenever  the  signal  in  the  primary  channel  is  unaccept- 
ably degraded.  This  situation  corresponds  to  the  case  of  non- 
overlapping deep  fade  loci  when  the  contour  diagrams  for  each 
diversity  channel  are  superimposed.  For  the  most  part,  a reason- 
able choice  of  diversity  operating  parameters  results  in  just 
this  condition.  This  is  shown  in  Fig.  1.6  and  consists  of  the 
20  dB  fade  contour  of  Fig.  1.4  (diversity  channel  #1)  with  the 
20  dB  fade  contours  of  a space  diversity  channel  with  receiving 
antenna  located  30  feet  below  the  primary  receiving  antenna. 

There  is  no  overlap  of  the  20  dB  fade  contours;  i.e.,  although 
some  are  close,  as  with  the  two  pairs  near  20  feet,  simultaneous 
center-of-band  fading  below  20  dB  will  not  take  place.  Note,  how- 
ever, for  a carrier  frequency  (about  8.415  GHz)  half  way  inbetween 
the  two  nodes  at  22  feet,  one  space  diversity  channel  may  exper- 
ience a 20  dB  fade  at  the  lower  end  of  the  band  (14  MHz  assumed) 
and  the  other  may  experience  a 20  dB  fade  at  the  high  end.  One 
concludes,  overall,  however,  that  simultaneous  center-of-band 
deep  fading  on  both  diversity  channels  is  unlikely,  as  it  should 
be.  1_17 


RELATIVE  LAYER  HEIGHT  (FT) 
(ABOVE  HIGHEST  TERMINAL) 


GROUND  RANGE  = 30  MILES 

REFRACTIVE  INDEX  GRADIENT  = -420  Nu/km 

LAYER  THICKNESS  ->62.6  FEET 

Figure  1.6  20-dB  Fade  Contours  on  Both  Diversity  Channels 

(Contours  are  numbered  according  to  diversity  channel) 
Receiver  #1  is  80'  above  transmitter;  Receiver  is 
55'  above  transmitter 
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Space  diversity  operation  is  illustrated  in  Fig.  1.7  for  a 
particular  link  of  interest,  the  Swingate-Houtem  link  spanning 
the  English  Channel.  The  desired  correspondence  of  nulls  on  one 
channel  with  high  signal  levels  on  the  other  is  clear  from  the 
figure.  A detailed  discussion  is  provided  in  Section  4. 

Relative  to  frequency  diversity,  simple  rules,  based  for  the 
most  part  on  the  expression  that  locates  the  deep  fade  frequencies 
(see  Section  2.3.3),  can  be  used  to  determine  the  "optimal"  sepa- 
ration of  carrier  frequencies  in  specific  situations.  As  a numer- 
ical example,  our  expressions  indicate  a frequency  separation  of 
460  MHz  for  the  physical  situation  used  to  generate  the  plot  in 
Fig.  1.8.  The  plot,  which  illustrates  the  frequency  selectivity 
across  the  whole  band  of  permitted  carrier  frequencies  on  the 
Swingate-Houtem  link,  shows  the  reasonableness  of  this  choice;  the 
choice  corresponds  well  with  the  distance  between  peak  and  null. 

To  generate  rules  of  general  utility,  statistical  modeling  is 
required. 

With  the  purpose  of  examining  a physical  mechanism  that  in 
some  circumstances  will  result  in  loss  of  bit  count  integrity,  we 
now  turn  to  a more  microscopic  view  of  the  deep  fade  loci.  Our 
concern  here  relates  to  the  following  question;  what  are  the 
effects  of  diversity  switching  between  space  diversity  channels 
as  a layer  moves  through  a deep  fade  locus? 

We  consider  three  simple  diversity  rules: 

1.  Switch  to  other  channel  if  current  channel  fades 
to  depth  greater  than  30  dB  and  other  channel  is 
stronger . 

2.  Switch  to  other  channel  if  current  channel  fades 
to  depth  greater  than  20  dB  and  other  channel  is 
stronger. 

3.  Switch  to  other  channel  if  stronger  than  current 
channe 1 . 

Recall  from  our  discussion  of  Section  1.3.2  that  the  largest 
excursions  of  group  delay  tend  to  occur  during  the  deepest  fades. 
This  implies  that  diversity  switching  from  a diversity  channel 
with  a deep  fade  on  it  to  another  channel  with  little  or  no  fade 
is  equivalent  to  switching  from  a channel  with  a large  group 
delay  to  a channel  with  very  small  group  delay.  Because  of  the 
rapid  switching  on  digital  links,  there  can  be  a large  sudden 
discontinuity  in  delay  and  hence  a loss  of  bit  count  integrity. 

For  concreteness  we  look  at  the  fade  in  Fig.  1.9  and  assume 
upward  layer  motion*  and  initial  operation  of  diversity  channel 
#1.  In  Fig.  1.10  we  have  presented  the  group  delay  plot 
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CARRIER  FREQUENCY  = 4.808  GHz 


Figure  1.7 


Theoretical  Space  Diversity  Fade  Record  on 
Swingate-Houtem  Link 
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Figure  1.8  Fade  Level  Vs.  Frequency  Over  Total  Band  for 
Horizontal  Layer  on  Swingate-Houtem  Link 
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Figure  1.9  Fade  Depth  Vs.  Layer  Height  on  Hohenstadt-Zugspitze 
Link 
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(ABOVE  TRANSMITTER) 

LAPSE  RATE  - 400  Nu/km 

LAYER  SLOPE  - .009 

CARRIER  FREQUENCY  - 8.301  GHz 


Figure  1.10  Group  Delay  Vs.  Layer  Height  on  Hohenstadt-Zugspitze 
Link 
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corresponding  to  Fig.  1.9  and  marked  off  the  points  at  which 
switching  takes  place  according  to  the  diversity  rules  listed 
above.  For  rule  1 we  get  a discontinuity  of  177  ns.  For  rule  2 
a 50  ns  discontinuity  occurs,  whereas  for  rule  3 it  can  be  shown 
on  an  expanded  scale  that  the  discontinuities  are  considerably 
less  than  10  ns. 

Relative  to  a symbol  length  (72.3  ns)  the  group  delay  dis- 
continuity that  occurs  when  rule  3 is  used  is  extreme.  In  view 
of  the  close  correspondence  between  modem  tracking  and  group 
delay  (Section  4.2)  it  can  be  concluded  that  switching  during  a 
fade  to  30  dB  (or  any  deeper  fade)  will,  in  some  circumstances, 
result  in  loss  of  bit  count  integrity  on  this  link. 

1.3.4  Modem  Performance  Degradation 

Performance  in  the  presence  of  representative  multipath 
conditions  has  been  evaluated  for  two  modem  designs:  a baseband 
modem  designed  for  existing  FM  radios,  and  a phase-continuous 
FSK  modem  which  interfaces  to  existing  radios  at  IF.  Both 
modems  achieve  a packing  density  of  2 Bps/Hz;  the  latter  is  a 
higher  performance  modem  in  the  sense  that  it  requires  a smaller 
Etj/No  to  achieve  a given  error  rate. 

For  each  of  27  multipath  conditions  calculated  for  the 
Hohenstadt-Zugspitze  link  performance  of  the  modems  was  calcu- 
lated at  three  fade  depths:  10,  20,  and  30  dB.  As  the  depth  of 
the  fade  increases,  frequency-selective  effects  become  more 
pronounced.  In  other  words,  fading  will  be  "flat"  (that  is 
undistorted)  for  smaller  fade  depths  and  frequency-selective 
(that  is  distorted)  for  larger  fade  depths.  At  some  fade  depth 
distortion  will  reach  the  point  where  the  eye-pattern  closes  for 
certain  symbol  sequences;  when  this  happens  it  will  be  impos- 
sible to  obtain  a satisfactory  error  rate  regardless  of  the 
value  of  E^/No  available  to  the  receiver.  This  behavior  may  be 
seen  qualitatively  in  Tabic  j.-1  which  summarizes  baseband  modem 
performance  for  each  of  the  channel  conditions  considered,  or 
in  Table  1-2  which  app’^ies  to  the  IF  modem.  (In  these  tables 
the  three  layer  heights  which  were  considered  for  each  combina- 
tion of  refractive  index  gradient  and  layer  tilt  slope  are 
lumped  together  for  simplicity.)  Selecting  a gradient  and 
slope,  one  can  look  across  horizontally  and  compare  performance 
at  fade  depths  of  10,  20,  and  30  dB.  For  relatively  flat  chan- 
nels the  transition  from  flat  fading  to  frequency-selective 
fading  occurs  at  very  deep  fades,  while  for  more  selective 
channels  the  transition  occurs  for  shallower  fade  depths.  Since 
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TABLE  1-1 

lUALITATIVE  FADING  BEHAVIOR  OF  BASEBAND  MODEM 
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TABLE  1>2 

QUALITATIVE  FADING  BEHAVIOR  OF  IF  MODEM 


fade  depth  is  quantized  to  10  dB  increments  in  this  table,  the 
transition  from  flat  fading  to  complete  degradation  is  usually 
abrupt.  Individual  performance  curves  on  which  the  tables  are 
based  may  be  found  in  Section  3.2. 

Summarizing,  our  evaluation  of  SNR  degradation  of  modem 
performance  due  to  frequency-selective  fading  reveals  that  the 
performance  losses  can  vary  widely  with  change  of  physical 
conditions.  In  all  cases  the  SNR  degradation  increases  as  the 
fading  on  the  link  becomes  more  frequency-selective.  Due  to 
the  increased  frequency-selectivity  of  some  physical  situations 
(see  Tables  1-1  and  1-2)  complete  degradation  can  occur  in  a 
20  dB  fade.  Increased  frequency-selectivity,  of  course,  cor- 
responds, to  increased  delay  spreads  of  the  ray  paths,  and 
the  trend  toward  greater  frequency-selectivity  corresponds  gener- 
ally with  increases  in  lapse  rate  of  refractive  index  gradient  and , 
ad-hoc,  increases  in  layer  tilt.  (As  an  example  we  point  out 
that  layer  slope  variations  from  .007  to  .010  in  Table  1-1, 
corresponding  to  a refractive  index  gradient  of  400  Nu/km, 
are  accompanied  by  delay  spreads  of  the  layer  rays  ranging 
from  about  2.4  ns  to  4.4  ns.)  The  parallel  trends  of  increased 
frequency-selectivity  and  increased  SNR  degradation  due  to 
signal  distortion  are  clear  from  the  tables.  (See  Section  3 
for  detailed  discussion.) 

Comparing  Tables  1-1  and  1-2  it  may  be  seen  that  the 
effects  of  multipath  on  the  two  modems  are  remarkably  similar. 
However,  a careful  examination  of  the  tables  or  of  the  figures 
in  Section  3 (Figs.  3.6  through  3.26)  shows  that  the  FSK  modem 
is  somewhat  less  susceptible  to  multipath  than  the  baseband 
modem. 
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SECTION  2 


CHANNEL  MODELING 


In  this  section  we  describe,  analyze,  and  examine  the 
causes  and  effects  of  frequency-selective  fading  due  to 
refractive  multipath  on  microwave  line-of-sight  links.  A major 
goal  of  the  effort  described  here  is  to  determine  the  eventual 
effects  on  modem  performance.  Because  of  its  separate  import- 
ance, however,  discussion  of  this  topic  is  deferred  until 
Section  3. 

In  Section  2,1  below  we  describe  the  propagation  channel. 
The  viewpoint  in  Section  2.1  is  primarily  a physical  one.  We 
discuss  the  effects  of  weather  on  the  formation  of  layers  and 
emphasize  the  importance  of  certain  layer  parameters. 

In  Section  2.2  we  describe  our  analysis  of  the  propagation 
channel.  This  includes: 

1)  A description  of  the  link  and  the  setting  up 
of  the  problem. 

2)  Solution  of  the  problem  using  geometrical 
optics  with  emphasis  placed  on  analyzing 
the  situation  wherein  the  layer  is  below 
the  highest  terminal. 

3)  Inclusion  of  layer  tilt  into  the  model. 

4)  Lapse  rate  requirements  for  the  support  of 
multipath . 


It  is  from  the  results  of  this  analysis  (yielding  all  of  the 
ray  path  delays  and  amplitudes)  that  one  can  compute  any  of  the 
system  parameters  of  interest. 
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In  Section  2.3  we  describe  the  system  parameters, fade 
depth  and  group  delay,  and  present  in  general  terms  a dis- 
cussion of  the  P and  T parameters  that  play  a pivotal  role  in 
both  the  representation  of  signal  distortion  and  the  modem 
performance  calculations  of  Section  3.  Additionally,  we  des- 
cribe in  this  section  a significant  relation  between  the  fre- 
quencies at  which  the  fade  depth  and  group  delay  extrema  occur. 

Two  important  communication  channels  of  interest,  the 
Hohenstadt-Zugspitze  link  and  the  Swingate-Houten  link  are 
analyzed  in  detail  in  Section  2.4.  Fading  characteristics  of 
the  channel  are  depicted  via  plots  of  the  multipath  delay  and 
amplitude  profiles  and  plots  of  fade  depth  and  group  delay  vs. 
layer  height  and  frequency.  Discussion  of  diversity  effects 
is  deferred  until  Section  4. 

Certain  propagation  parameters  relating  to  the  Hohenstadt- 
Zugspitze  link  have  been  used  in  the  modem  performance  calcula- 
tions presented  later  in  Section  3.  The  way  in  which  these 
parameters  have  been  selected  in  order  to  reduce  them  to  a 
manageably  small  set  is  described  in  Section  2.4.3. 

In  Section  2.5  we  present  a brief  review  of  the  litera- 
ture. The  emphasis  here  is  on  relating  our  computed  results 
to  the  available  data  that  has  been  accumulated  by  others. 

! 

[ 2.1  Description  of  Propagation  Channel 

i 

i In  this  section  we  describe  the  physical  character  of  the 

LOS  propagation  channel.  The  basic  nature  of  multipath  fading 
! on  LOS  links  is  described  in  Section  2.1.1.  Section  2.1.2 

describes  the  effect  of  weather  conditions  and  the  physical 
mechanism  governing  the  formation  of  layers.  S -action  2.1.3 

addresses  some  questions  related  to  layer  height  and  layer 
I thickness. 
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2.1.1  Morphology  of  Refractive  Multipath 

In  this  section  we  describe  the  mechanism  of  refractive 
multipath  and  the  physical  conditions  which  cause  it. 

Generally,  it  can  be  said  that  refractive  multipath 
arises  because  of  the  bending  of  electromagnetic  rays  in  an 
inhomogenous  medium.  The  equations  of  geometric  optics 
indicate  that  electromagnetic  rays  travel  in  straight  lines 
in  a perfectly  uniform  atmosphere.  In  the  troposphere, 
however,  the  atmosphere  is  often  not  uniform;  in  fact,  be- 
cause of  its  sandwiched  location  between  the  earth  and  the 
stratosphere,  there  may  be  rapid  changes  in  temperature, 
pressure,  and  vapor  pressure  with  height.  These  gradients, 
in  turn,  give  rise  to  a gradient  in  refractive  index.  In 
fact,  an  acceptable  order  of  magnitude  of  the  variation  M 
of  the  index  of  refraction  is  obtained  from  the  formula, 

^«/ip-^+7Ar  (2.1) 

where  hP  is  the  pressure  variation  in  millibars,  AT  is  the 
variation  of  the  temperature  in  degrees  Kelvin,  and  Ar  is 
the  variation  of  the  mixing  ratio,  defined  as  the  mass  of 
water  vapor  (in  grams)  mixed  with  1 kg  of  dry  air  [2.1]. 
Measurement  of  each  of  the  constitutive  parameters  in  (2.1) 
is  sometimes  used  for  indirect  measurement  of  N.  (See,  for 
example,  the  description  of  radiosonde  measurements  in  [2.2].) 

In  its  most  commonly  occurring  state,  the  troposphere 
is  characterized  by  a mild  refractive  index  gradient;  the 
refractive  index  decreases  by  about  40  Nu/km.  In  such  a 
case,  we  say  that  the  lapse  rate  s£  is  given  by 

s£  = 40  Nu/km  (2.2) 

(This  is  the  numerical  value  used  for  the  "standard"  or 
"direct  path"  atmosphere  in  the  remainder  of  this  report.*) 

On  some  occasions,  however,  anomalous  weather  conditions  can 
exist  which  result  in  the  formation  of  refractive  layers. 


It  will  be  recalled  that  N units  are  just  the  usual  dimen- 
sional units  of  refractive  index  multiplied  by  10^  (see 
[2.3,  p.  3-11]). 
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when  this  happens,  there  is  the  usual  "standard"  lapse  of 
40  Nu/km  until  some  unpredictable  height,  referred  to  as  the 
"layer  height",  is  reached,  at  which  point  a much  stronger 
lapse  (e.g.,  300  Nu/km)  occurs  until  the  "top"  of  the  layer 
is  reached.  At  this  point  the  milder  lapse  characterizing 
the  standard  atmosphere  takes  over. 

The  vertical  region,  characterized  by  the  larger  lapse, 
is  referred  to  as  the  "layer".  The  extent  of  this  region  is 
referred  to  as  the  "layer  thickness".  The  vertical  point 
(measured  with  respect  to  some  convenient  reference]^  at  which 
the  stronger  larger  lapse  begins,  is  referred  to  as  the 
"layer  height". 

There  are  many  complications  to  the  picture  described 
above.  Not  the  least  of  these  is  the  abundance  of  layer 
parameters  which  have  not  been  measured  extensively.  Dis- 
cussion of  this  topic,  because  it  impacts  our  computer  model, 
is  deferred  until  Section  2.1.3.  Another  complication  arises 
because  the  troposphere  may  not  separate  in  simple  fashion 
into  well-defined  layers  of  constant  lapse  rate.  Additionally, 
the  troposphere  may  vary  linearly  and/or  nonlinear ly  in  the 
horizontal  direction.  Although  the  possibility  of  such  con- 
ditions actually  occurring  can  in  no  way  be  ruled  out,  con- 
struction of  a model  based  on  such  considerations  would  be 
of  little  use  in  view  of  the  all-important  fact  that  there 
would  be  no  available  empirical  inputs  to  the  model.  Addi- 
tionally, the  use  of  the  refractive  index  gradient  as  a valid 
description  of  tropospheric  propagation  is  widely  accepted. 

In  fact,  the  entire  worldwide  effort  to  measure  refractive 
index  gradient  is  predicated  on  the  assumption  that  a linear 
or  near- linear  dependence  of  this  parameter  on  height  is  a 
dominant  influence  on  tropospherically-propagated  waves.  If 
this  were  not  so,  the  vast  array  of  data  collected  (see  an 
example  in  Figure  2.2)  would  be  of  little  use.  The  fact  that 
the  premise  is  supported  by  the  stable  geometry  of  large  air 
masses  and  frontal  inversions  is  discussed  in  Section  2.1.2. 

The  basic  laws  of  electromagnetic  wave  propagation 
indicate  the  propagation  of  waves  in  curved  paths  when  a 
refractive  index  gradient  is  present.  In  fact,  when  the 
gradient  is  a constant,  the  rays  travel  in  circular  arcs 
[2.4].  For  point-to-point  LOS  propagation  on  the  earth's 
surface,  these  arcs  may  be  concave-upward  or  concave- 
downward  depending,  respectively,  on  whether  the  refractive 
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index  increases  or  decreases  with  height.  A negative  refrac- 
tive index  gradient  is  by  far  the  more  common;  thus,  a ray 
typically  tends  to  travel  in  concave -downward  arcs.  This 
does  not  mean  that  the  ray  necessarily  tends  to  bend  back  to 
earth;  the  ray,  though  concave -downward,  may  have  less  curva- 
ture than  the  earth's  surface  and  effectively  bend  away  from 
the  earth.  A point  of  importance  in  this  connection  is  the 
following:  even  though  a ray  path  may  be  concave-downward, 

as  illustrated  in  Figure  2.1(a),  a plot  of  ray  height  vs. 
distance  x along  the  path  may  exhibit  an  upward  concavity 
as  illustrated  in  Figure  2.1(b).  From  an  engineering  point 
of  view,  this  latter  representation  is  the  one  of  importance. 
Analytically,  it  is  achieved  via  an  "earth-flattening"  coor- 
dinate transformation  [2.3].  The  path  resulting  from  this 
transformation,  like  that  in  Figure  2.1(b),  shall  be  called 
the  virtual  path.  Using  the  convention  that  a negative  radius 
indicates  an  upward  concavity  and  a positive  radius  indicates 
a downward  concavity,  we  have  from  [2.3]  that  the  virtual 
radius  in  miles  is  given  by 


0.62139  x 10^ 
=£  - 156.78 


(2.3) 


where  ^ is  the  lapse  rate  of  refractive  index  gradient  in 
Nu/km,  i.e.,  is  the  rate  at  which  the  refractive  index  de- 
creases with  height. 

From  the  expression  for  virtual  radius  of  curvature  in 
(2.3),  we  have  constructed  Table  2-1.  In  this  table  we  use 
the  descriptive  terms  "subrefraction",  "super-refraction", 
etc.,  that  have  been  used  extensively  in  the  literature. 

Each  category  corresponds  to  a range  of  lapse  rate,  and  we 
have  conformed  with  [2.5]  and  [2.6]  in  our  definitions  ... 
with  one  exception.  Because  of  a gap  that  would  occur  be- 
tween «£=0  and  id -100,  we  have  bestowed  the  title  "normal- 
refractive"  on  layers  having  lapse  rates  in  that  range. 

Our  major  concern  in  this  report  is  with  layers  in  the 
last  category.  These  are  the  "ducting"  layers  that  give 
rise  to  refractive  multipath.  Extensive  programs  of  measure- 
ment of  lapse  rates  which  have  been  and  are  being  carried 
out  on  a worldwide  basis  indicate  that  ducting  layers  can 
and  do  occur  with  sufficient  frequency  to  warrant  concern. 
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(a)  Actual  Path 

Concave -downward  Ray  Path  Above  Earth 


(b)  Virtual  Path 

Plot  of  y vs . X,  where 

y s height  of  ray  path  above  earth 
X = distance  along  earth  from 
center  point 


Figure  2.1  Comparison  of  Actual  Path  and  Virtual  Path 
Geometries 
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This  was  clearly  evident  in  the  early  comprehensive  data  of 
Bean,  et  al.  [2.5]  and  from  several  measurement  programs  that 
have  been  described  since  [2.6],  [2.7]. 

A wide-ranging  collection  of  probability  distributions  is 
gathered  together  in  [2.6].  One  of  particular  interest,  since 
it  relates  to  our  calculations  on  the  Hohenstadt-Zugspitze  link 
[ see  Section  2.4.1]  is  the  data  for  Stuttgart,  Germany,  shown  in 
Figure  2.2.  This  data  shows  that  .1%  of  the  observations  made 
in  Stuttgart  indicated  that  the  lapse  rate  exceeded  300  Nu/km. 
As  pointed  out  in  [2.8],  this  data  is  justifiable  extrapolated 
to  the  Hohenstadt-Zugspitze  link.  Unfortunately,  the  informa- 
tion in  Figure  2.2  does  not  constitute  a comprehensive  data 
source.  The  reason  for  this  is  that  measurements  were  only 
taken  two  times  a day  (at  the  same  time  of  day)  for  several 
months.  If  nothing  more,  the  data  certainly  indicates  the 
occasional  occurrence  of  large  lapse  rate,  whereas  some  places 
in  the  world  never  show  lapse  rates  in  excess  of  150  Nu/km. 

(See  data  for  Chita  in  the  Soviet  Union.)  In  such  regions 
the  phenomena  of  ducting  rarely,  if  ever,  occurs.  In  other 
regions,  it  appears  to  occur  with  substantial  regularity. 

2.1.2  Weather  Effects  and  Formation  of  Layer 

Refractive  layers  cannot  occur  unless  there  are  gradients 
in  pressure  and/or  temperature  and/or  moisture  content  of  the 
atmosphere.  From  a consideration  of  all  the  available  measure- 
ments, Craig,  et  al.  [2.9]  have  adopted  the  following  simple 
expression  for  numerical  use; 

(n-l)xlO^  = §(Pd'^^'^¥)  “ t(p‘^¥) 

where 

n = the  index  of  refraction  of  moist  air 
c = 79  K/mb 
b = 4800°K 

. o 

T = temperature  in  K 
p = total  pressure  in  mb 

e = the  partial  pressure  of  water  vapor  in  mb 

When  n is  replaced  by  the  modified  refractive  index,  the 
left-hand  side  of  (3.4)  becomes  a quantity  that  is  called 
the  refractive  modulus  M. 
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PERCENT  OF  TIME  ORDINATE  IS  NOT  EXCEEDED 


Note: 


k = -R/a  where  R is  the  radius  of  curvature 
in  text  and  a is  the  earth's  radius. 


Figure  2.2  Distribution  of  Radiosonde  Data  for 

Stuttgart,  Germany  (Federal  Republic) 
From  [2.6] 
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EFFECTIVE  EARTH  RADIUS  FACTOR,  k 


Pressure  varies  so  slowly  with  height  that  its  contribu- 
tion to  (2.4)  can  be  discounted.  On  the  other  hand,  the 
humidity  variation  is  an  important  part  of  The  water 

vapor,  evaporating  from  the  earth,  tends  to  diffuse  in  a 
vertical  upward  movement,  distributing  iteself  over  the  entire 
atmosphere.  In  the  absence  of  acmos; '^eric  turbulence,  the 
diffusion  proceeds  very  slowly  and,  because  of  this  stability, 
e may  have  a very  strong  vertical  gradient. 

The  separate  dependence  of  M on  temperature  change  and 
change  in  vapor  pressure  is  illustrated  in  Figure  2.3.  We 
see  that,  in  the  absence  of  any  vertical  gradient  of  vapor 
pressure,  a temperature  inversion  (temperature  increasing 
with  height)  tends  to  cause  M to  increase  less  rapidly  with 
height  than  in  the  standard  case.  If  the  temperature  inver- 
sion is  strong  enough,  it  causes  an  actual  decrease  of  M with 
height.  In  the  absence  of  any  vertical  temperature  gradient, 
the  vertical  M gradient  varies  in  direct  proportion  to  the 
variation  in  the  vertical  gradient  of  vapor  pressure. 

There  are  a multitude  of  meteorological  conditions  that 
can  result  in  a formation  of  a super-refracting  layer,  i.e., 
that  can  result  in  M inversions.  These  are: 

(1)  Heating  from  below 

(2)  Cooling  from  below 

(3)  Subsidence  inversions 

(4)  Frontal  inversions 

Craig  and  Montgomery  [2.10]  have  presented  a graphic 
description  of  "heating  from  below",  one  of  the  more  common 
causes  of  layer  formation.  Following  their  scenario,  we 
examine  the  sequence  of  events  occurring  on  a calm  clear 
morning  as  a result  of  the  ground  being  warmed  by  solar 
radiation.  Near  the  initially  cold  ground  there  is  an 
initially  stable  air  mass  but,  as  the  ground  heats  up, 
packets  of  air  warmed  by  contact  with  the  ground  begin  to 
rise  convectively . These  packets  accelerate  at  first  and 
then  begin  to  declerate  at  a level  that  is  well-defined  from 
thermodynamic  considerations,  i.e.,  "there  is  a definite  top 
to  the  convection".  The  overall  result  of  heating  from  below 
is  the  formation  of  three  layers: 

(1)  A narrow  surface  layer  near  the  ground  charac- 
terized by  strong  gradients  but  slight  turbulence 
(air  packets  are  just  starting  their  acceleration) 
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(2)  A thick  central  layer  characterized  by  strong 
turbulence  due  to  rapid  convection  of  the  air 
packets 

(3)  An  upper  layer  that  is  stable. 

The  ground  layer  which  extends  to  no  higher  than  about 
50  feet  is  unlikely  to  have  degrading  effects  on  links  other 
than  those  that  have  very  low-lying  antenna  systems.  The 
second  turbulent  layer  constitutes  a well-mixed  homogenous 
layer  and  is  characterized  by  the  "standard"  refractive 
index  gradient  (about  -40  Nu/km) . The  highest  layer  is 
stable,  but  is  comm.only  characterized  by  a lapse  of  water 
vapor  concentration.  To  quote  Craig  and  Montgomery  [2.10] 

"The  stable  layer  is  therefore  usually  superstandard."* 

This  observation  derives  from  the  data  of  Figure  2.3(b). 

The  features  that  the  above  description  have  in  common 
with  our  computer  model  illustrated  later  in  Figure  2.5  — 2.7 
should  be  apparent. 

In  their  description  of  heating  from  below,  the  authors 
go  on  to  discuss  the  horizontal  variations  in  the  boundary 
heights  between  the  layers  due  to  wind  movement  of  cool  air 
masses  over  warm  masses  of  water.  This  could  result  in 
tilting  of  the  layer,  as  we  discuss  in  Section  2.2.3. 
Additionally,  it  is  easy  to  see  that  layer  tilt  could  arise 
from  variability  of  earth  moisture  content  over  a link.  A 
possibility  not  mentioned  in  the  reference  cited  above  would 
be  a regular  horizontal  variation  in  refractive  index  profile 
over  a long  link  due  to  the  different  degrees  of  heating  that 
the  link  submits  to  as  the  sun  rises  above  it.  This  mechanism 
of  layer  tilt  is  particularly  reasonable  in  view  of  the  fact 
that  heating  from  below  is  usually  an  early  morning  phenomenon. 

One  interesting  feature  of  the  scenario  discussed  above 
is  the  "fully-developed  turbulence"  in  the  intermediate 
homogenous  layer.  This  layer,  which  is  called  the  "direct 
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'^Meteorologists,  for  purposes  of  easy  reference,  often  break 
down  the  several  situations  we  have  categorized  in  Table  2-1 
into  standard,  substandard,  and  superstandard.  By  "sub- 
standard" they  refer  to  the  condition  wherein  the  refractive 
index  increases  more  rapidly  with  height  than  in  the  standard 
case.  By  superstandard  they  refer  to  a greater  rate  of 
decrease . 
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path'"  or  "pseudo-direct  path"  in  our  model  (see  Section  2.2.1) 
has  not  been  investigated  as  a turbulent  (scintillating)  channel. 
As  mentioned  above,  the  region  below  the  stable  layers  causing 
refractive  multipath,  though  characterized  on  the  average  by  the 
lapse  rate  for  the  "standard"  atmosphere,  may  exhibit  strong 
turbulence.  Turbulence,  or  randomness  of  the  refractive  index, 
has  long  been  a topic  of  interest  for  experts  on  atmospheric 
propagation.  The  general  consensus,  however,  is  that  the  rms 
values  of  amplitude  fluctuations  are  so  small  as  to  not  warrant 
undue  consideration  as  far  as  communication  system  performance  is 
concerned;  typical  estimates  [2.3]  are  that  the  rms  value  of  the 
randomly  fluctuating  phasor  representing  the  turbulence  is  20  dB 
down  from  the  direct  path.  Consider  the  implications  of  this 
report,  however;  classical  interference  of  multipath  rays  can, 
with  a high  enough  probability  to  warrant  concern,  yield  a 
resultant  phasor  (signal  level)  that  is  20  dB  down  from  the 
direct  path  signal  level.  Now,  consider  a random  phasor  (rms 
value  20  dB  down  from  the  direct  path)  that  is  added  to  the 
phasor  representing  the  direct  path.  By  superposition  this 
phasoi  can  be  added  to  the  resultant  phasor  of  the  classical 
multipath  interference.  We  find  ourselves  in  the  position  of 
adding  a random  phasor  with  rms  value  20  dB  down  from  the  direct 
path  to  a small  "deterministic"  phasor  of  approximately  the  same 
amplitude.  The  effects  can  be  dramatic.  Whereas,  turbulence 
can  normally  be  neglected,  a deep  fade  condition  due  to  refrac- 
tive multipath  can  greatly  magnify  its  effects.  It  is  recommended 
that  this  condition,  probable  in  a meteorological  context,  be 
investigated  more  fully. 

The  phenomenon  of  subsidence  is  also  known  to  be  a cause  of 
superstandard  refractive  conditions;  this  phenomenon  is  the  sink- 
ing or  lowering  of  large  masses  of  air  from  high  to  low  levels. 

For  reasons  described  by  Craig  [2.11],  this  condition  tends  to 
result  in  stable  layers  below  which  exists  a region  wherein  air 
is  mixed  by  turbulence  or  convection.  The  base  of  the  layer 
(subsidence)  is  often  characterized  by  a subsidence  inversion 
wherein  the  humidity  rapidly  decreases  with  height;  this  latter 
situation  arises  because  the  air  above  the  inversion,  having 
descended  from  higher  levels,  is  usually  very  dry.  Because  of 
the  temperature  increase,  v\?hich  is  always  present,  and  the 
humidity  decrease,  which  is  frequently  present,  there  is  prac- 
tically always  a superstandard  M gradient  in  a subsidence 
inversion . 

Fronts  can  have  important  effects  on  temperature,  humidity, 
and  M profile,  but  because  their  effects  are  more  variable  than 
those  of  subsidence,  meteorologists  are  reluctant  to  make  general 
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statements  about  the  effects  they  will  have  on  propagation.  A 
front  (or  frontal  surface)  is  the  boundary  between  two  large  air 
masses  and,  because  the  overriding  air  mass  is  the  warmer  one, 
there  is  always  a stable  layer  and  often  a temperature  inversion. 
Unlike  the  interaction  between  two  large  air  masses  in  a sub- 
sidence inversion,  the  overriding  air  is  moist.  Some  information 
is  available  on  slopes  of  frontal  surfaces,  however  (see  Section 
2.2.3  below).  The  same  physical  mechanisms  governing  the  tilt 
of  fronts  govern  the  tilt  of  subsidence  boundaries  as  well  and, 
in  lieu  of  more  definitive  data,  have  been  used  in  our  computer 
model . 

2.1.3  Additional  Layer  Parameters 

In  this  section  we  discuss  two  additional  layer  parameters, 
the  layer  thickness,  and  the  layer  height. 

Dependence  of  multipath  fading  on  layer  thickness  has  an 
"on  or  off"  character  to  it . A layer  is  either  thick  enough  to 
trap  a given  ray  or  is  so  thin  that  the  ray  escapes  from  the  top 
of  the  layer  and  bends  away  from  the  earth.  A record  of  fade 
level  vs.  layer  thickness  (all  other  parameters  constant)  would 
display  sharp  discontinuities  due  to  this  phenomenon. 

Some  measurements  of  layer  thickness  have  been  made.  Statis 
tical  distributions  of  layer  thickness,  based  on  measurements 
over  a three-year  period  in  the  Arctic  (Fairbanks,  Alaska),  in  a 
temperature  climate  (Washington,  DC),  and  in  a tropical  maritime 
climate  (Swan  Island),  are  reported  by  Bean  and  Dutton  [2.2]. 

This  data  is  reproduced  in  Figure  2.4.  Note  that  one  can  make 
the  rough  conclusion  from  this  figure  that  layers  with  thickness 
in  excess  of  280  meters  were  observed  TU  of  the  time  in  the 
temperature  region. 

Beckmann  and  Spizzichino  [2.1]  have  also  presented  data  on 
layer  thicknesses.  Their  rough  deductions  from  measurement  of 
the  cutoff  wavelength  indicate  that  layers  with  thicknesses  in 
excess  of  100  m may  be  present  more  than  107o  of  the  time. 

Unfortunately,  there  does  not  seem  to  be  data  available  on 
the  joint  statistics  of  layer  thickness  and  refractive  index 
gradient;  i.e.,  we  cannot  gain  any  insight  into  the  frequency  of 
occurrence  of  thick  layers  with  large  lapse  rates. 

Layer  height  is  another  parameter  for  which  the  data  is  very 
limited.  In  fact,  because  there  seems  to  be  no  useful  data  avail 
able  for  layer  heights,  we  have  tended  to  view  this  parameter  as 
an  independent  variable  against  which  to  plot  the  results  of  our 
calculations.  There  are  reasons  to  suspect  that  the  height  of  the 
layer  is  a highly  variable  dynamic  quantity.  In  fact,  Bullington 
[2.12]  has  come  up  with  estimates  of  vertical  velocity  for  layers 
(a  few  feet  per  minute).  2-14 
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2 . 2 Analysis  of  Propagation  Channel 


In  this  section  we  describe  our  analysis  of  refractive 
multipath  on  LOS  links.  The  analysis,  of  course,  draws  on 
earlier  analytical  approaches  presented  in  a previous  CNR  report 
[2.13].  Comparisons  of  the  approach  to  the  problem  presented 
here  and  in  the  earlier  report  will  manifest  obvious  similarities 
in  the  general  way  that  the  problem  is  approached.  The  differ- 
ence, here,  is  that  we  have  analyzed  the  situation  wherein  the 
interface  of  the  atmospheric  layer  and  the  "standard”  atmosphere 
is  at  a vertical  location  below  the  highest  antenna  terminal. 

This  is  an  important  addition  in  view  of  the  fact  that  there 
are  links  in  the  DEB  that  have  vertical  separation  distances 
of  a few  thousand  feet  between  the  transmit  and  receive  antennas 
(e.g.,  Hohenstadt-Zugspitze  link)  and  that  typical  layer  heights 
on  such  links  fall  below  the  highest  antenna. 

A further  augmentation  of  the  model  has  been  achieved  by 
incorporating  layer  tilt;  this  addition  as  described  earlier 
can  have  great  impact  on  the  fading  character  of  signals  on 
some  links,  especially  the  longer  ones. 

The  link  model  is  set  up  and  described  in  Section  2.2.1. 

The  solution  of  the  problem,  including  calculation  of  direct 
path  and  multipath  quantities,  is  described  in  Section  2.2.2. 

The  manner  in  which  the  model  is  adapted  to  include  layer  tilt 
is  described  in  Section  2.2.3.  In  Section  2.2.4  we  show  how 
the  analytical  results  can  be  used  to  determine  bounds  on 
refractive  index  gradient,  i.e.,  we  determine  how  large  a lapse 
rate  of  this  all-important  parameter  is  required  in  order  that 
multipath  be  supported  on  a given  link. 

The  analysis  described  in  this  section  is  the  mathematical 
framework  for  the  CNR  software  that  has  been  used  to  generate 
all  of  the  fading  characteristics  of  the  channel  presented  in 
this  report  and  the  propagation  parameters  (ray  path  delays  and 
amplitudes)  used  as  input  values  to  the  modem  performance  soft- 
ware described  in  Section  3. 


2.2.1  Description  of  Link 


The  manner  in  which  refractive  layers  are  generated  by 
atmospheric  conditions  was  described  in  Section  2.1.2. 
Throughout  the  remainder  of  this  report,  it  is  assumed  that 
there  exists  a stable  layer  of  finite  thickness  at  the  link 
location.  The  bottom  interface  of  this  layer  can  lie  above 
or  below  the  highest  terminal,  as  illustrated  in  Figure  2.5. 

In  Figure  2.5(a)  we  see  the  situation  wherein  the  layer  inter- 
face is  above  the  highest  receiver;  this  is  subsequently 
referred  to  as  condition  "A".  It  will  be  recalled  that  this 
situation  has  been  analyzed  before  with  some  computations 
appearing  in  [2.3]  and  [2.13].  Because  of  limitations  in 
application  to  links  with  extreme  geometries  (like  some  of 
those  in  the  DEB),  it  was  deemed  necessary,  however,  to  per- 
form an  analysis  of  the  situation  shown  in  Figure  2.5(b). 

In  this  figure  we  see  the  lower  layer  interface  positioned 
below  the  receiver.  We  call  this  condition  "B".  The  dif- 
ference between  this  situation  and  the  previous  model  developed 
at  CNR  is  that  the  receiver  is  located  within  the  layer. 

Expanded  views  of  the  situations  illustrated  in  Figure 
2.5  are  shown  in  Figures  2.6  and  2.7.  The  first  of  these, 
which  appeared  in  an  earlier  report  [2.13],  illustrates 
condition  "A";  the  layer  interface  is  above  the  receiver. 

The  second  illustrates  condition  "B"  which  will  be  given  a 
detailed  treatment  in  Sections  2.2.2  and  2.2.3. 

The  primary  difficulty  in  the  analysis  (as  with  the  pre- 
vious analysis)  relates  to  determining  the  ray  transmit  and 
receive  angles,  0^  and  0q,  respectively.  The  nature  of  the 
medium  is  such  (constant  refractive  index  gradient)  that  once 
these  two  angles  are  determined,  the  single  corresponding 
multipath  ray  path  is  completely  defined;  the  associated  delay 
and  power  loss  can  then  be  quickly  computed.  A multiplicity 
of  ray  paths  is  possible  because  of  the  multiplicity  of  solu- 
tion pairs  0j  and  0j^  and  because  of  the  variable  number  n of 
L2  segments*.  Unfortunately,  the  fact  that  the  current  model 
(condition  "B")  places  the  receiver  in  a region  characterized 
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case  of  interest  is  n = 0.  For  any  given  link 
the  case  n > 0 requires  a more  extreme  refractive 
of  the  layer  than  the  n = 0 case. 
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Figure  2.6  Geometry  of  Single  Stratified  Layer  Above  Terminals,  Condition 
From  [2.13] 


by  a different  refractive  index  than  that  characterizing 
the  region  surrounding  the  transmitter  results  in  some  com- 
plications of  the  analysis.  We  have  been  able  to  solve  the 
problem  exactly,  however;  the  details  of  the  solution  are 
provided  in  Section  2.2.2. 

In  terms  of  the  parameters  and  geometries  illustrated  in 
Figure  2.7,  it  can  be  simply  stated  that  the  CNR  softv^are, 
based  on  analytical  results  presented  in  this  section,  com- 
putes the  amplitude  and  delay  for  each  of  the  ray  paths 
supported  by  the  propagation  channel.  It  is  in  terms  of 
these,  of  course,  that  complete  characterization  of  the 
communication  channel  is  achieved. 

Before  proceeding  with  the  analysis,  it  is  important  to 
separate  and  identify  the  input  parameters.  This  is  required 
because  of  the  wide  variety  of  different  results  obtainable 
with  different  input  parameters  and  the  importance  of  carefully 
selecting  the  input  parameters  in  accordance  with  available 
measurements.  The  list  of  input  parameters  is  provided  in 
Table  2-2;  we  have  categorized  the  parameters  according  to 
whether  they  are  link  design  parameters  (man-made)  or  layer- 
related  parameters.  We  have  also  categorized  them  according 
to  whether  they  are  dynamic  or  static. 

The  initial  analytical  treatment  assumes  a perfectly 
horizontal  refractive  layer;  the  tilted  layer  case  can  be 
converted  to  this  simpler  problem  via  the  coordinate  trans- 
formation discussed  in  Section  2.2.3.  Thus,  results  in  this 
section  will  not  depend  on  layer  slope  S. 

Additionally,  it  will  be  noted  that  there  are  no  fre- 
quency variables  in  Table  2-2.  Frequency  is  totally  a system 
parameter  (see  Section  2.3)  and,  even  though  all  the  system 
distortion  parameters  have  a strong  dependence  on  frequency, 
it  plays  no  role  in  the  determination  of  the  ray  path  delays 
and  amplitudes. 

2.2.2  Solution  of  Problem  Using  Geometrical  Optics 

The  equations  of  geometrical  optics  are  perfectly  suited 
for  use  in  determining  the  delays  and  amplitudes  characteri- 
zing LOS  links.  It  will  be  recalled  from  electromagnetic 
theory  that  wave  propagation  is  governed  by  the  equation  of 
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the  eikonal  (geometrical  optics  limit)  when  two  basic  conditions 
hold: 

(1)  The  index  of  refraction  does  not  change  appre- 
ciably in  a distance  equal  to  the  wavelengths 
of  radiation. 

(2)  In  a distance  equal  to  the  wavelength  of  radia- 
tion, the  fractional  change  in  the  spacing 
between  adjacent  rays  must  be  small  compared 

to  unity. 

Clearly,  geometrical  optics  becomes  more  applicable  as 
one  goes  to  higher  frequencies.  On  the  basis  of  condition  (1), 
a LOS  link  operating  at  a carrier  frequency  f ^ 1 GHz  was 

judged  to  be  safely  modelable  using  the  geometrical  optics 
approximation  [2.3].  On  this  basis  it  is  deemed  thoroughly 
applicable  to  the  4-  and  8-GHz  links  treated  in  this  report. 

It  is  not  as  easy  to  determine  whether  geometrical  optics 
applies  using  condition  (2).  Work  done  by  Wong  [2.14],  however, 
indicates  that  in  the  presence  of  atmospheric  layers  the  rays 
travel  in  well-defined  bundles.  Within  each  of  these  bundles, 
condition  (2)  is  satisfied  at  all  but  a finite  number  of  caustic 
points.  At  these  caustics,  which  are  formed  by  the  intersection 
of  adjacent  rays,  the  ray  theory  can  be  slightly  modified  [2.15], 
[2.16]  to  provide  approximate  field  intensity.  In  the  other 
regions  where  the  separate  ray  bundles  intersect,  the  relative 
amplitudes  and  phases  of  each  bundle  can  be  calculated  from 
ray  theory,  since  each  satisfied  condition  (2).  These  regions 
are  known  as  multipath  regions,  interference  regions,  and 
radio  anti-holes. 

The  use  of  ray  techniques  has  long  been  a commonly-accepted 
practice  for  analysis  of  microwave  line-of -sight  links  and  is 
used  throughout  this  report  without  further  discussion. 

Before  commencing  with  the  analysis,  it  is  instructive 
to  note  two  interesting  and  important  features  of  the  geometri- 
cal approach: 

(1)  The  time  it  takes  a wavefront  to  pass  from  point  a 
to  b along  a path  is  proportional  to  the  integral 
of  refractive  index  along  the  path  between  the 
points ; i . e . , 
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where  n is  the  refractive  index  and  Cq  is  the 
free  space  speed  of  light. 

(b)  When  the  refractive  index  of  a medium  changes 
linearly,  i.e.,  when  its  gradient  G = vn  is  a 
constant,  the  radius  of  curvature  of  each  ray 
path  in  the  medium  is  given  by: 
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Equation  (2.5),  representing  delay  as  a path  integral,  is  the 
basic  means  whereby  the  path  delays  on  LOS  links  are  calculated. 
One  need  only  know  the  refractive  index  and  the  path  geometry. 

If  G is  a constant,  the  value  of  n at  any  point  is  determined 
from  knowledge  of  hL  and  G in  Table  2-2.  Determination  of  the 
ray  path  geometry  is  a more  difficult  task,  however,  and  is  a 
main  focus  in  the  remainder  of  this  section.  We  simply  point 
out  here  that  (2.6)  is  a crucial  aid  in  determining  the  path 
geometry;  it  indicates  that  the  radius  of  curvature,  like  the 
gradient,  is  a constant.  Hence,  the  ray  paths  are  segments 
of  circles,  an  important  simplifying  feature. 

The  analytical  approach  to  this  problem  can  be  logically 
divided  into  two  efforts:  one  is  directed  at  determining 
delays  and  amplitudes  on  the  direct  path;  the  other  is  directed 
at  determining  delays  and  amplitudes  of  the  multipath  rays. 

The  distinction  between  the  direct  path  and  the  multipath 
rays  is  important.  The  direct  path  is  that  path  that  exists 
between  transmitter  and  receiver  when  there  is  no  refractive 
layer  present.  A normal,  or  "average",  atmosphere  is  assumed. 
This  path  can  alternatively  be  called  the  "free-space"  or 
"clear  atmosphere"  path,  but  we  make  no  use  of  such  nomen- 
clature. The  multipath  rays,  of  course,  are  those  rays  that 
are  refracted  by  the  atmospheric  layer.  Fortunately,  the 
direct  path  yields  simple  analytical  expressions  for  delay 
and  amplitude,  and  is  always  used  as  a reference  in  this 
report.  Thus,  rather  than  referring  to  a multipath  ray  with 
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delay  500,002  ns,  vis-a-vis  a direct  path  with  delay  500,000  ns, 
we  refer  to  the  multipath  ray  as  having  delay  2 ns.  A similar 
statement  can  be  made  regarding  the  amplitude.  Summarizing, 
we  merely  state  that  henceforth  in  this  report  the  path  ampli- 
tude a and  path  delay  i of  a multipath  ray  refer  to  the  relative 
quantities  defined  below.  For  the  ray  (out  of  a total  of 

N rays ) , 


Rel  abs  abs 

a . = a . /a, 

J J 1- 
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where  the  direct  path  is  taken,  as  always,  to  be  path  #1. 

Since  we  are  using  a and  4 without  the  superscripts  to  denote 
the  relative  quantities,  we  have  a,  = 1 and  = 0 for  the  direct 
path.  The  reader  will  note  that  this  convention  is  exactly  the 
same  as  considering  the  overall  channel  transfer  function  and 
normalizing  with  the  direct  path  channel  transfer  function. 


The  normalization  described  above  will  be  used  indepen- 
dently of  whether  or  not  a direct  path  between  transmitter  and 
receiver  exists.  That  both  cases  can  hold  should  be  obvious 
from  Figure  2.5.  In  part  (a)  of  the  figure,  a direct  path  to 
the  receiver  always  exists.  In  part  (b)  of  the  figure,  there 
is  no  direct  path.  In  either  case,  however,  all  quantities  will 
be  referenced  to  the  direct  path. 


All  of  the  calculations  in  this  report  make  use  of  the 
earth-flattening  coordinate  transformation  described  in  2,13, 
p.2-14].  From  the  cited  reference  we  have  taken  Figure  2.8 
which  summarizes  the  effects  of  the  transformation.  From  the 
figure  it  will  be  noted  that  a ray  whose  curvature  exactly 
matches  that  of  the  earth  has  an  infinite  radius  of  curvature 
in  the  f lattened-earth  coordinate  system,  whereas  a perfectly 
straight  ray  appears  as  a circular  segment  with  concave  side 
upwards  in  the  flattened  earth  coordinate  system. 


For  the  direct  path,  the  atmosphere  is  assumed  to  have 
its  standard  value  G = -39  Nu/km.  In  the  flattened  earth  coor- 
dinate system,  this  corresponds  to  a concave -upward  circular 
segment  with  radius  of  magnitude  R^ = 5280  miles.  Generally, 
the  radius  of  curvature  is  related  to  the  gradient  via  the 
relation 
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R = actual  radius  of  curvature 
^mod  ~ modified  radius  of  curvature 
a = radius  of  earth 


ray  R=a 

Earth 


, ray  R , = “ 

n-profile  mod 

m-prof ile 

mil  II I III  II 1 1 1 

Flattened  Earth 


ray  R = oo 


Earth 


n-prof ile 
m-prof ile 


I I I I I I I 


Flattened  Earth 


m = n+z/a 
M = N+z/a 


Figure  2.8  Effect  of  Earth-Flattening  Coordinate  Transformation 
of  Ray  Trajectories  (From  [2.13]) 
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R 


(2.9) 


- 0.62139  X 10^ 
G + 156.78 


It  will  be  noted  that  this  expression  yields  an  infinite  radius 
and  a change  of  sign  at  the  transition  point: 


G = -156.78 


(2.10) 


It  is  convenient  to  describe  this  behavior  in  terms  of  the 
lapse  rate  of  refractive  index  gradient,  which  is  just  the 
algebraic  negative  of  G.  Were  the  lapse  rate  to  rise  from 
its  "standard-atmosphere"  value  of  39  Nu/km  to  some  large  but 
realizable  value  (e.g.,  300  Nu/km),  the  ray  path  would  first 
be  concave -upward , reach  a point  at  approximately  157  Nu/km 
where  it  becomes  perfectly  flat,  and  then  become  concave- 
downward  bending  the  rays  back  toward  the  earth. 

The  earth- flattening  coordinate  transformation  has  the 
effect  of  modifying  the  refractive  index  n,  resulting  in  the 
"modified  refractive  index"  m.  The  relationship  between  the 
two  quantities  is  given  in  Figure  2.8,  and  is  discussed  at 
length  in  [2.13,  Section  2. 1.1. 2]. 

2. 2. 2.1  Calculation  of  Direct  Path  Quantities 

Like  the  refractive  index,  the  modified  refractive  index 
varies  linearly  with  height  and  can  be  simply  represented  as 


m(z)  = m^  + gQ  z 


(2.11) 


where  z is  the  height  above  some  conveniently  chosen  zero  point, 
mQ  is  the  value  of  the  modified  refractive  index  at  that  point, 
and  gQ  is  the  direct  path  gradient  of  the  modified  refractive 
index . The  radius  of  curvature  is  given  by 
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Thus,  from  our  standard  atmosphere  assumption  for  the  direct 
path  (Gq  = "39  Nu/km) , gQ  and  Rq  are  immediately  determinable 
via  (2.9)  and  (2.12). 

The  direct  path  geometry  is  illustrated  in  Figure  2.9; 
the  ray  path  curvature  has  been  grossly  exaggerated  to  make 
it  visible.  From  the  geometry,  we  see  that  any  point  along 
the  ray  path  positioned  such  that  its  radius  arm  is  at  an 
angle  a with  the  vertical  is  at  a height  z given  by 


z = h^  - (R^  cos  a - R^  cos  a^)  (2.13) 

Substituting  in  (2.11), 

m = m^  + g^h^  + cos  - cos  a (2.14) 

The  time  delay  is  proportional  to  the  integral  of  m over  the 
path,  i.e.,  over  the  circular  arc.  The  differential  arc 
length  is  given  by 


dcj  = R^  d a (2 . 15) 

and  mg  + ggh-j-  « 1 

Performing  the  simple  integral  and  dividing  by  c,  the  speed 
of  light  in  vacuo,  we  obtain 


(sin  cxj,  + sin  o^^) 


(1  + cos  a^)(a^  + o^) 


(2.16) 


A companion  expression  for  (2.16)  represents  power  on  the 
path  as  a function  of  the  angles  otj,  and  This  can  be 

derived  by  following  exactly  the  same  procedure  described  in 
[2.3]  and  [2.13]  but  omitting  the  simplifying  assumption 
that  sets  cos  Ofrp=  1.  As  a result,  we  get  the  power  factor* 


The  square  root  of  this  quantity  is  the  ray  path  amplitude 
that  we  have  repeatedly  referred  to. 


2-28 


Figure  2.9  Direct  Path  Geometry  for  LOS  Link 

(Rj--Ro) 
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^ (sin  "tp' 


(2.17) 


9 = 


1 

L sin  a 

R 


The  problem  of  calculating  otj,  and  remains.  This  is 
a circumstance  that  will  characterize  the  multipath  analysis 
performed  later,  and  has  been  a standard  approach  of  ours  in 
solving  problems  of  this  type:  first  solve  for  the  amplitudes 
and  delays  in  terms  of  the  transmit  and  receive  angles;  then 
solve  for  those  angles. 

This  last  step  is  relatively  simple  here.  We  note,  from 
the  geometry  in  Figure  2.9,  that  the  following  relations  hold: 


sin  = - 


sin 


cos 


Of, 


R 


cos  (Xj, 


(2.18) 


Squaring  both  sides  of  each  of  the  equations  in  (2.18), 
adding,  and  then  solving  the  resulting  quadratic  equation  for 
sin  gives  sin  Oij,  = Sj  and  sin  where  S.|,  and 

quantities  worth  defining  separately,  are  given  by 


- 1 

(2.19) 


where  r is  the  slant  range, 

r = + (^)^  (2.20) 
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and  ^ is  the  vertical  terminal  difference. 


(2.21) 


In  terms  of  these  quantities,  we  have  the  direct  path  delay 
and  amplitude  given,  respectively,  by 


l!i 

dL 


(2.22) 


The  amplitude  can  be  determined  directly  in  terms  of  the 
geometrical  parameters,  h-j.,  h^,  etc.  Substituting  the  ex- ^ 
pression  for  and  into  the  expression  for  we  obtain, 
after  some  manipulation,  the  result 

■ li 


a,  = — 

1 r 


1 + 


I A ^ 

[ZK^jhehy 


1 - 


,2R, 


' r 
i2R 


f/J 


(2.23) 


This  is  quite  reasonably  approximated  by  the  expression, 


2 

1 

V 

L ^ 

r 

1 r ^ 

1 + 

L A 

(2.24) 


which  is  the  expression  we  have  used  in  our  programs. 

It  should  be  pointed  out  that  (2.24)  takes  into  account 
larger  values  of  th  than  the  expression  used  to  generate  re- 
sults in  an  earlier  report  under  another  contract.  The 
earlier  approximation,  0£^®=>L~^,  which  was  appropriate  for 
earlier  applications,  is  inappropriate  for  the  extreme 
geometries  of  the  links  encountered  in  the  current  study. 
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2. 2. 2. 2 Calculation  of  Multipath  Quantities 


Since  the  direct  path  delays  and  amplitudes  are  directly 
calculable  from  the  expressions  just  derived,  we  now  take  the 
second  step  toward  determining  the  relative  multipath  delays 
and  amplitudes  in  (2.7)  and  (2.8);  we  now  calculate  the  delays 
and  amplitudes  of  the  rays  refracted  by  the  atmospheric  layer. 

As  discussed  before,  we  are  interested  in  both  of  the  situations 
illustrated  in  Figure  2,5.  Condition  "A"  in  the  figure,  which 
is  *lustrated  in  more  detail  in  Figure  2.6,  has  been  treated 
in  earlier  reports.  We  will  treat  it  briefly  at  a later  point 
in  this  section,  but  our  treatment  will  be  limited  to  presenting 
(not  rederiving)  the  earlier  results  in  the  formalism  adopted 
here  and  to  a few  comments  on  the  advantages  of  looking  at  the 
problem  in  this  new  way.  Of  particular  interest  in  this  regard 
is  an  inquality  that  we  generate  in  Section  2.3.3  which  demon- 
strates the  limitations  of  the  condition  "A"  results  when 
trying  to  apply  them  to  links  with  a large  vertical  separation 
between  transmitter  and  receiver,  e.g.,  the  Hohenstadt-Zugspitze 
link  in  Germany. 

Our  primary  interest  here  is  in  condition  "B"  of  Figure 
2.5  which  is  illustrated  in  much  more  detail  in  Figure  5.7. 

This  is  the  case  wherein  one  of  the  terminals  is  in  the  layer, 
a likely  situation  if  there  is  a significant  difference  in  the 
terminal  heights. 


To  analyze  this  case,  we  refer  to  the  geometry  and  nomen- 
clature of  Figure  2.7.  The  problem  is  primarily  one  of 
analyzing  the  system  geometry  to  determine  path  delays  and 
power  loss  in  terms  of  the  transmit  and  receive  angles  9q  and 
or,  equivalently,  in  terms  of  the  normalized  parameters 
kj  and  kj^,  which  are  defined  according  to 


sin  0 


R 

0 


2R 


(2.25) 


sin 


(2.26) 


where  L is  the  ground  range  between  transmitter  and  receiver, 

R is  the  radius  of  curvature  of  the  ray  in  the  layer,  and  R£ 
is  the  radius  of  curvature  of  the  ray  in  free  space.  The  radii 
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are  related  to  the  modified  refractive  index  gradients  gQ  of 
the  free  space  below  the  layer  and  g^^  of  the  layer  through 


1 

and 

1 

^1  ■ R 


(2.27) 


(2.28) 


They  are  related  to  the  unmodified  refractive  index  gradients 
through  the  earlier  expression  in  Figure  2.8.  In  terms  of  gQ 
and  g]^  and  the  reference  level  refractive  index,  mQ,  it  is 
helpful  to  define  a free  space  index: 


Mp  = ">0  + gQ^ 


(2.29) 


and  an  index  applicable  to  the  layer  given  by, 


Ml  “ '"o 


(2.30) 


We  now  turn  to  the  geometry  of  Figure  2.7.  A simple 
examination  gives  the  following  results: 


T 

= R^  sin  6^  - Rj  sin  9 

(2.31) 

= 2(R  + R^)  sin 

(2.32) 

= R sin  9j^  + R sin  0Q 

(2.33) 

The  section  L2  can  be  repeated  n times,  resulting  in  a total 
ground  range , 


L = + n L2  + 


(2.34) 
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Keeping  in  mind  our  goal  of  expressing  all  quantities 
in  terms  of  0^  and  Gq,  it  is  easy  to  determine  from  (2.31) 
and  (2.34)  that 


L - R sin  6^  + R sin  6^ 

sin  0,  = — ^ 

1 (2n  + l)(R  + R^) 


(2.35) 


Examination  of  Figure  2 . 7 also  yields  two  relations  that  will 
be  essential  in  the  later  development.  These  are: 

‘Y  " ■*'  ^f  ■ ^f  ®1  (2.36) 

and 

lY  = h^  - (r  cos  Sq  - R cos  0^)  (2.37) 


We  recall  that  the  delay  for  each  multipath  ray  is  found  simply 
by  performing  the  integral  of  refractive  index  along  the  path. 
The  calculation  is  facilitated  by  noting  the  useful  relations 
listed  below  which  give  the  variable  height  z of  any  point 
along  the  path. 

• In  Region  1: 

z = h^  + R^  cos  0Q  - R^  cos  0^  (2.38) 

• In  Region  2*  : 

z = hj^  + R(cos  0 - cos  0^)  (2.39) 

% 

• In  Region  2": 

z = h^  + R^(cos  0Q  - cos  0)  (2.40) 
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In  Region  3; 


z = 


R 


- (R  cos  0Q  - R c 


os 


e) 


(2.41) 


Substituting  these  expression  in  (2.29)  and  (2.30)  and 
then  performing  the  angular  path  integrations  of  refractive 
index  over  each  of  the  ray  segments  in  Figure  2.5  and  dividing 
by  c , we  obtain 


+ (2n+  I) 


(R+Rf) 


(MpoSi-sm 


(2.42) 


where 


is  defined  as, 


(2.43) 


Thus,  from  (2.35),  (2.42),  and  (2.43),  we  see  that  the  delay 
to  be  associated  wit|^  any  multipath  ray  depends  on  the  transmit 
and  receive  angles  0q  and  0q,  respectively. 

T R 

Having  determined  delay  as  a function  of  0q  and  0q,  we 
now  turn  to  the  calculation  of  power.  This  result  is  obtained 
by  considering  the  fundamental  intensity  law  for  ray  optics 
as  described  earlier.  The  result  is  given  by 


cos  0 


9 = 


0 


d 0 


0 


L sin  0 


R dL 
0 


(2.44) 


where  9 is  the  power  factor,  i.e.,  the  ratio  of  received  to 
transmitted  power. 

This  completes  the  representation  of  delay  and  power 
factors  as  functions  of  transmit  and  receive  angles.  Before 
writing  these  results  in  terms  of  the  more  useful  dimensionless 
parameters  in  (2.25)  and  (2.26),  we  note  that  gQh.j,  in  (2.43)  is 
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is  an  extremely  small  quantity.  In  fact,  examination  of  (2.42) 
and  (2.43)  reveals  that  we  can  replace  mQ  + gQh,p  with  unity 
without  incurring  error.  Thus,  using 


MpQ  = 1 + cos 


0 


(2.45) 


and  (2.35)  in  (2.42),  and  (2.25)  and  (2.26)  in  (2.42)  and 
(2.44),  we  obtain 


for  the  power  factor,  and 


(2.46) 


^f)„.  -l/'V 

TT- 


+ (2n+  1) 


for  the  delay.  In  this  expression. 


k.  + y 


C2n+ 1)(R  + Rj)  ( ■ 2 


(2.48) 


is  just  sin  in  (2.35). 

We  have  now  arrived  at  the  same  point  in  the  multipath 
analysis  that  we  were  at  in  the  direct  path  analysis  after 
(2.16)  and  (2.17).  It  remains  to  determine  the  transmit  and 
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and  receive  angles  (equivalent  to  determining  Icp  and  k„)  . This 
was  quite  simple  for  the  direct  path,  but  is  more  complicated 
here.  As  in  the  earlier  case,  however,  we  start  by  appealing 
to  the  geometry. 

From  (2.36)  and  (2.37)  we  have  the  constraining  equations. 


h^  = h^  - R^(cos  0^  - cos  6q)  (2.49) 

h^  = h^  + R (cos  0^  - cos  0Q  j (2.50) 

T R 

which  must  be  solved  simultaneously  for  0q  and  0q,  i.e.,  for 
k.j,  and  k^.  The  solution  is  facilitated  by  using  the  expansion 

cos  0j^-cos  0Q  = -^(sin^  0^  - sin^  0^  ^ + -I-  (sin^  0^^  - sin^  0^)  + ... 

(2.51) 

In  this  expansion  we  note  that  the  ratio  of  the  magnitudes 
(last  term  to  first)  is  given  by 


(2.52) 


which  will  be  entirely  negligible  for  the  small  angles  occurring 
on  LOS  links.  Consequently,  we  can  drop  all  terms  but  the 
first,  getting 

\ \ ~ "T  ( ®l)  (2.53) 

tY  = h^  + |(  sin^  0^  - sin^  8j  (2.54) 

Substituting  from  (2.25)  and  (2.26),  and  using  (2.35),  we 
obtain,  after  some  manipulation,  the  two  equations 
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1 


- 4 - A 

k^p2  - (k-1,2  .S 

n 


where 


R, 


P = 


r \ 


+ k. 


k = 


R 


'] 

I 8r 

- (h^-h^)  + p(h_^-Y)j-^ 


8r^ 

cy-V  -"%-‘y>|-t 

Li 

8r^ 


and 


^ (2n+l)^  (R  + Rf)^ 


(2.55) 

(2.56) 

(2.57) 

(2.58) 

(2.59) 

(2.60) 

(2.61) 


Tedious  manipulation  of  (2.55)  - (2.61)  indicates  that  the 
set  of  solutions  to  (2.55)  and  (2.56)  can  be  found  by  solving 
a quartic  equation  for  y and  then  using  either  (2.55)  or 
(2.56)  to  determine  x.  The  quartic  equation  is  given  by 


where 


4 3 2 

y + ay  + by  + cy  + d = 0 


4(l  + p^G)  / (P^G^  - 1) 


(2.62) 


b = 


[4G(p^-l)  - 2p^G^(P+2)  + (P-Q)]/[g(P^G^  - 1)] 

= 4[(Q-P)  - P^G(P+2)]/[g(P^G^  - 1)]  ^ ‘ ^ 

= [p^G(P+2)^  + A(P-Q)]/[g(p^G^  - 1)] 

2-38 


and 


P 


Q 

G 


B (A  + S) 
n 

4 


B (A  - S) 
n 

4 


A 

2 


2 

P (B  - 
n 

2p 


1) 

2“ 


1 


(2.64) 


As  it  happens,  quartic  equations  are  the  highest-order 
equations  for  which  general  solution  formulas  exist.  The 
existence  of  these  formulas  indicate  that  the  problem  is,  in 
essence,  completely  solved  by  (2,57)  - (2.63),  and  does  not 
require  recourse  to  iterative  numerical  procedures.  Addi- 
tionally, it  provides  us  with  a means  (by  putting  constraints 
on  the  coefficients)  of  determining  inequalities  which  govern 
the  existence  of  multipath  on  any  given  link.  This  topic  is 
discussed  in  Section  2.2.4. 


It  is  of  interest  to  examine  the  dependence  of  (2.61)  on 
n.  Recall  from  Figure  2.7  that  n can  be  thought  of  roughly 
as  the  number  of  times  the  ray  reenters  the  layer  after  its 
initial  entry.  For  any  given  link,  some  consideration  of 
Figure  2.7  reveals  that  increases  in  n require  greater  bending 
properties  of  the  layer,  i.e.,  larger  lapse  rates  of  the  re- 
fractive index  gradient.  From  our  earlier  discussion  in 
Section  2.1,  we  know  that  the  larger  lapse  rates  are  less 
likely  to  occur  in  nature.  In  fact,  it  has  been  our  experience 
in  examining  several  links  that  large  lapse  rates  are  required 
in  order  to  produce  multipath  in  the  minimal  reentry  case 
when  n = 0.  Thus,  we  are  already  at  the  low  end  of  the  proba- 
bility scale  when  we  consider  the  n = 0 case  for  practically 
all  links  of  interest. 


It  should  be  pointed  out  that  there  are  some  situations, 
themselves  improbable,  where  n>0  could  be  of  importance. 

This  situation  arises  when  the  layer  is  so  configured  that 
the  receiving  and  transmitting  terminals  are  both  very  close 
to  the  layer  interface.  A limiting  case  occurs  when  the  layer 
interface  coincides  exactly  with  receiver  and  transmitter. 

In  this  case,  the  ray  can  skip  in  and  out  of  the  layer  an 
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Formulas  for  solving  (2.62)  and  (2,65)  are  given  in  [2.  ]. 

Note  that  the  existence  of  the  two  simultaneous  equations, 
(2.55)  and  (2.56),  which  relate  and  k^,  can  now  be  used  to 
derive  an  analytical  expression  for  the  power  factor  9 in 
(2.46)  by  allowing  evaluation  of  the  derivative.  Manipulating 
(2.55)  and  (2.56),  we  obtain 


(2.67) 
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r 


where 


(2.68) 


and 


k 


T R 


(2.69) 


The  analytical  development  is  now  complete.  The  transmit 
and  receive  angles  can  be  determined  by  solving  the  quartic 
(2.62)  or  the  cubic  equation  (2.65);  and  by  substituting  the 
solution  back  into  the  constraining  simultaneous  equations 
(2.55)  and  (2.56).  With  a set  of  values  so  determined,  the 
amplitude  and  the  delay  for  the  corresponding  path  can  then 
be  calculated  from  (2.46)  and  (2.47). 

2.2.3  Inclusion  of  Layer  Tilt  into  Model 

The  analytical  results  presented  in  the  previous  sections 
correspond  to  the  physical  situation  wherein  the  refractive 
layer  is  perfectly  horizontal  with  respect  to  the  ground. 

This  is  not  always  the  case.  That  layers  are  subject  to 
tilting  has  been  reported  in  the  literature.  The  experimental 
data  is  limited,  but  these  limitations  are  not  due  to  absence 
of  the  phenomena;  they  exist  primarily  because  of  the  diffi- 
culty in  measuring  layer  tilt.  It  is  clear  that,  if  measure- 
ments of  layer  tilt  with  concurrent  measurements  of  refractive 
index  gradient  are  to  be  made,  one  must  have  the  services  of 
at  least  three  neighboring  radiosonde  stations  working  in 
close  conjunction  with  one  another.  Because  of  this,  data  on 
layer  tilt  is  extremely  sparse.  Gossard  [2.17]  has  presented 
some  limited  estimates  that  derive  from  his  study  of  refractive 
layers  perturbed  by  waves.  Because  of  the  indirect  manner  in 
which  these  estimates  were  derived,  they  have  not  been  deemed 
suitable  for  inclusion  here.  If  one  is  willing  to  make  the 
assumption  that  the  refractive  layers  derive  entirely  from 
frontal  inversions*,  one  can  use  the  data  of  Kerr  [2.18]  who 


Frontal  inversions  are  known  to  be  a causative  mechanism, 

but  are  not  thought  to  be  the  sole  mechanism  (see  Section  2.1.2). 
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has  pointed  out  that  wave  fronts  have  tilt  slopes  ranging 
between  0.0033  and  0.010,  corresponding  to  an  angular  range 
of  0.190  to  0.573  degrees.  For  lack  of  more  definitive  data, 
this  is  the  data  that  we  have  used  in  our  calculation  when 
deemed  appropriate . 

Layer  tilt  is  less  important  on  the  shorter  links  than 
on  the  longer  links.  Recollection  of  results  (condition  "A” 
only)  presented  in  an  earlier  report  [2.13]  and  some  simple 
geometrical  reasoning  indicates  why  this  is  so.  Earlier 
results  indicated  that  multipath  was  only  important  when  the 
layer  was  in  close  proximity  with  both  terminals.  When  the 
layer  moved  away  from  one  terminal  or  the  other,  the  multipath 
spread  would  decrease.  For  our  purposes  here,  we  simply  note 
that  the  multipath  condition  depends  sensitively  on  the  posi- 
tion of  the  layer  relative  to  both  terminals.  Tilt  of  a 
refractive  layer  over  a short  link  does  little  to  change  these 
relative  positions.  Over  a long  link,  however,  the  effect  can 
be  dramatic.  For  more  details,  the  reader  is  referred  to  our 
coverage  of  the  100-mile  Hohenstadt-Zugspitze  link  discussed 
in  Section  2.4.1. 

The  basic  work  on  determination  of  delays  and  amplitudes 
has  already  been  done,  and  the  phenomenon  of  layer  tilt  has 
been  incorporated  into  our  model  via  some  straightforward 
coordinate  transformations.  For  the  most  part,  these  coordi- 
nate transformations  reflect  the  changes  in  relative  posi- 
tioning of  the  terminals  and  layer  interface  due  to  layer  tilt, 
as  discussed  in  the  previous  paragraph. 

The  tilted  layer  configuration  is  illustrated  in 
Figure  2.10.  Were  the  tilted  planar  interface  to  stretch 
to  the  ground,  the  line  of  intersection  between  it  and  the 
ground  plane  would  be  perpendicular  to  the  vertical  plane 
containing  the  transmitter  and  receiver.  (Consideration  of 
a more  complicated  geometry  at  this  time  would  only  confuse 
the  issues.)  The  layer  is  seen  to  tilt  at  an  angle  0 relative 
to  the  earth.  It  is  reasonable  to  assume  that  the  tilt  applies 
to  the  refractive  index  profile  just  below  the  layer  interface 
as  well  as  i^  does  to  the  layer  proper.  As  a consequence  of 
this,  one  is  able  to  effect  a coordinate  transformation  that 
ultimately  enables  computation  of  the  multipath  profiles  for 
various  values  of  layer  height  hj^,  refractive  index  gradient  G, 
and  layer  tilt  slope  S.  As  inputs  to  the  called  subroutines, 
however,  one  uses  the  transformed  quantities 


w 
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Figure  2.10  Geometry  of  Layer  Tilt 
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0 


where 


cos 


h ' = r sin  JS 
R s 


L'  = r cos  P 
s 


s > R 


and 


^ = tan 


(2.70) 


(2.71) 


It  is  interesting  to  note  that  fading  on  some  links  (e.g., 
the  Hohenstadt-Zugspitze  link)  would  be  unexplainable  via  the 
mechanism  of  refractive  multipath  in  atmospheric  layers  were 
it  not  for  the  fact  that  layers  sometimes  tilt.  As  it  happens, 
the  layer  tilt  angle  is  a variable  of  major  importance  for  links 
with  extreme  geometries. 

2.2.4  Bound  on  Refractive  Index  Gradient 

To  determine  the  appropriateness  of  utilizing  the  condi- 
tion "A"  results  derived  under  an  earlier  contract  [2.13],  it 
was  decided  early  on  in  this  program  to  find  out  if  the  geometry 
of  longer  links  (with  a large  vertical  separation  between 
transmitter  and  receiver  terminals)  implied  a requirement  of 
extremely  large  lapse  rate  in  order  for  multipath  to  be  sup- 
ported on  the  link.  As  we  demonstrate  below,  this  turned  out 
to  be  the  case,  a fact  that  was  the  prime  motive  for  developing 
a set  of  equations  that  governed  condition  "B”.  Because  of 
time  limitations,  our  investigation  was  limited  to  the  case 
wherein  the  layer  is  horizontally  situated  relative  to  the 
earth . 

The  current  development  has  the  advantage  of  providing 
a method  of  approaching  the  governing  simultaneous  equations 
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in  such  a way  as  to  determine  the  realness  (or  complexness) 
of  the  roots  to  the  system.  Realness  of  roots  is  required, 
of  course,  in  order  to  have  multipath.  The  conditions 
governing  the  realness  of  the  roots  are  identical  to  the 
conditions  which  indicate  whether  or  not  multipath  exists  on 
the  link.  It  is  felt  that  a powerful  approach  is  available 
here  for  determining  the  ranges  of  parameter  values  which  allow 
multipath  on  the  link.  Unfortunately,  we  have  not  been  able  to 
pursue  this  question  in  depth.  The  development  here  is  limited 
to  determining  the  bounds  on  lapse  rate  under  condition  "A"; 
this  course  was  pursued  for  the  reasons  stated  above.  One 
germane  conclusion  here  is  that  lapse  rates  exceeding  those 
occurring  in  nature  are  required  if  multipath  (under  condition 
"A")  is  to  be  supported  on  some  of  the  "extreme”  geometry  links, 
e.g.,  the  Hohenstadt-Zugspitze  link. 

Our  development  is  initiated  by  applying  a change  of 
variables  (that  preserves  the  realness  or  complexness  of  the 
roots)  to  the  set  of  simultaneous  quadratic  equations  governing 
propagation  under  condition  "A".  We  apply  the  change  of 
variables 


X -f  y 


X - y 
^ ~ 2 


to  the  system  of  equations  (see  [2.13,  p.C-1]) 


to  get 


2 2 ^ 
X - y - A 


n 


uv 


(u  + v)^  - ^ (u  - 1)^  + (u  - v)^  - -S 
“n 


(2.72) 


(2.73) 

(2.74) 


(2.75) 

(2.76) 
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In  these  expressions  we  have 
definitions : 

“n  - 

A = 8 (hR  - h^)/L^ 

S = 16  (h  - h)/L^ 


sed  the  condition  "A” 

(2.77) 

(2.78) 

(2.79) 


where  Rf  is  the  radius  of  curvature  below  the  layer,  h is  the 
layer  height,  and  h is  the  arithmetic  average  of  receiver 
height  hj^  and  transmitter  height  h^. 

We  are  interested  in  the  conditions  that  must  be  met  in 
order  that  there  be  a pair  or  real  roots,  x and  y,  i.e.,  a pair 
of  real  roots,  u and  v.  Assuming  a real  solution  for  v,  we 
must  also  have  a real  solution  for  u.  Evaluating  the  discri- 
minant of  (2.76),  we  obtain  the  inequality 


n 

With  the  aid  of  the  trivial  inequalities 


2 

V 2 0 

B 2 1 

n 


we  get 


B 

n 


which. 


after  more  manipulation,  yields 


B 

n 


1 


(2.80) 


(2.81) 


(2.82) 


(2.83) 
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Substituting  (2.77)  and  (2.79)  into  this  expression,  some 
further  manipulation  gives 


(2.84) 


Thus,  in  terms  of  the  geometrical  parameters  of  the  link  (L  and 
h)  and  the  physical  parameters  of  the  direct  path  atmosphere 
(R£)  and  the  layer  (R  and  h) , we  have  determined  a bound  on 
n,  the  number  of  excess  passages  through  the  layer  (see 
Figure  2.6). 


In  order  for  multipath  to  exist  on  the  link  at  all,  the 
right-hand  side  of  (2.84)  must  be  equal  to  or  exceed  zero.  We 
get,  after  some  manipulation,  the  inequality 


Q(R)  ^ 0 (2.85) 

where 

Q(R)  = (16A)R  - (2.86) 

and  A=  h-h.  This  expression  increases  monotonically  with  R. 
Thus,  for  a positive  root  Rq.  of  the  quadratic  equation  Q(R)  =0, 
the  inequality 


R ^ R^ 


(2.87) 


must  hold  if  the  inequality  (2.85)  is  to  hold  (which  it  must 
if  multipath  is  to  exist  on  the  link);  the  inequalities  (2.83), 
(2.84),  (2.85),  and  (2.87)  are  all  inequalities  that  must  be 
satisfied  for  multipath  to  exist  on  the  link.  Thus,  (2.87) 
provides  us  with  the  upper  bound  on  R,  the  radius  of  curvature 
of  the  ray  in  the  layer.  Solving  for  the  positive  root  of 
(2.86),  wie  get 
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1 


R s:  R 


max 


(2.88) 


where 


R = R^ 
max  f 


8 R^  A 


- 1 


(2.89) 


The  ray  cannot  curve  in  an  arc  to  less  a degree  than  indicated 
by  (2.88)  and  (2.89)  if  multipath  is  to  exist  on  the  link.  Thus, 
the  bending  requirement  of  the  layer,  if  it  is  to  support  multi- 
path,  is  firmly  established.  The  reader  will  note  that  this 
puts  a bound  on  refractive  index  gradient;  i.e.,  the  lapse  rate 
of  the  layer  must  exceed  some  value  corresponding  to  (2.89)  if 
multipath  is  to  exist  on  the  link.  The  relationship  betwen  R 
and  the  lapse  rate  af  is  given  by 


R = 


.62139  X 10 
156.78 


(2.90) 


where  R is  in  miles  and  is  in  Nu/km.  When  =£  is  in  excess  of 
156.78  Nu/km,  it  is  clear  from  this  relation  that  R decreases 
with  increasing  lapse  rate.  From  (2.88)  we  see  that  R obtains 
a value  small  enough  to  support  multipath  on  the  link  when  the 
inequality 


2 


.62139  X 10 


R. 


yrr 


8 Rj^  A 


+ 156.78 


(2.91) 


is  satisfied. 


Curves  illustrating  the  lapse  rates  required  to  allow 
multipath  on  different  links  are  shown  in  Figure  2.11.  We  see 
that  the  lapse  rate  requirement  rises  sharply  as  A increases, 
i.e.,  as  the  layer  rises. 

The  curves  in  Figure  2.11  have  important  implications 
relative  to  the  limitations  of  the  condition  "A”  model.  Recall 
that  the  model  requires  that  the  layer  be  above  the  highest 
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terminal  (chosen  without  loss  of  generality  to  be  the  receiver). 
Hence , 


A = h 


h + h„  + h„ 

' ‘'r  ■ ““2— 


(2.92) 


Thus,  the  lapse  rate  requirement  for  any  given  link  is  seen  to 
increase  as  - h.j.  increases.  In  fact,  the  limitations  of  the 
condition  "A”  model  can  be  determined  by  replacing  A with 
(hR-hT)/2  in  Figure  2.11,  e.g.,  for  hj^  - h^  = 200  on  a 200-mile 
link,  we  see  that  a lapse  rate  of  680  Nu/km  is  required  in 
order  for  the  condition  "A"  LOS  multipath  model  to  give  rise 
to  multipath.  In  view  of  data  in  [2.  ] and  [2.  ],  this  is 

an  extreme  requirement. 

For  the  100-mile  Hohenstadt-Zugspitze  link  in  Europe,  the 
requirement  is  even  more  extreme.  Examination  of  the  link 
geometry  gives 


= 3511.16  feet 


(2.93) 


which,  when  used  in  (2.91),  indicates  that  a lapse  rate  in 
excess  of  872  Nu/km  is  required  for  the  earlier  model  to 
produce  multipath  on  the  link.  Clearly,  the  condition  "A" 
multipath  model  is  unsuitable  where  such  extreme  vertical 
separations  of  transmitter  and  receiver  are  encountered. 

On  a worldwide  basis,  the  measurement  of  refractive  index 
lapse  rate  has  been  going  on  for  several  years.  As  far  as 
we  have  been  able  to  determine  from  our  study  of  the  literature, 
the  inequality  in  (2.91)  constitutes  the  first  and  only 
"existence”  condition  (in  terms  of  refractive  index  lapse 
rate)  for  multipath  fading  on  LOS  links.  Since  it  is  not 
one  of  the  stated  goals  of  this  contract,  time  limitations 
have  prevented  us  from  attempting  to  determine  the  "existence" 
condition  for  condition  "B"  of  Figure  2.5.  It  is  expected 
that  the  analysis  would  be  more  difficult,  but  would  proceed 
along  the  same  lines  developed  here. 


Other  "existence"  conditions  related  to  different  quantities, 
e.g.,  path  length,  layer  thickness,  have  been  presented  in 
earlier  CNR  reports  [2.13]. 
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It  is  of  interest  to  determine  the  lapse  rate  requirement 
V7hen  n,  the  number  of  excess  passages  through  the  layer, 
exceeds  zero.  Such  a condition  corresponds  exactly  with  the 
condition  that  the  total  number  of  ray  paths  is  greater  than 
three,  as  discussed  in  connection  with  (2.62)  through  (2.65). 
An  analysis  similar  to  that  described  above  indicates  that  the 
lapse  rate  required  to  result  in  n excess  passages  is  given  by 


> 


n 


,62139  X 10 


8R^A 


6 


(n  + 1) 


+ 156.78 


n > 1 


Since  the  right-hand  side  of  this  expression  must  always 
exceed  156.78  (if  multipath  is  to  exist  at  all)  we  see  that 
a larger  lapse  rate  is  required  if  extra  passages,  and  hence 
extra  ray  paths,  are  to  exist.  Since  large  lapse  rates  are 
invariably  less  probable  than  small  ones,  we  conclude  that  the 
existence  of  more  than  three  paths  is  less  probable  than  the 
existence  of  just  three.  In  all  but  a few  anomalous  situa- 
tions the  difference  in  probabilities  will  be  substantial. 
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2 . 3 System  Parameters 


In  this  section  we  categorize  and  describe  the  various 
system  parameters  of  importance  under  this  study.  All  of 
these  parameters  derive  from  the  time-varying  channel  transfer 
function  T(f,t)  and  serve  to  bridge  the  gap  between  the  physical 
parameters  describing  the  propagation  channel  and  the  quanti- 
tative measures  of  system  performance.  Theoretically,  one 
need  only  know  T(f,t)  to  evaluate  system  performance.  In 
practice,  one  (or  both)  of  two  conditions  obtain; 

(1)  Not  enough  information  about  the  propagation 
channel  is  available  to  precisely  reconstruct 
T(f,t); 

(2)  In  the  event  T(f,t)  is  available,  it  is 
difficult  to  use  it  in  its  unmodified  state. 

One  is  forced  to  consider  finite  parameter  characterizations 
of  the  channel.  The  whole  topic  of  models  for  random  time- 
varying  channels  (including  reductions  to  finite  parameter 
models)  has  been  discussed  comprehensively  by  Bello  [2.19]. 

In  this  report  we  represent  signal  degradation  in  terms  of 
the  T-parameters  and  P-parameters , a development  motivated 
primarily  by  condition  (2)  above.  The  T-parameter  represen- 
tation is  of  a general  nature,  whereas  the  P-parameter  charac- 
terization arises  because  of  the  emphasis  here  on  angle 
modulation . 

Two  of  the  parameters  mentioned  above,  because  they  have 
classically  played  a dominant  role  in  providing  a physical 
description  of  the  channel  and/or  in  providing  a means  whereby 
one  can  gauge  system  performance,  are  given  a separate  treat- 
ment in  Section  2.3.1.  These  are  Tq  and  P]^ ; the  magnitude  of 
Tq  is  known  classically  as  the  fade  level  (at  band  center)  and  the 
negative  of  P^  is  known  as  the  group  delay  (at  band  center). 

Because  characterization  of  frequency-selective  fading 
generally  requires  a greater  number  of  parameters  than  the  two 
mentioned  above.  Section  2.3.2  is  devoted  to  a general  dis- 
cussion of  the  T-  and  P-parameters.  Their  importance  relative 
to  the  distortion  of  angle-modulated  signals  is  emphasized. 

In  Section  2.3.3  we  return  to  a closer  look  at  fade  depth 
and  group  delay;  we  use  the  accepted  form  for  the  channel 
transfer  function  on  LOS  refractive  multipath  links  [see  (9)] 
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and  examine  the  joint  behavior  of  these  parameters  during  a 
deep  fade.  We  make  the  important  observation  that  the  largest 
excursions  in  group  delay  can  only  occur  during  a deep  fade  . . . 
an  observation  that  will  take  on  much  greater  importance  when 
we  later  consider  diversity  switching  in  Section  4,  Addition- 
ally, we  derive  an  expression  whereby  one  can  determine  those 
frequencies  at  which  the  deepest  fades  occur. 

2.3.1  Fade  Depth  and  Group  Delay 

A time-varying  communication  channel  can,  in  general,  be 
represented  by  a transfer  function  T(f,t)  which  provides  a joint 
depiction  of  the  frequency-selective  and  time-selective  behavior 
of  the  channel.  In  this  section  we  discuss  two  of  the  most 
important  parameters  associated  with  T(f,t).  These  are: 

(1)  the  fade  level  which  is  just  the  magnitude  of  T(f,t);  and 

(2)  the  group  delay  which  is  the  frequency  derivative  of  the 
phase  of  T(f,t).  The  quantity  !T(f,t)l,  when  evaluated  at  one 
particular  frequency  (e.g.,  the  carrier  frequency),  is  the 
commonly  referred-to  fade  record,  a function  of  time  giving 
rise  to  variations  in  the  received  signal  level.  When  evaluated 
at  one  particular  time  instant,  a "snapshot"  of  the  frequency 
transfer  function  is  obtained.  In  this  section  we  are  inter- 
ested in  characterizing  the  frequency-selective  behavior  of 

the  channel  and,  for  that  reason,  will  concentrate  on  the 
latter  point  of  view. 

For  our  purposes,  it  is  of  interest  to  consider  the  channel 
transfer  function  written  in  the  general  form, 

T(f,t)  = A(f,t)  (2.94) 


At  a particular  frequency  and  time,  we  refer  to  the  fade  level 
L(f,t)  = lT(f,t)l  (2.9 


which  is  just  the  channel  transfer  function  magnitude, 
this  is  given  by 

L^g(f,t)  = 20  log^Q'T(f,t) 


2-3  1 


b)  Equivalent  System 


Figure  3.3  Block  Diagram  Showing  Signal  Processing  Elements 
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(3.46) 


cp(t)  = <p'  + f (2i  - 3)  § 

As  the  conversion  between  reference  frames  is  always  a known 
phase  shift,  this  apparent  source  of  confusion  should  cause  no 
conceptual  difficulties.  Thus  we  let  0(t)  denote  the  (phase) 
response  of  the  BPF  to  the  transmitted  phase  <p(t) . The  relation- 
ship between  these  quantities  is 


0(t)  = tan 


h(t)  * sin  o(t) 
h(t)  * cos  (p(t) 


(3.47) 


where  h(t)  is  the  low  pass  equivalent  impulse  response  of  the  BPE 
(This  is  valid  under  the  assumption  that  h(t)  is  real,  corres- 
ponding to  a symmetric  bandpass  filter.)  The  (phase)  response  of 
the  channel  ■y(t)  may  in  turn  be  written  in  terms  of  9(t)  using 
the  p series  discussed  earlier: 


y(t)  = e(t)  + p^0(t)  + P2e(t)  + j P2e^(t)  (3.48) 


In  the  computer  evaluation  of  performance,  values  of  0(t) 
are  obtained  by  numerically  performing  the  convolutions  in  the 
numerator  and  denominator  of  Eq.  (3.47)  for  the  phase  waveform 
of  Eq.  (3.46).  Derivatives  of  0(t)  required  for  calculating  y(t) 
from  Eq.  (3.48)  are  evaluated  by  standard  numerical  procedures. 
That  is , 


e(t) 


6(t  +A)  - 0(t  - A) 
IL 


(3.49) 


e(t)  = eCt  + Al-..2i(t)  + 0.(t-j^  (350 

If  multipath  had  no  effect  on  the  extraction  of  reference  by  the 
receiver,  the  relationships  presented  so  far  would  be  sufficient 
to  calculate  modem  performance.  However,  this  is  not  the  case 
and  we  must  investigate  the  effect  of  multipath  on  both  the 
derived  phase  reference  and  the  derived  symbol  clock. 
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The  PLL’s  shown  in  Fig.  3.3  function  as  narrow  band  filters 
which  pass  only  a discrete  component  at  the  appropriate  fre- 
quency. Since  the  bandwidths  of  these  loops  are  usually  very 
much  narrower  than  the  system  bandwidth,  we  will  neglect  all 
noise  effects  and  model  them  as  having  zero-bandwidth.  Thus  we 
have  the  following  transfer  functions  for  the  zonal  filters  at 


SIq  and  8f^ 


HQ(f)  = 6(f  - SIq) 


H^(f)  = 6(f  - 8f^) 


(3.51) 


(3.52) 


Here  we  are  assuming  a center  frequency  8f^.  The  only  component 
in  the  output  of  the  eighth  law  device  which  will  have  any  effect 
is  also  that  at  81^,.  Neglecting  gain  factors,  this  component 
will  have  complex  envelope 


,(t)  - 


(3.53) 


Proceeding  in  a straightforward  manner,  the  respective  filter 
outputs  become 


and 


* j[8y(T)  - 2ff8fQT]  I j2ff8f 

J e dr  e 


» j[8r(T)  - 2ff8f  t]  I j2Tr8f  t 

f e drf  e 


(3.54) 


(3.55) 


Since  no  amplitude  effects  are  retained  by  the  PLL,  these  expres- 
sions may  be  further  simplified  to 


j2»8fQt  + 


(3.56) 
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A 


and 


j2ir8f^t  + 


(3.57) 


where  the  phase  shifts  are  given  by 


'ITq  = tan 


J sin  [SyCr)  - 2ff8fQT]  dr 
-1  

00 

J cos  [8y(T)  - 2ir8fQT]  dr 


(3.58) 


= tan 


•1  - 


J sin  [8y(T)  - InSf^r]  dr 


J cos  [8y(T)  - 2iT8f^T]  dr 


(3.59) 


Now  we  must  describe  how  the  four  integrals  appearing  in 
Eqs.  (3.58)  and  (3.59)  may  be  computed  numerically.  For  concrete- 
ness, we  consider  that  in  the  numerator  of  Eq.  (3.57),  although 
all  four  are  handled  in  an  identical  manner.  First,  we  break 
the  interval  of  integration  into  one  baud  intervals: 


/■kT  -t-  7 T„ 

® / 2 g 

r cos  [8y(t)  - 2r8f,t]  dt  = E / 

" kX-I+r 

2 g 


cos  [8y(t)  - 2ff8fQt]  dt 


(3.60) 


where  Tg  is  the  group  delay  of  the  receiver  bandpass  filter  H(co) . 
This  filter  is  designed  with  minimal  distortion;  in  particular, 

over  the  interval  (kT  ” ^ “rg*  kT  -I-  ^ -t-  Tg)  the  filter  output 

phase  waveform  will  depend  only  on  the  pair  of  frequencies  trans- 
mitted during  the  intervals  (kT-T,  kT)  and  (kT,kTH-T)  and  the 
transmitted  phase  at  time  kT.  It  can  easily  be  shown  that  all  of 
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the  eight  equally  spaced  phases  allowed  at  kT  give  the  same  value 
for  the  integrand  in  question.  Thus,  the  integral  in  Eq.  (3.60) 
will  take  on  exactly  16  distinct  values  corresponding  to  the  16 
possible  combinations  of  preceding  and  succeeding  frequencies. 

In  considering  Eq.  (3.60)  we  must  recall  that  the  summation  over 
an  infinite  number  of  T sec  intervals  is  a mathematical  fiction 
arising  from  the  fact  that  we  have  assumed  zero  bandwidth  for 
PLL.  In  fact,  the  summation  will  include  a finite  but  very  large 
number  of  T sec  integrals.  In  this  summation,  each  of  the  16 
distinct  values  for  the  T sec  integral  will  occur  with  equal 
probability,  assuming  a random  data  sequence.  The  number  of 
terms  in  the  summation  will  be  large  enough  that  we  may  neglect 
variations  of  the  number  of  each  distinct  type  about  its  expected 
value.  Thus  the  final  value  of  the  summation  will  be  propor- 
tional to  the  average  of  the  16  possible  values  for  the  basic  T 
sec  integral.  Thus,  in  calculating  Eq.  (3.60)  we  will  simply 
calculate  and  sum  the  16  distinct  values  taken  on  by  the  integral 
appearing  on  the  right-hand  side.  This  involves,  of  course,  cal- 
culating the  actual  value  of  y(t)  in  the  presence  of  multipath 
using  Eqs.  (3.46),  (3.47)  and  (3.48)  for  each  of  the  16  fre- 
quency combinations.  In  this  manner  the  Integrals  appearing  in 
Eqs.  (3.58)  and  (3.59)  are  numerically  evaluated  (to  within  a 
proportionality  constant)  and  ^'q  and  are  calculated. 

Having  obtained  'I'q  and  we  note  that  phase  error  in  the 
symbol  reference  is  'J/q  - so  that  the  actual  sampling  instant 
is  given  by 


t 

s 


'ko  - 'I', 


2it 


(3.61) 


Thus,  due  to  the  timing  error,  y(0)  in  Eqs.  (3.41)  and  (3.44) 
must  be  replaced  by  y(ts)  in  calculating  the  error  probability. 
However,  an  additional  term  arises  from  the  derived  phase 
reference: 


(3.62) 


where  'I'n  is  simply  Replacing  t by  tg  [from  Eq.  (3.61)] 

and  fQ  oy  its  value  - 3/16T  gives  the  phase  of  the  reference  at 
the  sampling  instant: 
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(3.63) 


‘^ref 


Now  we  may  rewrite  Eqs.  (3.41)  and  (3.44)  to  include  the  effect 
of  the  derived  phase  and  symbol  references: 


It 


•ft  + (k  - j - 1)6^  - y(t  ) + 


ref 


(3.64) 


(3.65) 


Actual  error  rates  are  calculated  from  (3.45)  as  before. 
3.2  Receiver  Performance  Over  Multipath  Channels 


3.2.1  Channel  Selection 

Section  3.1  above  has  described  the  analytical  basis  for  the 
computer  software  which  has  been  developed  to  calculate  perform- 
ance of  the  two  modems  of  interest  in  the  presence  of  multipath. 
In  this  section  we  present  results  obtained  with  this  software 
for  the  Hohenstadt-Zugspitze  link.  In  obtaining  this  data,  it 
was  desired  to  consider  a wide  range  of  multipath  conditions 
which  could  exist  on  this  link.  This  has  involved  two  processes 
of  selection  which  will  be  described  before  proceeding  to  the 
actual  performance  results.  The  first  selection  process  is  that 
of  selecting  a set  of  physical  channel  parameters  for  which  to 
calculate  the  path  structure.  Given  the  path  structure,  in  the 
form  of  amplitudes  and  relative  delays,  one  must  then  calculate 
the  parameters  which  characterize  the  communication  channel  at 
various  fade  depths . 

The  multipath  structure  of  the  Hohenstadt-Zugspitze  has  been 
discussed  in  detail  in  Section  2.4.  In  selecting  channel  condi- 
tions, we  have  taken  sample  points  from  each  of  the  multipath 
profiles  illustrated  in  Figs.  2.16  to  2.18.  The  layer  height 
samples  are  given  approximately  by 
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h-  = h . + n 

L min 


h - h . ' 
max  min 


n = 1,2,3 


(3.66) 


where  hmin  is  the  smallest  layer  height  resulting  in  multipath 
for  a given  layer  tilt  and  hmax  is  the  largest.  This  scheme  is 
illustrated  in  Fig.  3.4  where  we  have  numbered  the  cuts  according 
to  the  value  of  n in  (3.66).  This  numbering,  accompanied  by  the 
exact  layer  heights  involved,  also  appears  on  the  error  rate 
curves  to  be  presented  later  in  this  report.  For  the  refractive 
index  gradient  the  values  -400,  -350,  and  -300  Nu/km  are  used; 
these  are  in  reasonable  accord  with  relevant  data  [2.6].  For 
each  value  of  refractive  index  gradient  the  layer  slopes  given 
in  Table  3-1  have  been  used  in  the  calculations.  Not  all  combin- 
ations of  refractive  index  gradient  and  layer  slope  result  in 
multipath,  however.  All  of  the  multipath  profiles  resulting  from 
those  combinations  that  do  support  multipath  have  been  illustrated 
earlier  in  this  report  and  have  been  sampled  in  the  manner  illus- 
trated in  Fig.  3.4.  This  results  in  consideration  of  27  distinct 
channel  conditions.  This  approach  via  the  propagation  software 
provides  us  with  the  discrete  set  of  amplitudes  and  delays  that 
have  been  used  in  the  simulations  described  below. 


Given  a set  of  path  delays  and  amplitudes  obtained  in  the 
manner  described  above,  phases  must  be  determined  for  each  path 
which  will  result  in  a deep  fade.  Referring  to  Section  2.3.2,  it 
has  been  shown  that  for  a set  of  N amplitudes  (aj^)  and  delays 
(ij^)  the  first  three  terms  in  the  frequency  power  series  model 
for  the  channel  transfer  function  are  given  by 

N 

Tp  = L 0.  (3.67) 

i=l  ^ 


N 

T,  = L |3.tj  (3.68) 

^ i=l  ^ 


T„  = ^ L 8.T??  (3.69) 

^ ^ i=l  ^ 1 

where  is  the  complex  amplitude 
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(3.70) 


-j2trf4.  je. 

“ “i®  ""  “i® 

and  7)1^  is  the  delay  after  subtraction  of  the  mean  delay 

T?i  = - €o  (3.71) 

The  value  of  Tq  gives  the  fade  depth  at  the  center  of  the  frequency 
band.  For  a deep  fade  the  value  of  Tq  will  be  quite  small  com- 
pared to  the  magnitude  of  the  vectors  /S^.  This  situation  is 
illustrated  in  Fig.  3.5  for  the  case  N = 3 which  is  of  interest 
here.  From  this  figure  it  is  clear  that  when  Tq  is  small,  small 
variations  in  either  the  magnitude  or  angle  of  the  3i's  can 
produce  significant  variation  in  both  the  magnitude  and  angle  of 
Tq.  On  the  other  hand,  these  small  variations  in  the  3j[*s  will 
not  greatly  affect  T^  and  T2  as  given  by  Eqs.  (3.68)  and  (3.69). 
Based  on  these  considerations,  the  following  procedure  has  been 
used  in  calculating  Tq,  T]^,  and  T2  for  a set  of  gains  and  delays: 

1)  Find  the  values  of  0£  which  give 


2)  Use  these  values  of  9^  to  calculate  T^^  and  T2. 

3)  Select  the  magnitude  of  Tq  to  give  the  desired 
fade  depth. 

4)  Allow  the  angle  of  Tq  to  vary  over  the  range 
(0,  2ff) . For  each  value  in  this  range  the 
coefficients  p^,  P2,  and  pj  may  be  calculated 
and  the  software  used  to  evaluate  error  prob- 
ability. The  error  probabilities  obtained  in 
this  way  are  then  averaged. 


3.2.2  Baseband  Modem  Performance 


For  each  of  the  27  multipath  conditions  performance  of  the 
baseband  modem  was  calculated  at  three  fade  depths:  10,  20,  and 
30  dB.  As  the  depth  of  the  fade  increases,  frequency-selective 
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effects  become  more  pronounced.  In  other  words,  fading  will  be 
"flat"  (that  is  undistorted)  for  smaller  fade  depths  and 
frequency-selective  (that  is  distorted)  for  larger  fade  depths. 

At  some  fade  depth  distortion  will  reach  the  point  where  the 
eye-pattern  closes  for  certain  symbol  sequences;  when  this  happens 
it  will  be  impossible  to  obtain  a satisfactory  error  rate  regard- 
less of  the  value  of  E^/No  available  to  the  receiver.  This 
behavior  may  be  seen  qualitatively  in  Table  3-2  which  summarizes 
modem  performance  for  each  of  the  channel  conditions  considered. 

(In  this  table  the  three  layer  heights  which  were  considered  for 
each  combination  of  refractive  index  gradient  and  layer  tilt 
slope  are  lumped  together  for  simplicity.)  Selecting  a gradient 
and  slope,  one  can  look  across  horizontally  and  compare  perform- 
ance at  fade  depths  of  10,  20,  and  30  dB.  For  relatively  flat 
channels  the  transition  from  flat  fading  to  frequency-selective 
fading  occurs  at  very  deep  fades,  while  for  more  selective 
channels  the  transition  occurs  for  shallower  fade  depths.  Since 
fade  depth  is  quantized  to  10  dB  increments  this  table,  the 
transition  from  flat  fading  to  complete  degr,?dation  is  usually 
abrupt.  Individual  performance  curves  on  which  Table  3-2  is 
based  are  shown  in  Figs.  3.6  through  3.15.  (Figures  in  which 
fading  is  entirely  flat  or  completely  degraded  have  been  elimin- 
ated, so  that  only  ten  curves  rather  than  27  are  required.)  In 
each  case  the  horizontal  axis  is  the  received  E^/Nq  and  a flat 
fading  reference  curve  is  shown  with  a dashed  line.  Thus  the 
degradation  due  to  frequency-selective  fading  is  just  the  hori- 
zontal separation  between  the  curve  in  question  and  the  reference 
curve. 

3.2.3  IF  Modem  Performance 

The  performance  of  the  phase-continuous  FSK  modem  was  eval- 
uated for  the  same  set  of  conditions  for  which  the  baseband  modem 
was  evaluated.  In  making  these  calculations  a bit  rate  of 
27.648  Mb/s  was  assumed  for  equivalence  with  the  baseband  modem, 
although  the  breadboard  of  the  IF  modem  actually  operates  at  a 
bit  rate  of  27.275  Mb/s.  As  before,  flat  fading  performance  is 
shown  with  a dashed  line.  In  this  way  the  loss  in  performance 
due  to  frequency-selective  fading  may  be  compared  between  modems 
even  though  the  absolute  performance  of  the  two  modems  (in  flat 
fading)  is  considerably  different.  Table  3-3  presents  a summary 
of  performance  of  the  IF  modem  which  may  be  compared  with  the 
earlier  table  for  the  baseband  modem.  On  the  whole,  the  effects 
of  multipath  on  the  two  modems  are  remarkably  similar.  However, 
a careful  examination  of  the  tables  or  of  the  accompanying  figures 
(3.16  - 3.26)  shows  that  the  FSK  modem  is  somewhat  less  suscept- 
ible to  multipath  than  the  baseband  modem.  As  an  example.  Fig.  3.17 
shows  less  degradation  than  Fig.  3.8;  both  apply  to  the  same  channel. 
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Figure  3.6  Simulated  Performance  of  DAU  Modem 
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Figure  3.8  Simulated  Performance  of  DAU  Modem 
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Figure  3.11  Simulated  Performance  of  DAU  Modem 
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Figure  3.12  Simulated  Performance  of  DAU  Modem 
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Figure  3.19  Simulated  Performance  of  IF  Modem 
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Figure  3.23  Simulated  Performance  of  IF  Modem 
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Figure  3.24  Simulated  Performance  of  IF  Modem 
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Figure  3.26  Simulated  Performance  of  IF  Modem 
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SECTION  4 


EFFECT  OF  FREQUENCY -SELECTIVE  FADING  ON  DIVERSITY  COMBINING 


This  section  is  devoted  to  an  examination  of  frequency- 
selective  fading  in  a diversity  context.  In  Section  4.1,  the 
propagation  aspects  are  discussed.  The  emphasis  is  primarily 
on  space  diversity,  but  some  discussion  of  frequency  diversity 
(including  a rough  rule  of  thumb  for  choosing  the  frequency 
separation)  is  provided.  A series  of  plots  provides  a graphic 
view  of  diversity  operation  in  a refractive  multipath  environ- 
ment and  some  interesting  conclusions  about  the  advisability  of 
certain  diversity  combining  options. 

Certain  aspects  of  the  tracking  performance  of  modems 
during  diversity  operation  are  discussed  in  Section  4.2. 

4 . 1 Propagation  Effects  in  Diversity  Operation 

In  this  section  we  present  the  results  of  exercising  the 
LOS  software  developed  at  CNR  to  determine  the  effects  on  diver- 
sity operation  of  the  frequency-selective  behavior  of  refractive 
multipath.  In  view  of  the  complicated  dependence  of  refractive 
multipath  on  several  parameters  (compounded  in  a diversity 
context),  this  is  a task  of  considerable  scope.  Nevertheless, 
we  have  been  able  to  generate,  through  a series  of  plots,  an 
information  overview  of  diversity  operation  in  a refractive 
multipath  environment  and  reach  some  important  conclusions 
about  the  advisability  of  choosing  certain  options  for  diversity 
switching  rules. 

A first  step  on  the  road  to  obtaining  the  desired  overview 
is  to  develop  a more  general  view  of  refractive  multipath  on  a 
single  diversity  channel.  This  was  achieved  early  in  this 
program  via  the  fade  level  contour  diagram  illustrated  in 
Figure  4.1.  This  figure  was  pieced  together  from  12  computer 
plots  and  constitutes  a dual  presentation  of  frequency  and 
spatial  selectivity  (operating  frequency  on  the  horizontal  axis 
and  layer  height  on  the  vertical  axis).  The  layer  height  on 
the  vertical  axis  is,  here,  unlike  our  earlier  plots,  measured 
relative  to  highest  terminal  (receiver).  The  reason  for  this 
is  that  our  contour  plots  were  computed  before  the  model  modi- 
fication of  Figure  2.5(b)  which  takes  into  account  layer 


positions  below  the  highest  terminal.  Thus,  there  is  a gap  of 
80'  in  considered  layer  heights  corresponding  to  the  vertical 
region  between  transmitter  and  receiver. 

For  the  situation  described  by  the  specified  parameters 
below  the  figure  caption,  each  contour  in  Figure  4.1  depicts 
the  combination  of  layer  heights  and  frequencies  that  will 
result  in  a fade  to  the  depth  with  which  each  contour  is  desig- 
nated. In  our  initial  investigations  under  this  program,  the 
software  leading  to  Figure  4.1  was  generated  in  order  to  allow 
searches  for  deep  fade  loci  over  extensive  areas  in  layer 
height  and  frequency  space.  As  we  can  see,  the  deep  fade  loci 
are  few  and  far  between;  in  fact,  there  are  only  five  of  them 
appearing  in  Figure  4.1.  Additionally,  we  note  that  the  10-dB 
fade  contours  extend  over  a few  hundred  MHz , whereas  the  20-dB 
fade  contours  are  much  smaller.  This  is  in  keeping  with  experi- 
mental data  pointed  out  by  Kaylor  [4.1]  who  has  noted  that  the 
sharp  deep  fades  occur  during  shallower  fades  that  can  extend 
over  a band  approximately  400  MHz  wide. 

Figure  4.1  was  generated  using  operating  frequencies  in 
the  4-GHz  band.  We  have  performed  similar  computations  for 
the  8-GHz  band  because  of  the  importance  of  the  latter  on  mili- 
tary links.  The  results  are  illustrated  in  Figure  4.2  for 
diversity  receiver  #1  (80  feet  above  transmitter) , and  in 
Figure  4.3  for  diversity  receiver  #2  (55  feet  above  the 
transmitter) . 

An  overlay  of  the  two  figures  results  in  Figure  4.4  where 
we  show  the  20-dB  fade  only.  We  see  that  there  is  no  overlap 
of  the  20-dB  fade  contours.  Although  some  are  close,  as  with 
the  two  pairs  near  20  feet,  there  will  not  be  simultaneous 
center-of -the-band  fading  below  20  dB  on  the  two  diversity 
channels.  It  is  interesting  to  note,  however,  that  the  nodes 
at  22  feet  are  separated  approximately  by  15  MHz;  for  a carrier 
frequency  half  way  inbetween  8.415  MHz),  one  space  diversity 
channel  may  experience  a 20-dB  fade  at  the  lower  end  of  the 
band,  and  the  other  a ZO-dB  fade  at  the  upper  end  of  the  band. 
For  the  current  configuration,  the  effect  on  a 14-MHz  bandwidth 
signal  could  be  marginal.  At  greater  bandwidths , it  would 
become  important. 

The  overall  impression  of  examining  plots  like  those  pre- 
sented above  is  that  the  likelihood  of  simultaneous  fades  on 
both  diversity  channels  is  very  small  (as  it  should  be). 
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GROUND  RANGE  = 30  MILES 

REFRACTIVE  INDEX  GRADIENT  = -420  Nu/km 

LAYER  THICKNESS  > 62.6  FEET 

HEIGHT  OF  RECEIVER  #2  = 55  FEET  (ABOVE  TRANSMITTER) 
Figure  4.3  Deep  Fade  Loci  in  8-GHz  Band  (Receiver  #2) 
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GROUND  RANGE  = 30  MILES 

REFRACTIVE  INDEX  GRADIENT  = -420  Nu/km 

LAYER  THICKNESS  >62.6  FEET 

Figure  4.4  20-dB  Fade  Contours  on  Both  Diversity  Channels 

(Contours  are  numbered  according  to  diversity  channel) 
Receiver  #1  is  80'  above  transmitter;  Receiver  #2  is 
55'  above  transmitter 


In  the  earlier  part  of  this  program,  we  computed  plots  of 
fade  depth  vs.  layer  height,  fade  depth  vs.  frequency,  group 
delay  vs.  layer  height,  etc.,  for  the  physical  situation  used 
to  generate  Figures  4.1  through  4.4.  These  were  presented  in 
the  monthly  reports.  Because  an  extremely  large  number  of  com- 
putational results  can  be  generated  for  any  particular  link,  we 
have,  in  our  attempt  to  limit  the  number  of  curves  presented 
here,  limited  our  calculations  to  the  Swingate-Houten  link 
spanning  the  English  Channel  and  the  Hohenstadt-Zugspitze  link 
in  Germany.  The  curves  presented  below,  however,  are  further 
reductions  in  that  they  represent  "cuts"  of  the  contour  plots 
presented  above  (constant  frequency,  variable  layer  height  or 
constant  layer  height,  variable  frequency  curves).  Such  curves 
were  presented  earlier  in  Section  2,  but  are  examined  here  in 
the  context  of  diversity  operation  on  the  links. 

One  of  the  reasons  for  the  reduction  described  above  is 
that  generation  of  the  contour  plots  presented  above  is  a rela- 
tively time-consuming  and  expensive  procedure  for  locating  the 
deep  fade  loci.  The  other  is  that  we  have  been  able  to  devise 
analytical  techniques  for  reducing  the  amount  of  effort  in 
finding  the  deep  fade  loci.  The  main  analytical  result  in 
this  context  is  that  given  earlier  in  (2.156).  This  result 
states  simply  that,  given  the  amplitudes  and  delays  for  any 
probable  multipath  situation,  one  can  determine  the  frequencies 
at  which  the  deep  fades  occur.  Thus,  one  is  able  to  determine 
the  lines  in  the  two-dimensional  layer  height-frequency  space 
along  which  the  deep  fade  loci  lie.  These  lines,  as  determined 
from  (2.156),  are  presented  as  dotted  lines  in  Figure  4.5. 

They  are  seen  to  pass  exactly  through  the  deep  fade  loci. 


It  is  of  interest  to  examine  a simplification  of  (2.156) 
that  is  suggested  by  the  approximate  linearity  of  all  the  multi- 
path  delay  profiles  we  have  computed  under  this  program;  when- 
ever three  paths  are  present,  the  delay  difference  £0"  = t ^ - t ^ 
between  the  multipath  rays  follows  an  approximate  linear 
dependence  on  layer  height.  To  make  use  of  this,  we  reproduce 
(2.156)  here  for  convenience; 


+ cos 


f = 


“l  - (“2-^  “3) 


Za^OL^ 


Ztt  £sr 


ts 


n = 0, ±2, . . . (4.1) 


This  expression  gives,  for  any  layer  height,  the  frequencies 
at  which  the  deep  fade  loci  appear.  We  have  already  seen  in 
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Figure  4.5  Frequency  Constraints  on  Deep  Fade  Loci  (The  deep 
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the  contour  plots,  e.g..  Figure  4.5,  that  the  deep  fade  loci 
tend  to  cluster  in  groups  of  two.  These  clusters  correspond 
to  the  choice  of  n in  (4.1)  above.  We  now  show,  using  the 
linearity  of  the  multipath  delay  profiles,  that  the  frequency 
separation  of  the  clusters  at  any  layer  height  can  be  repre- 
sented in  terms  of  two  parameters  £0"  (0) , the  delay  separation 
at  some  convenient  reference  layer  height,  and  the  maximum 

layer  height  supporting  multipath.  After  examination  of  the 
multipath  delay  profiles  (e.g..  Figure  2.41),  it  is  a simple 
matter  to  determine  that 


(‘■max  A) 


(4.2) 


Thus,  the  clusters  occur  every  Af  Hz  where 


Af  = 


max 

h -h^ 
max  L 


1 

Ar(0) 


(4.3) 


For  the  data  used  to  generated  Figure  4.2  , ArwO.4  ns  and 
hmax“59.  So,  at  hj^  = 50'  the  deep  fade  loci  are  16.4  GHz 
apart.  The  smallest  value  is  at  hL  = 0.  There  the  loci  are 
2.5  GHz  apart. 

For  this  example,  the  implication  is  clear;  a plot  of  fade 
level  vs.  frequency  (constant  layer  height)  will  show  deep 
fades  only  once  in  extremely  wide  bands. 

It  is  of  more  interest  to  examine  the  frequency  width  of 
the  fade  contours  in  the  context  of  our  above  discussion. 
Consider  a physical  situation  wherein  the  maximum  layer  height 
supporting  multipath  is  given  by  h^  and  the  reference  height 
delay  difference  is  Ar2^(0);  compare  this  with  another  physical 
condition  of  the  link  that  is  characterized  by  corresponding 
values  h2  and  Ar2(0).  If  we  make  the  reasonable  (but  rough) 
assumption  that  the  fade  contour  widths  in  frequency  are  in  the 
same  proportion  as  the  frequency  separation  between  the  deep 
fade  loci,  we  have 


Recall  that  the  zero  reference  is  at  the  receiver  height. 
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(4.4) 


w (h  ) h /h  - h \ At  (0) 

where  w^  and  W2  are  the  contour  widths.  Thus,  of  course,  is  a 
measure  of  the  relative  frequency  selectivity  in  each  of  the 
fading  conditions.  Note  that  if  we  constrain  ourselves  to  the 
region  = hj.^  = h^^  > 0 where,  as  pointed  out  above,  the  fre- 
quency selectivity  is  greatest,  we  have  simply 

W2(0)  /5(r^(0) 

= ^2(0) 

Thus,  we  see  that  there  is  an  important  dependence  between  the 
frequency  selectivity  of  fades  and  the  delay  width  and  layer 
height  extent  of  the  multipath  delay  profiles.  Examination  of 
the  great  variety  of  multipath  delay  profiles  presented  in 
Section  2 indicates  that  there  can  be  a wide  variety  of 
frequency-selective  fading  on  LOS  links. 

Before  proceeding  with  presentation  of  our  calculated 
curves  for  two  special  links  of  interest,  we  point  out  two 
observations  about  the  contour  plots  above  that,  though  simple, 
deserve  comment. 

First,  the  frequency  separations  between  the  deep  fade 
loci  in  the  plots  does  not  depend  on  the  band  examined,  i.e., 
the  deep  fade  loci  occur  with  the  same  frequency  spacing  in  the 
4-GHz  band  as  they  do  in  the  8-GHz  band. 

Secondly,  the  relative  occurrence  vs.  layer  height  of 
the  deep  fades  does  depend  on  the  frequency  band.  Comparing 
Figure  4.1  and  Figure  4.2,  we  see  that  the  vertical  separation 
of  the  deep  fade  loci  in  Figure  4.1  (4-GHz  band)  is  approxi- 
mately twice  what  it  is  in  Figure  4.2  (8-GHz  band).  This 
dependence  can  be  demonstrated  in  a more  general  sense  by 
appealing  to  the  approximate  linearity  of  the  multipath  delay 
profiles  that  were  previously  used  to  generate  (4.2).  Substi- 
tuting (4.2)  in  (2.155),  we  obtain  an  expression  analogous  to 
(2,156)  that  yields  the  layer  height  separation  of  the  deep 
fade  loci.  Their  location  within  the  layer  height  band  allow- 
ing multipath  is  governed  by 
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wherein  we  see  that  the  layer  height  separation  of  the  deep 
fade  loci  is  inversely  proportional  to  f.  Thus,  spatial 
selectivity  becomes  more  important  at  the  higher  frequencies. 

We  now  present  the  results  of  computations  that  we  have 
done  for  two  links  of  interest:  the  Swingate-Houtem  link  over 
the  English  Channel,  and  the  Hohenstadt-Zugspitze  link  in 
Germany.  We  start  with  the  Swingate-Houtem  link.  This  link, 
as  depicted  in  Figure  4.6,  consists  of  several  possible  diver- 
sity channels.  Because  of  time  limitations,  we  have  only  done 
calculations  for  the  two  emphasized  in  Figure  2.35.  These  con- 
sist of  one  transmitter  and  two  diversity  receivers  separated 
by  30  feet.  The  first  is  located  47  feet  above  the  transmitter, 
the  second,  17  feet  above  the  transmitter.  In  Figure  4.6  we 
illustrate  diversity  operation  on  this  link  when  the  lapse  rate 
is  300  Nu/km  and  the  carrier  frequency  is  4.4465  GHz.  This 
carrier,  one  of  the  operating  carriers  on  the  link,  was  not 
chosen  with  any  particular  interest  in  getting  deep  fades.  It 
is  clear  from  the  figure  that  the  choice  of  diversity  separa- 
tion has  resulted  in  anti-positioning  of  the  deep  fades,  i.e., 
one  can  very  roughly  say  that  the  deeper  fades  on  one  channel 
are  positioned  against  high  signal  levels  on  the  other.  The 
multipath  delay  profile  corresponding  to  this  physical  situation 
was  illustrated  earlier  in  Figure  2.37.  Comparing  the  two,  it 
can  be  seen  that  the  cutoff  points  correspond  and  that  signal 
enhancement  occurs  at  the  points  where  the  ray  paths  merge. 

Similar  dependencies  are  evident  when  a different  carrier 
frequency  is  chosen,  as  shown  in  Figure  4.7  for  a carrier  fre- 
quency of  4.808  GHz.  This  plot  shows  a fade  to  less  than 
-20  dB  at  65.9  feet  (as  a computer  printout  of  the  numerical 
values  shows). 

In  order  to  get  some  insight  into  frequency  diversity 
operation,  we  have  plotted  the  frequency  dependence  of  fade 
depth  and  group  delay  for  the  layer  height  of  65.9  feet  that 
gives  the  deep  fade  mentioned  above.  The  plots,  provided  in 
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Figure  4.6  Diversity  Fade  Levels  vs.  Layer  Height  on 
Swingate-Houtem  Link  (horizontal  layer) 
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Figure  4.7  Diversity  Fade  Levels  Vs.  Layer  Height  on 
Swingate-Houtem  Link 
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Figures  4.8  and  4.9  cover  the  whole  band  of  operating  frequencies. 
Consideration  of  (4.3)  and  Figure  2.37,  (which  yields  = 180 

feet  and  At(0)  1,7  ns)  indicates  that  the  nulls  will  occur  at 

a spacing  of  927  MHz.  Thus,  for  frequency  diversity  operation  a 
spacing  of  460  MHz  would  be  suitable.  This  contention  is  sup- 
ported by  the  figures  though  it  is  clear  that  a slightly  smaller 
spacing  would  not  result  in  serious  degradation. 

In  Figures  4.10  and  4,11  we  show  the  frequency  dependence 
across  the  band  when  a lower  lapse  rate  (200  Nu/km)  pertains. 
Because  of  the  contrast  between  the  two  multipath  delay  profiles 
(Figures  2.36  and  2.37)  we  expect  less  sharp  frequency  dependence 
because  of  the  lower  lapse  rate,  but  an  exact  comparison  using 
(4.4)  cannot  be  made  here  because  of  the  different  layer  heights. 


Curves  similar  to  those  given  earlier  in  Figures  4.6  and 
4.7  for  the  Swingate-Houtem  link  have  been  calculated  for  the 
Hohenstadt-Zugspitze  link.  One  such  curve,  depicting  fade  level 
vs.  layer  height  is  depicted  in  Figure  4.12  for  a lapse  rate  of 
400  Nu/km  and  a layer  slope  equal  to  .009.  (The  corresponding 
delay  profile  was  shown  earlier  in  Figure  2,18  of  Section  2.4.1.) 
Once  again  we  note  the  comforting  fact  that  the  nulls  on  one 
diversity  channel  tend  to  be  located  opposite  points  of  high 
signal  level  on  the  other.  The  effect  is  not  as  dramatic  as 
in  the  earlier  plots,  e.g..  Figure  4.7.  To  some  degree  the  fade 
records  on  each  of  the  branches  tend  more  to  follow  one  another 
than  in  the  previous  plots.  The  suspected  reason  for  this  is 
that  the  diversity  antenna  separation  is  only  about  17  feet. 

It  is  expected  that  more  suitable  choices  of  diversity  antenna 
separation  can  be  found. 

It  is  of  interest  to  compare  the  spatial  selectivity  of 
the  curves  in  Figure  4,12  (Hohenstadt-Zugspitze  link)  and 
those  in  Figure  4.6  (Swingate-Houtem  link).  Recall,  from  our 
discussion  in  connection  with  (4.6)  that  cycles  of  height 
variation  have  length 


A = 


h 

max 

fAT(O) 


(4.7) 


In  this  expression  the  frequencies  have  already  been  specified. 
The  quantities  hjnax  At(0)  can  be  determined  from  the  respec- 

tive multipath  profiles;  these  are  Figure  2. 18  for  the  Hohenstadt- 
Zugspitze  link,  and  Figure  2.37  for  the  Swingate-Houtem  link. 

When  this  is  done  we  find 
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Figure  4.8  Fade  Level  Vs.  Frequency  Over  Total  Band  for 
Horizontal  Layer  on  Swingate-Houtem  Link 
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Figure  4.9  Group  Delay  Vs.  Frequency  Over  Total  Band  for 
Horizontal  Layer  on  Swingate-Houtem  Link 
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Figure  4.11  Frequency-Dependence  Over  the  Total  Band  on  Swingate-Houtem  Link 


r 


o 


^ ^ 1 

lS:iS.C  1533.3  15'-.1.  7 1550..;  1555.  5 1555.  ~ 

^ E L H T I '/  E L P I E K H E 1 1.-’  H ■ C F T i 


(ABOVE  TRANSMITTER) 

LAPSE  RATE  = 400  Nu/km 

LAYER  SLOPE  = .009 

CARRIER  FREQUENCY  = 8.591  GHz 


Figure  4.12  Diversity  Fade  Levels  Vs.  Layer  Height  on 
Hohenstadt-Zugspitze  Link 
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= 4.1 


(4.8) 


i.e.,  the  spatial  variation  on  the  Hohenstadt-Zugspitze  link  is 
approximately  four  times  as  rapid  as  it  is  on  the  Swingate- 
Houten  link.  This  is  due  primarily  to  the  increased  value  of 
the  factor  f^rCO)  on  the  Hohenstadt-Zugspitze  link. 

Though  time  has  not  permitted  a precise  validation  of  the 
numerical  value  in  (4.8)  comparison  of  Figure  4.12  and  Figure 
4.6  over  a 50-foot  interval  (e.g.,  the  whole  abscissa  in  Figure 
4.12  and  120  - 170  feet  in  Figure  4.6)  indicate  clearly  that 
fading  on  the  Hohenstadt-Zugspitze  link  is  much  more  rapid  spa- 
tially. If  one  is  willing  to  assume  a rough  equivalence  of  the 
thermo -d3niamic  conditions  leading  to  layer  motion  on  each  link, 
the  natural  inference  is  that  fading  in  time  will  be  more  rapid 
on  the  Hohenstadt-Zugspitze  link.  No  equivalence  need  be 
assumed  when  comparing  different  physical  conditions  on  the  same 
link;  it  exists  a-priori. 

For  completeness,  we  present  plots  (Figure  4.13)  of  group 
delay  on  each  of  the  two  diversity  channels.  Note  that  an 
extremum  of  group  delay  occurs  at  the  same  layer  height 
(•*>  1526  feet)  as  the  deepest  fade  illustrated  in  Figure  4.12. 

We  now  turn  to  a more  microscopic  view  of  the  deep  fade 
loci.  Our  concern  here  relates  to  the  following  question;  what 
are  the  dominant  effects  on  diversity  operation  when  a layer 
moves  through  a deep  fade  locus  and  a diversity  switching  rule 
is  used?  The  primary  interest  here  is  in  examining  some  physical 
mechanisms  that  may  result  in  loss  of  bit  count  integrity. 

We  consider  three  simple  diversity  rules; 

1.  Switch  to  other  channel  if  current  channel 
fades  to  depth  greater  than  30  dB  and  other 
channel  is  stronger. 

2.  Switch  to  other  channel  if  current  channel 
fades  to  depth  greater  than  20  dB  and  other 
channel  is  stronger. 

3.  Switch  to  other  channel  if  stronger  than 
current  channel. 


4-20 


(ABOVE  TRA.NSMITTER) 

Lapse  Rate  " 400  Nu/km 

Layer  Slope  ” .009 

Carrier  Frequency  “ 8.591  GHz 


Figure  4.13 


Diversity  Group  Delays  Vs.  Layer  Height  on 
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Recall  from  our  discussion  of  Section  2.3.3  that  the  largest 
excursions  of  group  delay  tend  to  occur  during  the  deepest  fades. 
This  implies  that  diversity  switching  from  a diversity  channel 
with  a deep  fade  on  it  to  another  channel  with  little  or  no  fade 
is  equivalent  to  switching  from  a channel  with  a large  group 
delay  to  a channel  with  very  small  group  delay.  Because  of  the 
rapid  switching  on  digital  links,  there  may  be  a large  sudden 
discontinuity  in  delay  and  hence  a loss  of  bit  count  integrity. 

For  concreteness  we  look  at  the  fade  in  Figure  4.14  and 
assume  upward  layer  motion*  and  initial  operation  of  diversity 
channel  #1.  In  Figure  4.15  we  have  presented  the  group  delay 
plot  corresponding  to  Figure  4.14  and  marked  off  the  points  at 
which  switching  takes  place  according  to  the  diversity  rules 
listed  above.  For  rule  1 we  get  a discontinuity  of  177  ns. 

For  rule  2 a 50  ns  discontinuity  occurs,  whereas  for  rule  3 it 
can  be  shown  on  an  expanded  scale  that  the  discontinuities  are 
considerably  less  than  10  ns. 

Relative  to  a symbol  length  (72.3  ns)  the  group  Delay  dis- 
continutiy  that  occurs  when  rule  3 is  used  is  extreme.  Any 
conclusion  about  the  disastrous  effects  of  such  a circunstance 
should  be  tempered  with  recognition  of  the  fact  that  the  group 
delay  at  the  center  of  the  band  may  not  best  characterize  the 
overall  delay  of  the  signal.  This  is  because  group  delay  at  any 
one  particular  frequency  becomes  more  tenuous  as  a characterizer 
of  signal  delay  as  frequency-selectivity  increases.  For  this 
reason  statements  that  the  phenomenon  described  above  will  have 
disastrous  effects  on  bit  count  integrity  would  be  premature. 
Nevertheless,  this  phenomena  and  the  fact  that  serious  perform- 
ance degradation  can  take  place  during  20  dB  fades  in  situations 
characterized  by  a high  degree  of  frequency  selectivity  (see 
Section  3),  strongly  suggests  the  use  of  rule  3 above. 


Note  that  regularity  of  layer  motion  need  not  be  assumed  over 
a large  range.  In  fact,  to  think  of  the  abscissa  in  Figure 
4.14  as  a time  axis,  we  need  only  visualize  layer  motion  as 
being  regular  over  a very  small  range,  e.g.,  .3  feet. 
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4.2  Tracking  Performance  of  Modems 


The  modem  performance  calculations  presented  in  Section  3 
required  exact  evaluation  of  modem  tracking  behavior  under  a 
variety  of  multipath  conditions.  In  general,  these  results  show 
that  modem  timing  closely  follows  the  group  delay  of  the  channel 
at  the  carrier  frequency  except  under  conditions  of  extreme 
degradation.  This  is  particularly  true  of  the  baseband  modem. 
For  this  modem  no  significant  deviation  between  center-of- 
band  group  delay  and  modem  tracking  was  found  for  any  of  the 
channels  which  did  not  cause  complete  degradation.  For  channels 
with  an  irreducible  error  rate,  large  tracking  deviations  (as 
much  as  a few  tenths  of  a bit)  were  calculated.  This  is  not 
surprising  since  irreducible  error  rates  arise  from  eye  closure, 
and  eye  closure  for  some  data  sequences  will  shift  the  position 
of  the  mean  zero  crossing  which  is  used  by  the  tracking  loop. 

IF  modem  tracking  also  follows  the  center -of -band  group 
delay,  although  small  deviations  occur.  Figure  4.16  shows  both 
group  delay  and  modem  tracking  as  the  angle  of  Tq  is  varied  over 
360°.  (A  discussion  of  the  significance  of  this  procedure  may 
be  found  in  Section  3.2.1.)  These  results  are  for  the  channel 
conditions  of  curve  #1  of  Fig.  3.16,  although  similar  results 
were  seen  for  all  of  the  channels  analyzed. 
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The  group  delay  has  the  classical  definition 


(2.97) 


where  the  prime  denotes  differentiation  with  respect  to  f.  It 
will  be  helpful  for  us  to  consider  the  above  expression  in  some 
alternate  forms.  Some  elementary  manipulation  of  (2.94)  gives 


(2.98) 


which,  by  writing  T(f,t)  in  terms  of  its  real  and  imaginary 
parts , 


T(f,t)  = X(f,t)  + j Y(f,t) 


(2.99) 


yields 


: .v  ^ X»(f.t)Y(f.t)  - X(f.t)Y*(f.t) 

, , 9 


T (f,W 2 

S 27rL  (f,t) 


(2.100) 


The  above  expressions  constitute  the  definitions  and  for- 
malism that  shall  be  used  in  much  of  the  remainder  of  this 
report.  Some  specialization  of  these  results  is  in  order, 
however.  This  is  because  of  the  fact  that  refractive  multipath 
links  are  known  to  support  a finite  number  of  discrete  paths, 
as  discussed  in  Section  2.1.  This  allows  us  to  use  the  special 
form 


T(f,t)  = 


(t)  e 


-j2lTfT^(t) 


(2.101) 


where  there  are  N paths,  ajL(t)  being  the  time-varying  amplitude 
of  the  i^^  path  and  Tj^(t)  being  its  delay.  As  in  much  of  our 
following  work,  we  shall  keep  in  mind  that  we  are  dealing  with 
a channel  snapshot  and  no  longer  carry  along  an  explicit  depen- 
dence on  t in  our  notation.  Thus,  we  have 


2-54 


T(f)  = 


IN 


-jZJTfT. 


(2.102) 


The  expression  in  (2.102)  is  worth  considering  in  some 
slightly  altered  forms.  We  shall  have  occasion  to  refer  back 
to  these  in  later  sections  of  this  report. 

One  form  derives  from  considering  the  frequency  relative 
to  the  carrier  frequency,  i.e.,  by  considering  dependence  on  a 
relative  frequency 


f = f - f 


(2.103) 


This  has  the  effect  of  replacing  the  real  amplitude  in  (2.102) 
with  a complex  amplitude.  Substituting  (2.103)  in  (2.102) 
we  get 


■■  ■ I; 


-j27rf'T. 


(2.104) 


where 


-i2iTf  T . 

/*!  = e ^ - 


(2.105) 


It  will  also  be  of  interest  to  consider  the  effects  of 
adding  a fixed  delay  offset  to  each  of  the  ray  path  delays. 
Replacing  r with  t+Tq,  we  get 


T(f')  = e 


-j27rf  r„  N 


-jZirfjA  -j2iTf'(TQ-Hr^) 


(2.106) 


In  our  use  of  these  expressions  later  on,  we  will  not  carry 
along  the  prime  on  f;  the  meaning  will  be  clear  from  the 
context . 

Expressions  for  fade  level  and  group  delay  in  terms  of 
the  ray  path  delays  and  amplitudes  can  be  easily  derived 
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1 


using  (2.102).  From  (2.102)  and  (2.95)  we  get 


“i  “j 


cos  2vf(T^ 


(2.107) 


for  the  fade  level.  To  determine  an  expression  for  group  delay, 
we  use  (2.102)  and  (2.98)  and  obtain 


T 

g 


N 


Im{ 


Z 2 01.T  e 
i=l  ^ ^ 


-j2irfT^  \ 


N 


i=l 


-j2fffT. 


e 


(2.108) 


which,  after  some  manipulating,  gives 


T 

g 


a.  a.  T . 
i J 1 


cos  2l7f(T^ 


(2.109) 


This  completes  our  summary  of  the  analytical  results  for 
fade  level  and  group  delay. 

Some  insight  into  the  behavior  of  the  parameters  discussed 
in  this  section  can  be  gleaned  by  considering  the  illus'  n 
of  phasor  geometry  in  Figure  2.12.  The  figure  represt 
graphical  depiction  of  the  way  in  which  the  contributing  ^rms 
in  (2.104)  link  together  to  form  the  resultant  phasor  (boi.d-face 
arrow).  The  magnitude  of  this  phasor  is  the  fade  level,  and  it 
is  seen  in  the  figure  to  vary  substantially  across  the  2 A band, 
decreasing  at  f ^ - A to  half  the  value  it  had  at  f ^ + A.  In  con- 
junction with  this  variation,  the  resultant  traces  a phasor 
trajectory  in  the  xy  plane  and,  in  so  doing,  rotates.  It  is 
the  variation  of  the  central  angle  0 with  respect  to  frequency 
that  results  in  the  occurrence  of  group  delay,  as  is  clear  from 
(2.97). 

It  is  interesting  to  note  that  insights  into  the  statis- 
tical behavior  of  fading  can  be  garnered  from  Figure  2.12. 
Consider  a situation  wherein  the  x and  y components  of  the 
resultant  phasor  in  Figure  2.12  are  random  variables  with  a 
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probability  distribution  given  by  p(x,y) . About  this  situation 
we  can  make  the  perfectly  reasonable  assumption  that  no  pre- 
ference is  given  to  the  infinite  fade  point  at  the  origin,  i.e., 
that  there  are  no  singularities  of  p(x,y)  at  the  origin.  This 
implies  smooth  behavior  of  p(x,y)  in  the  neighborhood  of  the 
origin.  We  are  interested  in  the  probability  of  fading  below 
some  level  Lg;  the  result,  of  course,  is  just  the  probability 
that  the  x and  y components  are  such  as  to  put  the  resultant 
phasor  in  a circle  centered  at  the  origin.  If  the  fade  is  deep 

enough,  i.e.,  if  Lq  is  small  enough,  this  probability  is  given 
by 


P(LsLq)  « Lq  p(0,0)  (2.110) 


This  intuitive  but  correct  observation  corresponds  to  a result 
derived  by  Lin  [2.20]. 

2.3,2  Signal  Distortion  in  Multipath  Channels 

The  propagation  channel  is  a time-variant  linear  system. 
In  complex  notation,  the  input-output  relationships  are  fre- 
quently presented  with  the  aid  of  the  time-variant  impulse 
response  g(t,4)  or  the  time-variant  transfer  function  T(f.t). 
as  follows: 


w(t)  = J Z(f)T(f,t)eJ^’^^*^  df 
w(t)  - J z(t-0g(t,0  d^ 


(2.111) 


where  w(t)  is  the  output;  z(t),  Z(f)  are  the  input  process 
and  its  spectrum;  and  T(f,t),  g(t,4)  are  related  by  Fourier 
transformation 


T(f.t)  - J d4 

g(t,€)  - T(f.t)e-^^’^^*'  df 


(2.112) 
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In  calculating  modem  performance  we  are  generally  Interested 
In  situations  where  the  degree  of  frequency  selectivity  displayed 
by  the  channel  over  the  bandwidth  of  the  data  signal  is  moderate. 
This  is  particularly  true  for  the  two  high  packing  density  modems 
to  be  considered  later,  since  high  packing  density  is  achieved  at 
the  expense  of  sensitivity  to  noise  and  intersymbol  interference. 
Consequently,  a few,  or  frequently  only  the  first  three  terms  in 
a power  scries  expansion  of  the  transfer  function  suffice  to 
characterize  the  frequency  selectivity.  The  utility  of  such  a 
power  series  to  model  a time-variant  dispersive  channel  by  a chain 
of  differentiations  has  been  studied  by  Bello  [2.19].  Formally, 


T(f,t)  = S T (t)(2irjf)"  (2.113) 

n“0 

which  represents  the  channel  by  a sum  of  the  outputs  of  success- 
ively higher  orders  of  (bandpass)  differentiators,  with  the  nth- 
order  differentiator  being  multiplied  by  the  complex  function 


T„(t) 


1 

6'^(f.t) 

n!(2iTj)" 

J (-e)"g(t.«)  dc 


(2.114) 


Bello  [2.21]  has  shown  that  the  power  series  representation 
(2.113)  will  be  most  rapidly  convergent  when  the  expansion  is 
taken  for  a channel  in  which  a ’’mean*'  path  delay  has  been  removed. 
Assuming  this  mean  path  delay  is  £q,  the  input-output  represen- 
tation corresponding  to  (2.113)  (including  the  mean  path  delay)  is 
given  by 
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00 


(2.115) 


w(t)  = E T (t) 
n=0  " 


d z(t  - 4n) 


dt 


n 


This  parametric  model  is  shown  in  Fig.  2.13.  (This  chain- 
differentiator  model  is  called  the  f-power  series  model  in  [2.21] . 

When  the  frequency-selective  fading  in  the  channel  is  suf- 
ficiently small,  only  the  first  term  in  the  series  (2 . 115)  will  be 
sufficient  to  characterize  the  channel  output,  i.e.. 


w(t)  = Tjj(t)z(t  - Cq)  (2.116) 

where  z(t),  w(t),  are  the  complex  envelopes  of  the  channel  input 
and  output,  respectively.  Equation  (2.116)  maybe  recognized  as  a 
"flat  fading"  or  nonfrequency-selective  channel  model.  If  the 
first  two  terms  are  used. 


w(t)  = TQ(t)z(t  - ^q)  + T^(t)z(t  - €q)  (2.117) 


which  is  called  a "linearly  frequency-selective  fading"  channel 
[2.22]  since  it  corresponds  to  approximating  T(f,t)  by  a complex 
linear  term  in  the  frequency  variable.  One  may  continue  and 
define  a "quadratically  frequency-selective  fading"  channel: 

w(t)  = TQ(t)z(t  - 4q)  + T^(t)z(t  - 4q)  + T2V(t  - 4q) 

(2.118) 

This  is  the  representation  which  will  be  used  for  performance 
modeling . 

Given  any  channel  impulse  response  structure  one  may  compute 
the  parameters  (Tn(t)}  of  the  chain  differentiator  model.  We  are 
particularly  concerned  with  an  N-path  channel  for  which  the  time- 
varying  channel  transfer  function  is  given  by 
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Figure  2.13 


T(f,t) 


(2.119) 


N 


-j2n.(t)[f + f ] 


where  f^,  is  the  carrier  frequency  and  f is  the  frequency  devia* 
tion  from  the  carrier.  Since  the  existence  of  a mean  delay 
produces  a factor  exp  (-j2irf4Q)  in  T(f,t)  which  can  fluctuate^ 
with  f quite  rapidly,  it  is  desirable  to  expand  only  the  portion 
of  T(f,t)  which  does  not  include  this  factor.  Then 


-j2irfCo 

T(f,t)  = TQ(f,t)e  ^ 


(2.120) 


where  TQ(f,t)  is  the  time-varying  channel  transfer  function  after 
the  mean  delay  has  been  removed.  From  Eqs.  (2.119)  and  (2.120) 
(f,t)  for  the  N-path  model  is 


N -j2fff^4  (t)  - j2irf[4  (t) 

T„(f,t)  = S a. (t)e  ^ 10  (2.121) 

i=l 


Defining  a complex  gain  of  the  i'th  path  by 

-j2irf  ti(t) 

|3^(t)  = a.(t)e  ^ 


(2.122) 


and  a relative  delay  of  the  i'th  path  by 

\(t)  - {^(t)  - 


(2.123) 


then 


N -j2fffT;  (t) 

T (f,t)  = E i8.(t)e 
i=l  ^ 


(2.124) 
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It  can  be  easily  shown  from  Eq.  (2.121)  that  the  parameters  of 
the  model  are  given  by 


Tn(t)  = 


iliiZ 

n! 


N 

L 

i=l 


i3i(t)[T,^(t)]" 


(2.125) 


We  have  shown  above  how  distortion  of  the  complex  envelope 
z(t)  is  related  to  the  coefficients  of  the  chain  differentiator 
model  {Tn(t)]  and  how  these  coefficients  may  be  calculated  for 
an  N-path  channel.  We  now  find  expressions  for  the  corresponding 
distortion  of  the  modulating  waveform  when  z(t)  is  either  phase  of 
frequency  modulated.  In  either  case  we  may  write 


z(t) 


ej<P(t) 


(2.126) 


In  the  case  of  phase  modulation,  the  modulating  waveform  is  given 
directly  by  <p(t);  in  the  case  of  frequency  modulation  it  is  given 
by 


x(t) 


(2.127) 


Using  (2.127)  and  the  quadratic  channel  model  [Eq.  (2 . 118) ] we  find 
that  the  complex  envelope  of  the  channel  output  is  given  by 


jo(t  - 4q) 

w(t)  = e ]TQ(t)  + jT^(t)cp(t-4Q) 

+ T2(t)[j<p(t  - 4q)  - 0^(t  - ^q)] 


(2.128) 


For  the  LOS  channel  of  interest,  time  variations  of  the  channel 
are  much  slower  than  time  variations  of  o(t)  and  its  derivatives. 
Thus,  one  may  regard  the  channel  as  fixed  when  finding  the  dis- 
criminator output.  Writing  w(t)  in  the  form 


2-63 


w(t)  = A(t)TQ(t)e 


j[<p(t  - 4n)  + B(t)l 


(2.129) 


the  desired  results  for  output  angle  and  frequency  are 
0(t)  = ^(t  - + B(t) 


(2.130) 


y(t)  - i «t)  - xtt-ig)  B(t) 


(2.131) 


By  a simple  calculation  it  may  be  shown  that  for  small  distortion 
B(t)  is  given  by 

B(t)  = p^<p(t  - ^q)  + p^oVt  - 4q)  + i P30^(t  - 4q)  (2.132) 

where  the  p parameters  are  defined  by 

T^(t) 

Pi  " 


T2(t) 


^^2  TQ(t) 


T (t) 

T2(t) 

To(t) 

To(t) 

To(t) 

(2.133) 


In  this  result  derivatives  (and  their  products)  of  order  higher 
than  2 have  been  discarded,  maintaining  accuracy  consistent  with 
the  quadratic  channel  model.  Substituting  (2.133)  into  (2.130) 
gives  the  phase  modulation  model  directly: 


8(t)  = <p(t  - ^q)  + p^«p(t  - €q)  + - 4p)  + ^ p o^(t  - 4 ) 

(2.134) 
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and  substituting  into  (2.131)  with  the  help  of  (2.132)  gives 
the  frequency  modulation  model 


y(t)  = x(t  - 4q)  + p^x(t  - 4q)  + P2V(t  - 4q)  + P3x(t  - 4Q)x(t  - ^q) 

(2.135) 


These  models  will  be  the  basis  for  the  modem  performance 
analyses  presented  in  Section  3. 

2.3.3  Fade  Depth  and  Group  Delay  Extrema 

In  this  section  we  explore  and  answer  some  important 
questions  relating  to  the  frequency  dependence  of  fade  depth 
and  group  delay.  Two  questions  are  of  utm  ,st  importance: 

(1)  Given  a fixed  set  of  channel  parameters,  at 
what  frequencies  do  the  deep  fades  occur? 

(2)  What  is  the  joint  behavior  of  fade  level  and 
group  delay?  Do  large  excursions  of  group 
delay  usually  occur  during  deep  fades? 

The  first  question  is  of  importance  because  of  the  fact 
that  the  nonfading  situations  are  of  little  interest  to  us. 
Because  of  the  improbability  of  deep  fades,  it  is  helpful  to 
have  analytical  results  that  can  aid  one  in  locating  the  deep 
fade  loci.  These  are  presented  in  this  section. 

The  importance  of  the  second  question  relates  to  loss  of 
bit  count  integrity  during  diversity  switching.  Since  many 
diversity  schemes  involve  switching  during  deep  fades,  a large 
group  delay  excursion,  occurring  simultaneously  with  the  deep 
fade,  could  result  in  a large  rapid  discontinuity  occurring 
when  the  system  is  switched  to  the  other  channel.  In  this 
section  we  develop  a bound  on  the  group  delay  which  indicates 
that  the  largest  group  delay  excursion  will  only  occur  during 
deep  fades . The  importance  of  this  phenomenon  will  be  dis- 
cussed later  in  Section  4.1. 

Some  interesting  insights  into  the  problem  at  hand  are 
gained  by  considering  the  geometrical  picture  in  Figure  2.14. 
Here  the  transfer  function  at  f - f]_  is  depicted  as  the  sum  of 
its  component  phasors,  exp(-j2irf  ^) , etc. 
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clearly,  as  the  frequency  changes,  the  resultant  phasor  T(f) 
traces  out  a trajectory  in  the  XY  plane.  This  trajectory  is 
illustrated  by  the  dotted  line  in  Figure  2.14.  It  is  of  use 
to  define  a "velocity"  vector  that  describes  the  "motion"  of 
the  resultant  phasor  along  the  trajectory.  We  use  these  words 
in  quotation  marks  (to  be  discarded  in  future  discussion) 
because  we  are  not  dealing  with  true  velocity  and  motion  in 
the  sense  of  kinematics.  Dependence  here  is  on  frequency 
rather  than  time;  nevertheless,  the  concept  of  velocity  de- 
veloped here  is  a useful  one.  In  terms  of  the  transfer  function 
which  we  here  separate  into  its  component  parts  and  write  as 
a vector. 


T(f)  = X X + Y y 


(2.136) 


where  x and  y are  the  unit  vectors  in  the  X and  Y directions 
as  illustrated  in  the  figure,  the  velocity  vector  is  defined 
simply  as 


V = V X + V y 
X y 


V = X' (f)  V = Y' (f) 

X y 


(2.137) 


where  the  primes  denote  differentiation  with  respect  to  f . 

It  is  also  of  interest  to  refer  to  the  absolute  speed 
along  the  trajectory.  This  is  simply 


V 


+ v2) 

y' 


(2.138) 


The  velocity  vector  at  f = f 2 illustrated  in  Figure  2.14. 
The  angle  0 appearing  in  the  figure  is  the  same  as  that  referred 
to  in  Eq . (2 . 94)  . 

A simple  examination  of  the  geometry  in  Figure  2.14  and 
comparison  with  (2.97)  reveals  that  the  absolute  value  of  the 
group  delay  at  any  frequency  is  inversely  proportional  to  2ffL 
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X 


Figure  2.14  Phaser  Geometry 
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and  directly  proportional  to  the  absolute  value  of  that  com- 
ponent of  the  velocity  vector  tangent  to  a circle  centered  at 
the  origin  and  passing  through  the  point  (X,Y).  To  firm  this 
up  analytically,  we  define  a unit  vector  perpendicular  to  T(f), 
[i.e.,  tangent  to  the  circle  of  radius  L at  the  point  (X,Y)]; 


u 


t 


- Y X + X v 

Vx^  + Y^ 


It  is  easy  to  see  that  u^  has  the  desired  properties 


(2.139) 


“t  • “t  ‘ 1 


u^  • T(f)  ==  0 


(2.140) 


where  we  have  used  to  dot  product  of  vector  analysis.  It  is 
an  easy  matter  to  see  from  (2.100)  and  the  geometry  in 
Figure  2.14  that 


T 

g 


(2.141) 


which  is  verified  by  substituting  (2.137)  and  (2.138)  into  this 
expression  to  get  our  earlier  result 


^ X Y'  - X'  Y 
^ 2n(X^  + Y^) 


(2.142) 


We  will  return  to  this  formalism  later,  but  one  interesting 
observation  can  be  made  immediately.  Consider  a given  point 
(XqjYq)  and  a resultant  phasor  moving  along  a trajectory  passing 
through  the  point  (Xq,Yq)  in  such  a way  that  it  crosses  the 

fade  circle  of  radius  L=  (Xq  + Y^j^  . For  the  value  of  f at 
which  the  trajectory  crosses  the  fade  circle,  the  magnitude  of 
the  group  delay  is  always  less  than  the  value  it  would  have  if 
the  resultant  phasor  were  moving  (at  the  same  speed)  along  a 
trajectory  tangent  to  the  fade  circle.  This  is  a simple  con- 
sequence of  the  inequality 
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(2. 143) 


Iv 


a:  V 


and  we  shall  make  use  of  it  later. 


We  now  address  question  (1)  mentioned  at  the  beginning 
of  this  section,  and  determine  those  frequencies  at  which  the 
deepest  fades  occur.  From  (2.102)  it  is  easy  to  see  that  two 
necessary  conditions  for  the  occurrence  of  fades  as  deep  as, 
or  deeper  than,  that  corresponding  to  some  magnitude  squared 
value  l2  are; 


and 


2irfT . L 
1 


(2.144) 


t 

= 1 


a. 

1 


sin 


2vfT  . ^ L 
1 


(2.145) 


Thus,  assigning  a small  value  to  L places  a constraint  on  the 
values  of  f (for  a given  set  of  and  T^)  allowing  deep  fades. 

It  is  of  interest  to  examine  the  constraints  on  f when 
L = 0.  This  is  done  by  solving  the  system  of  equations 


Sa.  cos  2tj£t  . = 0 


Sa.  sin  2vfT . = 0 


(2.146) 

(2.147) 


This,  of  course,  is  the  constraint  condition  resulting  from 
directly  setting  T(f)  in  (2.102)  equal  to  0.  Generally,  it  is 
difficult  to  solve.  For  the  case  N = 3,  some  specific  conclu- 
sions can  be  made. 


Since  this  problem  is  by  no  means  one  of  the  major  goals 
of  this  program,  we  concentrate  then  on  a three  path  model. 

To  find  the  infinite  fade  points,  we  set 
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-jZnfT . 


a. 


0 


(2.148) 


Manipulation  of  this  expression  is  facilitated  if  we  factor  out 
the  first  term.  Since  a is  never  0,  the  constraint  (2.148) 
becomes  ^ 

-j2fffT’  -j2irfT' 

1+a^e  +a^e  =0  (2.149) 


where , for  j = 2 , 3 


= a. /a, 
3 1 


and 


(2.150) 


t! 

J 


(2.151) 


Separating  the  real  and  imaginary  parts  of  (2.149),  we  obtain 


1 + 2jrfT^  + cos  Zrrfr^  = 0 


(2.152) 


and 


02  sin  2iTfr2  + sin  277fT2  = 0 


(2.153) 


Squaring  and  adding  both  of  these  equations,  we  obtain 

1 - 


cos  2fff (t2  - T^) 


• \ = 


2 


(2.154) 
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or 


cos  2it£(t^  - T^)  = 


2 I 2^  2\ 

“l  - (“2 + “3) 


^“2  “3 


From  (2.155)  we  have 


f = 


^ -1 
+ cos 

[2  1 2 . 2\\ 
“1  • (“2+ “3^ 

^ ®2 

2ff(T3  - T^) 

(Tj-V 

where 


(2.155) 


(2.156) 


n = 0,  +1,  +2,  ...  (2.157) 


Thus,  it  is  clear  for  the  model  under  consideration  that  two 
deep  fades  may  occur  every  (.t^-T2)~^  Hz  in  frequency,  e.g., 
for  a delay  difference  on  the  order  of  one  nanosecond  there  can 
at  most  be  two  nulls  across  a 1-GHz  band.  We  point  out,  how- 
ever, that  a simple  examination  of  (2.155)  and/or  the  phasor 
geometry  indicates  that  the  constraint  on  f can  never  be 
satisfied  unless  the  conditions 


“2  ^ “3  ^ “1  (2.158) 

and 

Ia2-a3i  ^ (2.159) 

are  met . 

It  should  be  pointed  out  that  the  use  of  (2.156)  to  locate 
a deep  fade  does  not  guarantee  that  a deep  fade  will  actually 
occur  at  that  frequency;  it  only  means  that  deep  fades,  should 
they  occur,  will  occur  at  the  given  frequencies. 


I 


I 
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We  now  turn  from  a determination  of  the  deep  fade  loci 
to  an  examination  of  the  joint  behavior  of  the  fade  level  and 
group  delay.  This  is  in  response  to  question  (2)  at  the 
beginning  of  this  section;  our  analytical  goal  is  to  determine 
bounds  on  the  group  delay. 


From  (2  '8)  we  see  that 


g 


-j-T  |lm[T'(f)T*(f)|2 


(2.160) 


Using  (2.99),  we  obtain 


|T'(f)T*(f)l^  = L^[(x')^+(y')^] 


(2.161) 


With  the  definitions  for  v and  v given  in  (2.137)  and  the 
additional  definition  ^ 


/ 2 

V = yv^  + V 


(2.162) 


which  can  be  thought  of  as  the  "speed"  of  the  resultant  along 
the  phase  trajectory  in  Figure  2.14,  we  obtain  the  simple 
inequality 


1 |2  ^ V 

^ 2 2 
® 4ff  L 


Evaluation  of  v^  and  v^  according  to  (2.137)  gives 


(2.163) 


N 


V = - 2 

X 


V = - 2 

y 


n 7 a.  T . si 


” I;  “1  ^ 


. sin  2rrfT_ 


cos  2nfT 


(2.164) 


(2.165) 


2-72 


Application  of  Lagrange's  identity 


(2.166) 


gives 


% = <^’'>11;  “j 


2 2 

■7  sin  2iTfT. 
J J 


N N 


- ^ " °lc^j  (2.167) 


with  a similar  expression  for  v . Adding  the  two  expressions 
gives  ^ 


(2.168) 


where 


N Nr 

S Sh 


Q = (2ir)‘ 


r,  sin  lufT,  - a,  t . sin  lirfr 
k k k J : 


+ 2 2iTfT 


(2.169) 


By  doing  the  indicated  squaring  in  (2.169),  comining  terms,  and 
substituting  back  into  (2.168),  it  is  a simple  matter  to  see  that 


N N 


= (2ir)2  ^ a.a.  T.r.  cos[27rf(T^-Tj)] 


(2.170) 
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Thus,  from  (2.163)  the  group  delay  is  bounded  according  to  the 
following  inequality: 


T 

g 


N 

E 

i=l 


“i'"i' 


(2.171) 


We  have  derived  inequalities  by  other  means  than  that  used  here 
(e.g.,  by  applying  the  Schwarz  inequality  to  each  of  the  terms 
in  v^  = v^  + v^)  but  in  all  cases  these  inequalities  have  been 
found  to  be  weaker  than  (2.171). 

Examination  of  the  steps  leading  to  (2.171)  indicate  that 
the  bound,  if  it  is  to  be  achieved,  can  only  be  achieved  at 
the  fade  depth  minimum. 

As  a numerical  example  we  have  taken  some  computed  values 
of  amplitudes  and  delays  (typical  for  the  shorter  links).  These 
are; 


= 1 “2  ~ 1*18^  “3  = 1.196 

= 0 = 1.266  = 1.740 


(2.172) 


Across  a 30-MHz  band  centered  on  a carrier  frequency  of 
4.9888  GHz,  the  selected  set  of  amplitudes  and  delays  above 
gives  rise  to  a minimum  fade  level  given  by 


D = -43.457  dB 


Using 


we  get 

T 

g 


= 281.7  ns 


(2.173) 


(2.174) 


(2.175) 


max 
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as  the  maximum  group  delay  excursion  corresponding  to  the 
set  of  parameters  in  (2.172).  This  compares  favorably  with 
a computed  value  of  240  ns. 

The  example  leading  to  (2.175)  is  a worst-case  situation, 
of  course;  we  have  calculated  the  bound  at  the  frequency  which 
gives  rise  to  the  deepest  fade.  The  bound  in  (2.171)  is  useful 
in  other  respects;  e.g.,  it  indicates  the  relative  unimportance 
of  group  delay  excursions  during  shallower  fades.  As  a case  in 
point,  the  group  delay  upper  bounds  for  10  dB , 20  dB , and  30  dB 
fades  are  5.983  ns,  18.920  ns,  and  59.83  ns,  respectively. 
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2.4  Two  Specific  Communication  Channels  of  Interest 


Using  the  analytical  results  described  in  Section  2.2  ve 
have  computed  delays  and  amplitudes  on  two  specific  links  of 
interest,  the  Hohenstadt-Zugspitze  link  in  Germany  and  the 
Swingate-Houten  link  which  crosses  the  English  Channel  at  a 
point  near  Dover.  Additionally  we  have  calculated  the  dependence 
of  fade  depth  and  group  delay  on  frequency  and  layer  height  for 
these  two  links. 

Because  modem  performance  on  the  Hohenstadt-Zugspitze  link 
will  be  evaluated  at  a later  point  in  this  report  we  have  empha- 
sized calculation  of  the  multipath  delays  and  amplitudes  on  that 
link.  It  will  be  recalled  that  the  set  of  multipath  amplitudes 
and  delays  provides  a complete  characterization  of  the  link  from 
which  all  other  quantities  of  interest  can  be  calculated,  e.g., 
fade  depth,  group  delay,  and  distortion  parameters.  Most  sig- 
nificantly, the  delays  and  amplitudes  can  be  used  to  calculate 
error  probability  and,  hence,  determine  the  SNR  degradation  of 
the  frequency-selective  channel  relative  to  a flat-fading  chan- 
nel. The  results  of  calculations  of  the  SNR  degradation  are 
provided  in  Section  3. 

2.4.1  Hohenstadt-Zugspitze  Link 

The  Hohenstadt-to-Zugspitze  (H  Z)  and  Zugspitze-to- 
Hohenstadt  (Z  — H)  link  are  of  interest  because  of  their  rather 
extreme  geometries  relative  to  typical  earth-based  LOS  microwave 
links.  As  illustrated  in  Fig.  2.15  the  two  terminals  are  sepa- 
rated by  approximately  100  miles  with  the  terminal  at  Zugspitze 
elevated  over  a mile  relative  to  that  at  Hohenstadt.  (It  should 
be  noted  that  all  heights  are  referred  to  the  height  of  the 
primary  antenna  at  Hohenstadt.*)  The  intuitive  motion  that  such 
an  extreme  geometry  can  lead  to  extreme  fading  is  supported  by 
the  program  of  measurements  described  in  [2.8].  It  is  reported 
there  that  fading  over  the  H — Z path  is  greater  than  over  any 
other  tested  in  the  DEB  Stage  I system.  Therefore,  it  was  not 
surprising  for  us  to  encounter  more  extreme  values  of  multipath 
than  we  have  encountered  on  shorter  links  that  we  have 
investigated . 


The  "primary"  antenna  is  the  one  antenna  at  either  location 
that  is  enabled  to  both  transmit  and  receive. 


2-76 


7021.20 


Figure  2.15  Link  Geometry  for  Hohenstadt-Zugspitze  Link 


Before  describing  the  results  of  the  calculations  vie  should 
point  out  that  there  are  two  features  of  the  H -*  Z link  which, 
in  one  important  sense,  tend  to  balance  one  another  out.  These 
are  the  extreme  length  of  the  link  and  the  extreme  height  of  the 
receiver  relative  to  the  transmitter.  Simple  physical  consider- 
ations indicate  that  an  increased  ground  range,  all  other  things 
being  constant,  tends  to  increase  the  delay  spreads.  On  the 
other  hand,  elevating  the  receiver  and  holding  all  other  para- 
meters constant  can  lead  one  eventually  to  a situation  wherein 
there  is  no  multipath.  This  is  particularly  noticeable  when  one 
assumes  that  the  refractive  layer  generating  the  multipath  is 
perfectly  horizontal.  As  one  increases  the  receiver  height,  the 
elevation  angle  of  the  transmitted  electromagnetic  field  becomes 
so  great  that  the  refractive  layer,  in  order  to  support  multipath, 
must  have  a capability  to  bend  the  ray  that  exceeds  the  capabil- 
ities of  atmospheric  layers  occurring  in  nature.  In  fact,  we 
have  determined  that  a perfectly  horizontal  atmospheric  layer 
will  not  support  multipath  on  the  H -•  Z link  unless  its  lapse 
rate*  exceeds  800  Nu/km.  For  virtually  all  of  the  useful  quan- 
titative data  that  has  come  into  our  hands  (e.g. , [2.5],  [2.6], 
and  [2.7])  this  value  is  considerably  off  scale. 

Were  it  not  for  the  fact  that  layers  tilt  [2.17],  [2.18]  the 
fading  on  the  H — Z path  would  be  unexplainable  via  the  mechanism 
of  refractive  multipath  in  atmospheric  layers.  As  it  happens, 
the  layer  tilt  angle  is  a variable  of  major  importance  for  links 
with  the  extreme  geometry  of  the  H Z link. 

The  tilted  layer  configuration  was  illustrated  earlier  in 
Fig.  2.10  and  inclusion  of  layer  tilt  into  the  computer  model  was 
discussed  in  Section  2.2.3.  As  pointed  out  in  that  section  data 
on  layer  tilt  is  meagre.  Kerr  [2.18]  has  pointed  out  that  warm 
fronts  have  tilt  slopes  ranging  between  .0033  and  .010,  corres- 
ponding to  an  angular  range  .190  degree  to  .573  degree. 

We  have  used  this  range  of  slopes  in  our  calculations  choosing  a 
.001  slope  increment  to  cover  the  range  from  .003  to  .010.  The 
corresponding  angles  are  provided  in  Table  2-3. 

Results  are  shown  in  Figs.  2.16,  2.17,  and  2.18  for  three 
different  values  of  the  lapse  rate  of  refractive  index  gradient, 
300,  350,  and  400  Nu/km.  The  curves  in  these  figures  are  the 


By  lapse  rate  we  are  referring  to  the  lapse  rate  of  refractive 
index  gradient.  Algebraically  this  quantity  is  just  the  nega- 
tive of  the  refractive  index  gradient. 
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TABLE  2-3 


VALUES  OF  LAYER  SLOPE  USED  IN  CALCULATIONS 
AND  CORRESPONDING  LAYER  TILT  ANGLES 


Layer 

Slopes 

Layer 

Tilt  Angle 
(Degrees) 

.003 

.172 

.004 

.229 

.005 

.286 

.006 

.344 

.007 

.401 

.008 

.458 

.009 

.5156 

.010 

.573 
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Zugspitze  Link  for  a Lapse  Rate  of  300  Nu/km 


Lapse  Rate  = 350  Nu/km 


(su)  XB-[9a 
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Figure  2.17  Multipath  Delay  Profile  of  Diversity  Channel  No.  1 on  the  Hohenstadt-to- 
Zugspitze  Link  for  a Lapse  Rate  of  350  Nu/km 


Lapse  Rate  = 400  Nu/km 


o 


(su)  i^Piaa 
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Figure  2.18  Multipath  Delay  Profile  of  Diversity  Channel  No.  1 on  the  Hohenstadt-to- 
Zugspitze  Link  for  a Lapse  Rate  of  400  Nu/km 


multipath  delay  profiles;  they  represent  the  delay  of  each  ray 
path  supported  by  the  propagation  channel.  The  delays  are 
plotted  vs.  layer  height  and  are  measured  relative  to  the  delay 
of  the  direct  path  channel.  As  an  example,  note  from  Fig.  2.18 
that  a layer  with  slope  .010  supports  three  ray  paths  when  it 
is  at  a height  of  1500  feet.  For  two  of  these  paths  there  are 
delay  lags;  one  is  at  3 ns,  the  others  at  6 ns . For  the  other 
path  there  is  a delay  advance  of  about  -1  ns.  As  the  layer  moves 
up  this  path  delay  converges  on  zero,  reaching  that  point  when 
the  layer  rises  above  the  receiver  and,  by  so  doing,  allows  a 
direct  path  to  exist  between  transmitter  and  receiver.  For  this 
reason  we  refer  to  it  as  the  pseudo-direct  path.  As  the  layer 
rises  it  is  also  clear  from  the  figure  that  the  upper  two  ray 
paths  converge;  it  should  be  emphasized  that  the  point  of  con- 
vergence is  by  no  means  the  same  as  the  point  of  convergence 
between  the  direct  and  pseudo-direct  path.  Another  point  that 
should  be  emphasized  is  that  the  delay  profiles  presented  in  this 
report,  for  economy  of  presentation,  only  show  the  delays  that 
exist  in  a multipath  situation,  i.e.,  when  there  is  more  than  one 
delay. 

One  interesting  feature  of  each  of  these  curves,  e.g.. 

Fig.  2.17,  is  that  the  multipath  profiles  tend  to  open  up  or 
spread  as  the  slope  of  the  layer  increases.  This  is  because  the 
layer  is  approaching  the  position  of  maximal  multipath  which 
occurs  when  the  layer  and  the  straight  line  connecting  the  trans- 
mitter and  the  receiver  are  perfectly  parallel.  This  position  is 
never  quite  reached  here  since  the  straight  line  slope,  at  a value 
of  .013,  is  slightly  in  excess  of  .010,  the  largest  layer  slope 
of  interest. 

Another  feature  of  Figs.  2.16  to  2.18  is  that  increases  in 
lapse  rate  allow  for  the  occurrence  of  multipath  at  smaller  tilt 
angles.  This  is  not  surprising  in  view  of  our  knowledge  of  the 
general  fact  that  the  large  lapse  rates  tend  to  be  more  amenable 
to  support  of  multipath  than  the  smaller  ones.  An  extreme  case 
is  shown  in  Fig.  2.19  where  we  have  used  a lapse  rate  value  of 
500  Nu/km.  It  can  be  seen  that  multipath  support  exists  over 
layer  slopes  ranging  from  .005  to  .010  with  delay  spreads  some- 
times in  excess  of  20  ns.  Fortunately,  this  lapse  rate  is  not 
within  the  range  expected  for  the  H -•  Z link.  (See  comments  in 
[2.8]  and  data  on  page  107  of  [2,6].) 
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Figure  2.19  Multipath  Delay  Profile  of  Diversity  Channel  No.  1 
on  the  Hohenstadt-to-Zugspitze  Link  for  a Lapse 
Rate  of  500  Nu/km 


2-84 


The  dependence  on  lapse  rate  is  more  clearly  demonstrated 
in  Fig.  2.20  below  where  we  show  how  the  multipath  delay  profile 
changes  as  the  lapse  rate  is  increased.  From  an  extremely  small 
profile  when  the  lapse  rate  is  350  Nu/km  the  profile  is  seen  to 
evolve  to  a condition  wherein  there  are  multipath  delay  separa- 
tions in  excess  of  20  ns  when  the  lapse  rate  is  550  Nu/km.  For 
the  curves  in  Fig.  2.20  the  layer  slope  is  kept  fixed  at  a value 
of  .008. 

For  Figs.  2.16  through  2.20  it  should  be  pointed  out  that 
each  profile  in  each  of  the  figures  is  plotted  only  over  the 
range  of  layer  heights  supporting  multipath.  In  actuality  the 
single  negative  path  extends  indefinitely  to  the  left  whereas  the 
pseudo-direct  path  (see  Fig.  2.18)  extends  beyond  the  profile 
and  steadily  increases  until  it  becomes  zero.  This  is  the  point, 
of  course,  at  which  the  layer  rises  above  the  receiving  terminal 
and  direct  path  propagation  becomes  possible.  This  observation 
coupled  with  examination  of  Figs.  2.16  through  2.20  indicates 
that  the  refractive  multipath  model  would  have  been  totally 
inadequate  for  the  extreme  geometry  situation  encountered  on  the 
H -»  Z link  were  it  not  for  the  model  up-date  described  in  Section 
2. 2. 2. 2 and  illustrated  in  Fig.  2.5(b). 

When  the  delay  calculations  are  complemented  with  amplitude 
calculations  the  LOS  refractive  channel  characterization  is 
complete.  The  plots  of  amplitude  vs.  layer  height  do  not  have 
as  great  an  effect  as  their  delay  counterparts  do,  and  we  have 
restricted  the  presentation  to  the  amplitudes  for  a lapse  rate 
of  400  Nu/km  as  shown  in  Fig.  2.21.  The  layer  slopes  range  from 
.007  to  .010  as  in  Fig.  2.18.  As  the  tilt  increases  the  same 
horizontal  spreading  that  was  observed  in  the  delay  profiles  is 
observed  here.  It  should  be  noted,  however,  that  the  amplitude 
increases  for  the  smaller  tilt  angles.  This  is  because  the  delay 
profiles  are  "tighter"  at  the  smaller  tilt  angles,  as  pointed  out 
earlier.  This  means  as  the  layer  height  passes  through  that 
range  of  values  that  allows  the  existence  of  multipath,  the 
delays  and  amplitudes  of  each  of  the  ray  paths  are  closer  in 
value.  In  the  geometrical  sense  this  means  that  the  rays  are 
closer  to  one  another  in  space.  This  latter  is  a focussing 
condition  and  makes  for  higher  signal  amplitudes  on  each  path. 

Figure  2.21  is  rather  densely  packed.  The  nature  of  the 
amplitude  dependence  on  layer  height  is  more  clearly  displayed 
in  Fig.  2.22.  This  figure  shows  both  the  amplitude  and  delay 
profile  for  a lapse  rate  of  400  Nu/km.  We  can  easily  identify 
the  pseudo-direct  path  as  the  one  that  converges  on  zero  (layer 
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Figure  2.20  Multipath  Delay  Profile  of  Diversity 
Channel  No.  1 on  the  Hohenstadt-to- 
Zugspitze  Link.  Layer  slope  is  .008. 
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Figure  2.21  Multipath  Amplitude  Profile  of  Diversity  Channel 


No.  1 on  the  Hohenstadt-to-Zugspitze  Link  for  a 
Lapse  Rate  of  400  Nu/km 
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Figure  2.22  Multipath  Profiles  on  Hohenstadt-Zugspitze  Link  When  Lapse  Rate  is  400  Nu/km 
and  Layer  Slope  is  .010 


height  increasing)  in  the  delay  profile  and  converges  on  unity 
in  the  amplitude  profile. 

The  data  described  above  has  been  used  as  input  parameters 
for  the  modem  performance  evaluations  discussed  in  Section  3. 

For  easy  reference  the  method  of  choosing  the  parameters  and  the 
parameter  choices  are  provided  with  a separate  section  (Section 
2.4.3). 

With  the  delays  and  amplitudes  described  above  we  have 
calculated  (in  addition  to  the  performance  calculations  presented 
in  Section  3)  two  basic  fading  parameters,  the  fade  level  and 
the  group  delay.  We  now  turn  to  the  presentation  of  several 
curves  that  illustrate  the  basic  nature  of  the  fading. 

One  such  curve  is  presented  in  Fig.  2.23.  This  plot  illus- 
trates the  dependence  of  fade  level  on  layer  height  for  an 
arbitrarily  chosen  carrier  frequency  of  8.2  GHz.  This  figure 
should  be  examined  in  connection  with  the  .009  multipath  delay 
profile  of  Fig.  2.18.  Comparing  the  two  we  see  that  the  fading 
occurs  over  the  same  interval  of  layer  height  that  is  character- 
ized by  a multiplicity  of  rays  in  Fig.  2.18.  This  range  extends 
from  approximately  1520  feet  to  1920  feet.  At  the  higher  end 
of  the  range  there  is  the  focussing  condition  (enhanced  amplitude) 
that  we  have  mentioned  before.  In  the  middle  of  the  range  there 
are  three  fades  that  go  down  to  about  -10  dB.  The  fact  that 
deeper  fades  do  not  occur  is  not  surprising  in  view  of  the  dearth 
of  deep  fade  loci  evident  in  Fig.  1.1  and  our  arbitrary  choice 
of  carrier  frequency. 

From  the  figure  it  can  be  seen  that  the  signal  level  is  not 
zero  dB  outside  the  fading  interval.  The  reason  for  this  is 
that  the  one  remaining  propagation  path  still  intersects  the 
layer,  i.e.,  the  highest  terminal  is  still  immersed  in  the  layer. 
Thus  there  is  still  some  signal  enhancement;  at  the  higher  layer 
locations  propagation  takes  place  via  the  "quasl-direct"  path 
of  Fig.  2.5(b).  The  amplitude  decreases  at  the  higher  levels 
until  finally  the  layer  rises  above  the  highest  receiver, 
direct  path  propagation  is  allowed  to  take  place  and  the  signal 
amplitude  (always  measured  relative  to  the  direct  path)  becomes 
0 dB. 

Since  refractive  layers  are  thought  to  move  vertically 
[2.12],  Fig.  2.23  can,  to  some  extent,  be  viewed  as  a dynamic 
plot.  This  is  true  also  of  subsequent  plots  (depicting  fading 
behavior  as  a function  of  layer  height)  to  be  presented  in  this 
report.  In  fact,  it  is  of  interest  to  temporarily  assume  some 
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Figure  2.23  Fade  Level  Vs.  Layer  Height  on  Hohenstadt-to 
Zugspitze  Link 
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constancy  of  speed  and  direction  for  the  layer  and  view  the 
figures  as  plots  depicting  dependence  on  time.  Viewed  in  this 
light  Fig.  2.23  and  later  plots  bear  more  than  a superficial 
resemblance  to  fading  records  of  various  channel  parameters. 

The  computer  plot  of  Fig.  2.23  is  limited  in  resolution 
because  of  the  fact  that  it  was  generated  from  only  200  points. 
Magnification  of  the  plot  reveals  that  a fade  deeper  than  any 
one  of  those  depicted  in  Fig.  2.23  occurs  near  1900  feet.  The 
magnification  of  this  fade  is  shown  in  Fig.  2.24.  It  can  be 
seen  that  the  fade  extends  slightly  deeper  than  -20  dB. 

We  recall  from  the  discussion  of  Section  2.  that  the 
largest  group  delay  excursions  are  allowed  to  take  place  during 
deep  fades.  For  this  reason  it  is  of  interest  to  examine  the 
dependence  of  group  delay  over  the  same  interval  of  layer  height 
shown  in  Fig.  2.24.  This  illustration  is  provided  in  Fig.  2.25. 
We  see  that  the  extrema  of  group  delay  (about  100  ns)  and  fade 
level  coincide. 

As  we  discuss  in  Section  4 the  group  delay  in  Fig.  2.25  does 
not  apply  to  all  the  frequency  components  of  a signal  traversing 
the  refractive  channel.  Such  a consideration  raises  questions 
that  are  best  addressed  by  Figs.  2.26  and  2.27.  There  we  depict 
the  frequency  dependence  of  the  fade  level  and  group  delay  across 
a 30  MHz  band  centered  on  the  8.2  GHz  carrier.  There  we  see 
that  some  frequency  components  of  the  signal  are  hardly  attenu- 
ated at  all,  whereas  others  are  down  20  dB;  some  frequency 
components  are  hardly  time  shifted  at  all  whereas  others  arrive 
with  a 100  ns  time  shift.  It  is  comforting  to  note,  however, 
that  the  frequency  components  most  delayed  are  also  those  that 
are  most  attenuated.  Thus,  neglect  of  the  distortion  effects 
of  frequency  components  10  dB  down  at  the  receiver  is  equivalent 
to  neglect  of  frequency  components  delayed  by  more  than  about 
10  ns.  Though  such  considerations  are  of  interest  they  are  prob- 
ably best  addressed  through  the  T and  P parameter  formalism 
developed  in  Section  2.3.2  and  applied  extensively  in  Section  3. 
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Figure  2.24  Fade  Level  Vs.  Layer  Height  on  Hohenstadt-to- 
Zugspitze  Link 
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Figure  2.25  Group  Delay  Vs.  Layer  Height  on  Hohenstadt-to- 
Zugspitze  Link 
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Figure  2.26  Fade  Level  Vs.  Frequency  on  Hohenstadt-to- 
Zugspitze  Link 
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Figure  2.27  Group  Delay  Vs.  Frequency  on  Hohenstadt-to- 
Zugspitze  Link 
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When  the  analytical  results  of  Section  3 are  applied  to 
determine  the  frequency  (for  any  given  layer  height)  at  which 
the  deep  fades  occur,  much  more  extreme  conditions  can  be  found. 
Running  of  a computer  program  based  on  (2.156)  indicates  a deep 
fade  possibility  at  8.591  GHz.  As  demonstrated  in  Fig.  2.28  a 
fade  record  for  this  carrier  frequency  yields  a fade  exceeding 
30  dB  at  a layer  height  of  approximately  1526  feet.  A 
"magnification"  of  this  fade  is  provided  in  Fig.  2.29  and  a 
group  delay  plot  follows  in  Fig.  2.30.  The  frequency-selective 
behavior  is  illustrated  in  Figs.  2.31  and  2.32.  These  plots 
show  extreme  values  (approaching  400  ns)  for  group  delay. 

Things  can  be  made  to  look  even  worse.  Running  of  our 
programs  to  locate  the  deep  fade  loci  (in  frequency)  indicates 
that  a choice  of  carrier  of  8.301  GHz  can  result  in  a fade  in 
excess  of  40  dB  at  a layer  height  of  1606.5  feet.  This  is 
illustrated  in  Fig.  2.33.  The  corresponding  group  delay  (yield- 
ing extreme  excursions  up  to  750  ns)  is  illustrated  in  Fig.  2.34. 

The  numerical  values  presented  in  this  section  can  have  the 
appearance  of  sometimes  being  extreme.  This  was  something 
anticipated  at  the  onset,  however.  The  Hohenstadt-Zugspitze 
link  is  an  "extreme"  link.  It  is  longer  and  has  a greater 
terminal  separation  than  any  other  link  to  which  we  have  com- 
fortably applied  our  computer  model.  It  is  encouraging  to  note, 
however,  that  measurements  performed  on  the  link  also  give  rise 
to  "extreme"  empirical  data  [2.8];  they  indicate  a greater 
incidence  of  multipath  fading  than  any  other  of  the  DEB  links. 

Aspects  of  diversity  performance  on  this  link  are  deferred 
until  Section  4. 

2.4.2  Swingate-Houtem  Link 

In  this  section  we  present  the  results  of  using  the  LOS 
software  to  calculate  multipath  delay  and  amplitude  profiles 
and  to  investigate  frequency-selective  and  spatially-selective 
fading  characteristics  of  the  Swingate-Houtem  link.  This  link 
consists  of  a very  long  (55  mile)  line-of-sight  microwave  path 
across  the  English  Channel  from  Swingate  near  Dover,  England  to 
a site  near  Houten,  Belgium. 

Some  refractivity  data  for  this  link  has  been  presented  in 
[2.23].  The  data  was  taken  at  Shoeburyness  which  is  located 
approximately  50  km  northwest  of  Dover  along  the  English  coast 
and  indicates  that  the  choice  of  refractivity  gradient  -300  Nu/km 
that  we  have  used  in  many  of  our  calculations  below  is  entirely 
reasonable.  2-96 
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Figure  2.28  Fade  Level  Vs.  Layer  Height  on  Hohenstadt-to- 
Zugspitze  Link 
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Figure  2.29  Fade  Depth  Vs.  Layer  Height  (Expanded  View) 
on  Hohenstadt-Zugspitze  Link 
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Figure  2.30  Group  Delay  Vs.  Layer  Height  on  Hohenstadt-to 
Zugspitze  Link 
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Figure  2.31  Fade  Depth  Vs.  Frequency  on  Hohenstadt-to-Zugspitze 
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Figure  2.32  Group  Delay  Vs . Frequency  on  Hohenstadt-to- 
Zugspitze  Link 
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Figure  2.33  Fade  Level  Vs.  Layer  Height  on  Hohenstadt-to-Zugspitze  Link 
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Figure  2.34  Group  Delay  Vs.  Layer  Height  on  Hohenstadt-to- 
Zugspitze  Link 
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The  geometry  of  the  Swingate-Houtem  link  is  illustrated  in 
Fig.  2.35.  The  two  terminals  are  separated  by  55  miles.  The 
positions  of  the  two  top-most  receiving  terminals  relative  to 
the  "a"  transmitting  terminal  are  47  feet  (receiver  "A")  and 
17  feet  (receiver  "b").*  The  operating  frequencies  are  in  the 
4 GHz  band  extending  from  4.465  GHz  to  4.908  GHz. 

The  calculations  presented  in  the  remainder  of  this 
section  relate,  because  of  time  limitations,  to  the  two  top- 
most diversity  paths,  i.e.,  we  consider  the  "A"  transmitter  to 
"a"  receiver  path  (diversity  channel  #1)  and  the  "A"  transmitter 
to  "B"  receiver  path  (diversity  channel  #2) . 

Multipath  delay  profiles  (plots  of  multipath  delay  vs. 
layer  height)  for  the  Swingate-Houtem  link  are  illustrated  in 
Figs.  2.36  through  2.41.  Figures  2.36  through  2.38  show  the 
delays  arising  on  diversity  channel  #1  (47  feet  above  the 
transmitter)  and  Figs.  2.39  through  2.41  show  the  delays  arising 
on  diversity  channel  #2  (17  feet  above  the  transmitter).  For 
each  diversity  channel  the  figures  are  sequenced  in  such  a way 
as  to  show  the  effects  of  increasing  lapse  rate. 

In  Fig.  2.36  we  see  that  the  multipath  rays  are  very  close 
in  delay  with  two  extra  paths  being  supported  over  a very  small 
interval  of  layer  height  near  15  feet.  Multipath  is  supported 
for  layer  heights  up  to  no  more  than  about  40  feet.  The  lapse 
rate  used  to  calculate  the  profile  in  Fig.  2.36  is  200  Nu/km. 
This  low  value,  which  is  just  above  that  which  supports  ducting 
(157  Nu/km)  is  the  reason  for  the  tightness  of  the  profile. 

In  Fig.  2.37  we  see  the  effects  of  an  increase  in  lapse 
rate  to  300  Nu/km.  The  profile  spreads  significantly  and 
support  of  multipath  propagation  continues  to  a layer  height  in 
excess  of  180  feet. 

The  effects  of  a further  increase  in  lapse  rate  are  illus- 
trated in  Fig.  2.38.  The  trends  already  discussed  are  evident 
here  to  an  even  greater  degree.  Significantly,  a lapse  rate  of 
400  Nu/km  here  supports  multipath  up  to  340  feet. 


The  zero  reference  level  is  at  sea  level.  The  elevations  of 
the  two  antenna  sites,  not  given  in  Fig.  2.35,  are  provided  in 
[2.231  . 
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Figure  2.36  Multipath  Delay  Profile  on  Swingate-to-Houtem  Link  for  a Horizontal 
Layer 
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Figure  2.39  Multipath  Delay  Profile  on  Swingate-to-Houtem  Link  for  a Horizontal  Layer 
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Figure  2.40  Multipath  Delay  Profile  on  Swingate-to-Houtem  Link  for  a Horizontal  Layer 
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Similar  trends  are  evident  in  Figs.  2.39  - 2.41  which  show 
the  effects  on  diversity  channel  #2  (receiver  at  relative  height 
of  17  feet)  which  is  closer  to  the  transmitter  height.  Two 
basic  differences  are  observable,  however.  These  are: 

1)  The  "extra"  multipath  rays  exist  over  a more 
substantial  interval  of  layer  height  than  pre- 
viously. This  is  because  the  imaginary  straight 
line  connecting  transmitter  "A"  and  receiver  "B" 
is  more  nearly  horizontal  than  in  the  previous 
case.  This  situation,  as  discussed  earlier,  is 
more  amenable  to  the  existence  of  extra  paths. 

2)  The  multipath  support  disappears  at  a lower 
layer  height  than  in  the  previous  case.  This 
is  because  the  layer  separates  from  the  lower 
receiver  "B"  earlier  than  it  does  from  receiver 
"A". 

The  ray  path  amplitudes  for  the  situation  illustrated 
earlier  in  Fig.  2.37  are  presented  in  Fig.  2.42.  Comparing  the 
two  figures,  we  note  that  a focussing  condition  (Fig.  2.42) 
takes  place  as  the  ray  path  delays  converge  (Fig.  2.37). 
Additionally,  we  note  that  the  "extra"  paths  in  Fig.  2.37, 
though  negligible  in  delay  difference  relative  to  the  more 
enduring  multipath  rays  can  be  substantial  in  amplitude.  The 
implication  of  this  is  that  the  extra  short-lined  multipath 
rays  can  contribute  substantially  to  signal  enhancement  but  only 
rarely  play  any  significant  role  during  deep  fades.  In  fact, 
the  disappearance  of  these  extra  paths  results  in  discontinuities 
of  fade  level  as  we  will  point  out  later. 

Because  of  the  fact  that  the  amplitudes  seem  to  play  less 
of  a role  as  indicators  of  deep  fading  conditions  we  have  chosen 
not  to  display  the  amplitude  profiles  as  extensively  as  the 
delay  profiles  presented  earlier. 

All  of  the  above  curves  depict  phenomena  that  depend 
strictly  on  the  physical  character  of  the  link  and  its  geometry. 
None  of  them  depend  on  frequency.  We  now  turn  from  the 
physical  character  of  the  link  to  its  fading  character.  This 
latter,  of  course,  depends  on  frequency. 
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In  order  to  be  fair  and  not  over-emphasize  the  most  serious 
cases  of  signal  degradation  [obtained  by  using  (2.156)  to 
locate  the  frequencies  at  which  the  deepest  fades  occur]  we  have 
assumed  for  most  of  the  calculations  presented  in  the  remainder 
of  this  section  the  carrier  frequencies  specified  for  the 
Swingate-Houten  link.  As  an  example  we  present  Figs.  2.43  and 
2.44  which  depict  a fading  condition  existing  for  the  4.4465  GHz 
carrier  frequency.  We  point  out  that  the  asymmetry  of  the  group 
delay  plot  in  Fig.  2.44  is  a characteristic  that  we  have  observed 
repeatedly  in  our  calculations. 

For  the  operating  carrier  frequency  4.808  GHz  plots  of 
fade  level  vs.  layer  height  are  illustrated  in  Fig.  2.45  for 
diversity  receiver  #1  and  in  Fig.  2.46  for  diversity  receiver 
it^2.  The  increase  in  range  of  layer  heights  over  which  the 
multipath  exists  is  evident  as  the  lapse  rate  of  refractive 
index  gradient  increases.  Also  evident  in  these  figures  is  the 
signal  enhancement  that  takes  place  when  the  extra  paths  exist. 

The  disappearance  of  these  paths  at  the  higher  end  is  accomp- 
anied by  a sharp  rise  in  signal  amplitude.  This,  of  course,  is 
due  to  the  merging  of  the  extra  rays  which  results  in  a focussing 
condition  as  mentioned  before.  A similar  phenomena  involving 
the  dominant  rays  is  evident  in  each  figure  at  the  highest  end 
of  the  layer  height  interval  over  which  multipath  is  supported. 

In  nature  these  sharp  increases  may  be  dampened  somewhat  by 
concurrent  reductions  of  lapse  rate  as  the  layer  moves  upward. 

The  last  figure  of  the  sequence  described  above.  Fig.  2.46(c) 
indicates  a deep  fade  occurring  at  65.9  feet.  The  frequency 
dependence  of  this  fade  is  depicted  in  detail  in  Fig.  2.47. 

We  see  clearly  that  a very  large  fade  has  occurred  on  diversity 
channel  #2  whereas  diversity  channel  #1  is  in  a condition  of 
signal  enhancement.  All  of  the  frequency  components  in  a 14  MHz 
band  centered  on  the  4.808  GHz  carrier  fade  to  -16  dB  or  deeper. 

At  the  center  of  the  band  the  fade  level  is  approximately 
-30  dB.  This  is  the  nominal  fade  level  that  we  have  had 
occasion  to  refer  to.  At  adjacent  frequencies,  however,  the 
frequency  components  can  be  attenuated  to  as  far  down  as  -45  dB. 
Since  these  adjacent  frequencies  are  well  within  the  band  it  is 
a moot  question  whether  the  overall  signal  level  is  best  char- 
acterized by  the  nominal  fade  depth  at  the  center  of  the  band  or 
by  the  behavior  at  some  adjacent  frequency.  In  any  event  the 
fading  is  characterized  by  a higher  degree  of  frequency 
selectivity . 
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Figure  2.44  Fade  Level  Vs.  Layer  Height  for  Swingate-Houtem 
Link  (Horizontal  Layer) 
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Figure  2.45(a)  Fade  Level  Vs.  Layer  Height  for  Swingate-Houtem  Link  (Horizontal  Layer) 
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Figure  2.45  (c)  Fade  Level  Vs.  Layer  Height  for  Swingate-Houtem 

Link  (Horizontal  Layer) 
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Figure  2.46(b)  Fade  Level  Vs.  Layer  Height  for  Swingate-Houtem  Link 
(Horizontal  Layer) 
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Figure  2.46  (c)  Fade  Level  Vs.  Layer  Height  for  Swingate-Houtem 

Link  (Horizontal  Layer) 
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The  forementioned  frequency  selectivity  is  also  apparent 
in  the  plot  of  group  delay  vs.  frequency  in  Fig.  2.48.  We  see 
that  extreme  excursions  are  possible  within  the  band.  The 
extremely  large  group  delay  excursions  must  be  considered  in 
conjunction  with  the  earlier  fade  depth  vs.  frequency  curve. 
Doing  this  we  see  that  the  frequency  components  most  delayed 
are  also  those  most  attenuated.  Thus,  the  effects  on  signal 
distortion  at  the  receiver  are  not  clear.  Such  questions,  of 
course,  were  the  motivation  for  the  P and  T parameter  represen- 
tations discussed  earlier  in  Section  2.3.2. 

To  further  examine  the  extreme  values  that  can  be 
encountered  we  have  chosen  a carrier  frequency  just  in  excess 
of  that  used  to  generate  the  earlier  curves.  The  resulting 
dependence  of  fade  level  and  group  delay  on  layer  height  is 
depicted  in  Figs.  2.49  and  2.50  where  we  note  fade  depths  to 
-63  dB  and  group  delay  excursions  in  excess  of  2000  ns.  As 
pointed  out  in  connection  with  the  earlier  curves,  the  narrow 
frequency  domain  occupancy  of  the  most  extreme  excursions  and 
the  congruence  of  fade  depth  and  group  delay  excursions  may 
tend  to  negate  the  drastic  extremes  appearing  in  plots  of 
center-of-the-band  quantities. 

2.4.3  Propagation  Parameters  Used  for  Computations 
of  Modem  Performance 


In  order  to  compute  the  SNR  degradation  due  to  frequency- 
selective  fading  on  the  Hohenstadt-to-Zugspitze  link  we  have 
taken  sample  points  from  each  of  the  multipath  profiles  illus- 
trated in  Figs.  2.16  to  2.18,  The  layer  height  samples  are 
given  approximately  by 


h,  = h . + n 

L min 


max 


min 


n = 1,2,3 


(2.177) 


where  hj^^n  is  the  smallest  layer  height  resulting  in  multipath 
for  a given  layer  tilt  and  h^ax  is  the  largest.  This  scheme  is 
illustrated  in  Fig.  2.51  where  we  have  numbered  the  cuts  accord- 
ing to  the  value  of  n in  (2.177).  This  numbering,  accompanied 
by  the  exact  layer  heights  involved,  also  appears  on  the  error 
rate  curves  to  be  presented  later  in  this  report. 
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Figure  2.48  Group  Delay  Vs.  Frequency  for  Swingate-Houtera  Link 
(Horizontal  Layer) 
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Figure  2.49  Fade  Depth  Vs.  Layer  Height  on  Swingate-Houtem  Link 
(Horizontal  Layer) 
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Lapse  Rate  = -300 


Carrier  Frequency  = 4.81  GHz 

Figure  2.50  Fade  Depth  Vs.  Layer  Height  on  Swingate-Houtem  Link 
(Horizontal  Layer) 
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Hohenstadt-to-Zugspitze  Link 


Extensive  computer  simulations  of  the  DAU  modem  and  IF  modem 
operating  over  the  channel  described  above  have  been  carried  out 
using  the  simulation  software  developed  under  this  contract.  This 
software  is  based  upon  the  frequency  power  series  model  of  the 
multipath  channel  and  uses  the  distortion  parameters  P]^,  P2, 
and  p^  (see  Section  2.3.2)  to  characterize  the  distortion  in  the 
discriminator  output.  The  SNR  degradation  due  to  frequency- 
selective  fading  is  the  topic  of  Section  3. 

Before  presenting  performance  results  we  briefly  summarize 
the  process  of  identifying  channel  parameters  for  simulation. 
First,  for  the  geometrical  parameters  we  utilize  the  Hohenstadt- 
to-Zugspitze  specifications  described  earlier.  For  the  refrac- 
tive index  gradient  the  values  -400,  -350,  and  -300  Nu/km  are 
used;  these  are  in  reasonable  accord  with  relevant  data  [2.6]. 

For  each  value  of  refractive  index  gradient  the  layer  slopes 
given  in  Table  2-3  have  been  used  in  the  calculations.  This 
results  in  consideration  of  27  distinct  channel  conditions. 

Not  all  combinations  of  refractive  index  gradient  and  layer 
slope  result  in  multipath,  however.  All  of  the  multipath  pro- 
files resulting  from  those  combinations  that  do  support  multi- 
path  have  been  illustrated  earlier  in  this  report  and  have  been 
sampled  in  the  manner  illustrated  in  Fig.  ^.51.  This  approach 
via  the  propagation  software  provides  us  with  the  discrete  set 
of  amplitudes  and  delays  that  have  been  used  in  the  simulations 
described  in  Section  3. 

2 . 5 Comparisons  With  Experiment 

In  this  section  we  briefly  compare  our  computed  results 
with  experimental  data  and  theoretical  conclusions  of  other 
workers  in  the  field.  Most  of  our  comparisons  relate  to  work 
done  by  investigators  of  the  Bell  Telephone  Company. 

Our  comparisons  with  the  available  literature  reveal  no 
significant  differences  between  our  computed  results  and  exper- 
imental measurements — in  fact,  the  overwhelming  impression 
imparted  in  such  a comparison  is  that  our  results  show  the  same 
physical  dependencies  as  the  experimental  data. 

As  an  example  of  the  above  noted  correspondence  between 
experimental  results  and  our  computed  results  we  point  out  that 
measurements  reported  by  Babler  [2.25]  show  no  disagreement, 
in  terms  of  fundamental  numerical  values  with  our  computed 
quantities.  In  fact,  his  measured  plots  of  amplitude  vs.  fre- 
quency, bear  a striking  resemblance  to  the  ones  we  have  calcu- 
lated. Babler  reports  on  statistical  measurements  of  linear 
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and  quadratic  distortions  of  amplitude.  The  quantities  he 
deals  with  derive  strictly  from  consideration  of  signal  ampli- 
tude and  do  not  include  any  information  on  the  important  fre- 
quency variation  of  phase.  Similar  limitations  apply  to  the 
diversity  data  in  [2.26].  As  a continuation  of  the  forementioned 
program,  however,  the  limitations  were  removed  in  the  work  of 
Subramanian,  et  al.  [2.27]  who  present  fade  records  in  con- 
junction with  records  of  phase  difference.  These  records 
(Figures  5 and  6 in  [2.27])  bear  a close  resemblance  to  the 
plots  of  fade  level  and  group  delay  that  we  have  calculated  for 
links  of  comparable  length,  and  yield  similarly  close  numerical 
values .* 

The  authors  of  [2.27]  conclude  their  paper  with  the 
statement : 

"Although  simple  two-ray  models  can  be  made  to 
account  for  both  amplitude  and  phase  dispersion, 
the  complex  temporal  behavior  indicated  in  our 
results  along  with  that  reported  by  Babler  lead 
us  to  believe  that  the  multi-ray  (more  than  two 
rays)  phenomenon  is  the  cause  of  these  deep  fades." 

(A  two-ray  model,  of  course,  would  yield  a simple  sinusoidal 
dependence  on  frequency.)  The  above  "belief ,"  in  our  model,  is  a 
direct  consequence — one  direct  path  and  two  layer  paths  being 
the  most  probable  situation  for  condition  "A"  of  Fig.  2.5,  (see 
demonstration  of  this  at  the  end  of  Section  2.2.4),  and  one 
pseudo-direct  path  and  two  layer  paths  being  the  standard  case 
for  condition  "B",  [see  discussion  in  connection  with  the  cubic 
equation  (2.65)]. 

Relative  to  the  importance  of  frequency-selectivity  on 
microwave  LOS  links  we  note  that  Babler,  as  previously  pointed 
out  by  Kaylor  [2.28],  points  out  the  increased  importance  of 
frequency-selectivity  during  deep  fades.  This  is  something  that 
has  stood  out  clearly  in  the  analysis  (see  Section  2.3.3),  and 
in  virtually  all  of  the  computations  in  this  report. 


By  "phase  difference"  is  meant  the  difference  in  phase  between 
two  tones  .55  MHz  apart.  Though  this  resolution  is  suitable 
for  estimates  of  group  delay  (frequency  derivative  of  phase) 
on  the  link  measured  in  [2.27]  the  group  delay  curves  shown 
there  are  not  as  sharp  as  our  computed  ones. 
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Another  example  of  the  appropriateness  of  our  model  relates 
to  the  occurrence  of  signal  enhancement-  This  seems  not  to  have 
been  a major  concern  (justifiably  so)  of  the  BTL  workers. 
Naturally,  for  a phenomena  that  doesn't  result  in  signal  degrad- 
ation, there  is  no  profound  interest.  Nevertheless,  the  second- 
ary results  of  a useful  propagation  model,  though  of  little 
interest,  should  correspond  with  real  world  results.  As  pointed 
out  in  [2.29]  microwave  paths  that  support  the  formation  of 
ground  based  atmospheric  layers  may  cause  signal  enhancement  up 
to  30  dB.  We  have  observed  signal  enhancements  up  to  24  dB  in 
the  calculation  from  our  model  (see  Fig.  2.46)  and  were  it  not 
for  the  limited  resolution  of  these  plots  (there  was  no  practical 
motivation  for  looking  at  a finer  resolution)  the  enhancements 
would  be  slightly  larger.  Relative  to  occurrence  of  fades  we 
point  out  that  Barnett's  conclusion  [2.30]  that  the  probability 
of  deep  multipath  fading  is  linear  in  frequency  is  in  keeping 
with  ours  that  the  spatial  rate  of  fading  is  linear  in  frequency 
(see  discussion  in  Section  4.1  following),  i.e.,  that  the  higher 
frequencies  will  yield  proportionately  higher  a number  of  layer 
heights  at  which  deep  fading  will  occur. 

We  now  turn  to  the  work  of  Ruthroff  [2.31]  whose  attempt  to 
analyze  refractive  multipath  was  the  forerunner  of  our  own.  It 
is  easy  to  point  out  the  limitations  of  pioneering  work.  As  an 
example,  we  point  out  that  some  simplifying  assumptions  about 
link  geometry  were  made  (e.g.,  transmitter  and  receiver  antenna 
heights  were  assumed  to  be  equal).  Additionally,  and  perhaps 
more  importantly,  a determination  of  ray  path  amplitudes  was  not 
attempted  [an  absolute  requirement  of  system  function  modeling 
as  evident  from  Eq.  (2.119)].  It  is,  no  doubt,  because  of  this 
last  limitation  that  the  author  makes  the  following  erroneous 
supposition : 

"In  the  present  model  the  amplitudes  of  the  signals 
received  via  the  refracted  rays  are  nominally  equal 
to  the  amplitude  of  the  signal  received  in  the 
absence  of  the  refracting  layer." 

This  is  equivalent  to  assuming  the  ai  in  (2.119)  are  all  of 
value  unity.  We  have  found  this  not  to  be  the  case  and  have 
uncovered  substantial  physical  reasons  why  it  should  not  be  the 
case.  Our  analysis  indicates  that,  with  high  probability,  a 
total  of  three  interfering  rays  are  supported — two  within  the 
layer — a third  beneath  the  layer.  The  very  existence  of  the 
layer  multipath  rays  indicates  that  the  layer  is  in  a focussing 
condition  (i.e.,  there  is  an  intersection  of  ray  paths  at  the 
terminals) . This  strongly  indicates  a signal  amplitude  above 
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unity.  Our  calculations  reveal  invariably  that  this  is  the  case 
in  a multipath  fading  situation. 

Ruthroff  goes  on  to  say: 

"Since  the  probability  of  three  rays  of  equal  ampli- 
tude combining  to  form  a deep  fade  is  much  less  than 
the  probability  of  two  rays  of  equal  amplitude  com- 
bining to  form  a deep  fade,  we  conclude  that  the  dis- 
tribution of  attenuation  is  dominated  by  two-ray 
interference. " 

This  contradicts  the  conclusion  (cited  above)  of  other  workers 
who  felt  that  the  complexity  of  their  data  could  not  be  explained 
in  two-ray  terms,  and  our  detailed  analysis  which  indicates  the 
probable  existence  of  three  paths  (or  more) . 
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SECTION  3 

SNR  DEGRADATION  DUE  TO  FADING 


3 . 1 Receiver  Modeling  for  Performance  Evaluation 
3.1.1  Introduction 

In  Section  2.3  we  presented  the  frequency  power  series  (T 
parameters)  model  for  the  propagation  channel.  From  this  we 
derived  the  p-parameter  representation  for  receiver  output  wave- 
forms. This  representation  may  be  expressed  in  terms  of  either 
frequency  or  phase.  In  the  first  case  it  takes  the  form 


e(t)  = o(t)  + p^cp(t)  + P2tD(t)  + P2o(t)(p(t)  (3.1) 

where  (p(t)  is  proportional  to  the  input  frequency  and  6(t)  to  the 
output  frequency.  In  the  second  case  we  have 

e(t)  = o(t)  + p^o(t)  + P2(p(t)  + j p^o^(t)  (3.2) 

where  o(t)  is  the  input  phase  and  0(t)  the  output  phase.  In  this 
section  we  will  use  the  p-series  model  to  calculate  performance 
of  two  modems.  The  first  is  a baseband  modem  operating  with  con- 
ventional FM  radios;  for  this  application  the  form  of  Eq . (3.1) 
will  be  used.  The  second  is  a narrow  band  phase-continuous  FSK 
modem;  detection  in  this  modem  is  based  on  phase  measurement  so 
the  form  of  Eq . (3.2)  will  be  used. 

Before  proceeding  to  discuss  the  two  modems,  it  is  appropri- 
ate to  indicate  the  approach  which  will  be  used  to  evaluate  per- 
formance in  the  presence  of  distortion.  In  essence,  we  will 
consider  a sequence  of  symbols  consisting  of  the  symbol  to  be 
detected  plus  any  adjoining  symbols  which  contribute  to  inter- 
symbol interference.  For  each  of  the  many  possible  sequences, 
the  error  probability  will  be  evaluated  numerically.  The  final 
step  will  be  to  average  these  results  to  find  the  true  average 
error  probability.  This  procedure  is,  of  course,  different  from 
a true  simulation  in  which  samples  of  the  actual  waveform  are 
generated,  processed  through  a mathematical  model  of  the  channel 
and  receiver,  and  errors  counted.  The  disadvantage  of  the  latter 
approach  is  that  only  high  error  rates  may  be  obtained,  since 
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inordinately  long  runs  would  be  required  to  observe  low  error 
rates.  Thus,  the  approach  to  be  used  will  facilitate  calcula- 
tion of  low  error  rates  which  could  not  be  obtained  with  a true 
simulation.  This  is  a valuable  property  in  the  analysis  of 
line-of-sight  modems  where  the  performance  goal  is  a bit  error 
rate  of  5 x 10”  . 

3.1.2  Baseband  Modem 

3. 1.2.1  Operation  in  the  Absence  of  Distortion 

The  DAU  modem  supplies  a quaternary  baseband  waveform  for 
an  FM  transmitter.  This  waveform  consists  of  a sequence  of  band- 
limited  pulses  modulated  by  quaternary  symbols.  Referring  to 
Fig.  3.1,  we  may  write  the  baseband  waveform  as 


CD 

x(t)  = E a^h^(t-kT)  f’  - 


(3.3) 


In  terms  of  the  transmitted  energy  per  bit  and  symbol  dura- 
tion T,  the  FM  modulated  signal  in  complex  notation  may  be  written 


6 J x(t)  dr 


(3.4) 


The  basic  pulse  waveform  used  in  this  modem  is: 
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We  will  adjust  the  amplitude  a so  that  the  peak  instantaneous 
value  of  x(t)  (taking  into  account  pulse  overlap  for  the  worst- 
case  symbol  sequence)  is  unity  and  hence  the  peak  instantaneous 
frequency  deviation  will  be  6;  this  value  turns  out  to  be  a =0.685. 

At  the  receiver,  the  instantaneous  frequency  at  the  output 
of  the  discriminator  is  given  by  6x(t).  This  is  passed  through 
a receiver  filter  identical  (except  for  gain)  to  the 

transmit  filter: 


Hj^(w)  = ^ Hj^(aj) 


(3.7) 


In  cascade,  Hj^(oo)  and  Hj^(co)  give  the  so-called  "507o  roll-off" 
bandlimited  waveform: 


G(u:)  = H^(a;)H^(co)  = ( 


1 + sin  Iw'  T 


g(t)  = h^(t)  *\(t) 
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(3.9) 


In  a quaternary  modem  with  equally  likely  symbols,  the 
average  error  probability  may  be  written  as 

" 4 {error  Isj^  = -a  j ^ Pr  {error  Is^^  " ' f} 

+ ^ Pr  |errorlsj^  = + ^ Pr  { error  (sj^  = a|  (3.10) 
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Taking  into  account  the  modulation  index  6,  the  proper  decision 
thresholds  for  the  receiver  are  0 and  ± 26a/3.  Thus,  if  yj^  is 
the  sampled  receive  filter  output  due  to  signal  and  n^^  the  noise 
component,  the  probabilities  appearing  in  (3.10)  are  given  by 


Pr  {errorlSj^  - -a|  - Pr  | - -a|  (3.11) 

Pr  {error's^  . . . Pr  | y^^  + . . || 

+ "k  " °'®k  " - f ) (3.12) 


k ■=  °'=k  ■ 3 


Pr  I error  ^ ^ - Pr  | yj^  + n 

+ Pr  I y^  + n^^  > f } 


and 


Pr  {errorlsj^  = a}  Pr  { y^^  + n^  < = a} 


(3.13) 


(3.14) 


Recognizing  that  the  noise  will  be  zero-mean  and  Gaussian,  we  may 
combine  these  equations  as 
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where  '.^(x)  is  the  probability  of  a Gaussian  exceeding  x times 
the  variance: 
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Q(x) 


(3.16) 


Note  also  that  the  argument  of  y denotes  the  transmitted  s3anbol  Sj^ 

The  response  of  Eq.  (3.9)  has  the  Nyquist  property  of  no 
intersymbol  interference  when  sampled  correctly;  that  is, 


g(kT) 


k = 0 
k ^ 0 


(3.17) 


Thus,  in  the  absence  of  channel  distortion,  the  signal  component 
in  the  receiver  output  is  given  simply  by 


y(Sj^)  = 6sj^  (3.18) 

When  this  is  true  Eq . (3.15)  may  be  simplified  greatly;  noting 
that  all  the  arguments  on  the  right-hand  side  are  identical,  it 
becomes 


P 

e 


(3.19) 


In  calculating  the  noise  variance  at  the  receiver  filter 
output  we  will  use  the  quadrature  approximation.  Thus,  if  n(t) 
is  the  complex  envelope  of  the  noise  in  the  receiver  front-end, 
the  phase  noise  is  given  by 


Itn  [n(t)] 


(3.20) 


For  a two-sided  spectral  density  of  Nq/2  watts/Hz,  the  spectral 
density  of  the  complex  low-pass  waveform  n(t)  is  2No  and  that  of 
Im  (n(t)}  is  Nq  so  that  the  density  of  the  phase  noise  is 
NqT 


3-6 


(3.21) 


The  discriminator  output  is  proportional  to  instantaneous 
frequency: 


-I-  A. 

2ff  dt 


[n(t)l 


(3.22) 


and  will  have  spectral  density  given  by 


(211)2 


(3.23) 


Thus,  the  noise  variance  in  the  filter  output  is  given  by 
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2ff 


(3.24) 


Evaluation  of  this  integral  using  the  transfer  function  of 
Eq.  (3.7)  gives  the  following  result 


2 


a 


N 


0 


4(2TrT)^lTQl2E^ 


(3.25) 


Substituting  this  result  into  Eq . 


(3.19)  gives  the  result: 


0.0555  (2ff6T)^  'T„’^  E^. 

U b 

1/2 

j 

^0 

(3.26) 


In  evaluating  performance,  a value  of  6 must  be  selected  to 
approximate  the  value  used  in  the  actual  modem.  This  value, 
which  was  selected  to  produce  a 99?„  bandwidth  equal  to  half  the 
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bit  rate,  is  somewhat  unclear.  In  making  the  calculations  we 
have  assumed  a transmitted  bit  rate  of  27.648  M bit/sec  (the 
highest  rate  of  which  the  modem  is  capable)  and  a peak-to-peak 
deviation  of  8.1  MHz.  With  this  value  2it6T  = 1.85  and  Eq.  (3.26) 
becomes 


P 

e 


(0.190|TqI 


(3.27) 


-9  2 

This  gives  an  error  probability  of  5 x 10  at  \Tq  \ E^/Nq  >=*22.5 
dB  (where  (Tg  p E^/Ng  is  the  value  of  E^/Ng  in  a flat  fade),  which 
corresponds  well  with  both  other  predictions  and  measured  results 
(cf.  Fig.  19  of  [3.1]). 


3. 1.2. 2 Performance  in  the  Presence  of  Multipath 

We  will  now  consider  the  process  by  which  performance  is 
calculated  in  the  presence  of  multipath.  Our  basic  p-parameter 
model  relates  the  time  derivations  of  input  and  output  phase: 


e(t)  = o(t)  + p^o(t)  + P20(t)  + p^O  '(P  (3.28) 

We  let  f(t). denote  the  signal  component  of  the  discriminator 
output,  or  0(t)/2iT.  Making  use  of  the  transmitted  phase 
definition  of  Eq . (3.4),  the  p-parameter  model  becomes 

f(t)  = 6[x(t)  + p^x(t)  + P2x'(t)  + 2Tr6p2x(t)x(t)]  (3.29) 

Note  that  x(t)  is  given  by  Eq.  (3.3);  consequently  the  deriva- 
tives appearing  in  (3.29)  are  given  by 

x(t)  = L a^h^(t  - kT)  (3.30) 

and 

x’(t)  = r a,^h^(t  -kT)  (3.31) 

Referring  to  Fig.  j-1,  the  waveform  at  the  sampling  point  is  the 
result  of  passing  f(t)  through  the  receive  filter,  or 
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y(t)  = f(t)  *hj^(t)  (3.32) 

Using  Eqs.  (3.29),  (3.3),  (3.30),  (3.31),  and  (3.9)  this  becomes 
y(t)  = 6 L a^[g(t  -kT)  + p^g(t  -kT)  + P2g(t  -kT)] 

+ 6(2tr6)p3[(L  a^^h^(t  -kT))(L  a^h^(t  - 4T))]  *hj^(t) 

(3.33) 

The  first  term  in  this  expression  may  be  readily  calculated 
since  g(t)  and  its  derivatives  are  known.  Since  g(t)  falls  off 
rapidly,  only  a few  symbols  influence  y(t).  In  the  computer 
evaluation  of  performance  two  symbols  preceding  and  two  symbols 
following  the  one  being  detected  were  considered,  giving  4^  = 256 
individual  cases.  On  the  other  hand,  the  second  term  must  be 
calculated  by  numerical  convolution.  For  each  of  the  256  cases 
(representing  distinct  sequences  of  adjacent  symbols)  the  error 
probability  may  be  calculated  from  Eq . (3.13);  the  average  error 
probability  is  found  by  averaging  over  these  cases.  However, 
before  making  this  calculation  we  must  specify  the  instant  of 
time  kT  + tg  at  which  y(t)  is  sampled  to  detect  aj^.  This  problem 
is  considered  below. 

In  calculating  receiver  timing  we  assume  a tracking  loop 
which  responds  to  the  mean  zero  crossing  point  of  the  waveform 
y(t).  Consider  the  T second  interval  between  t = kT  and 
t = (k  + l)T.  (We  are  neglecting  delay  in  order  to  simplify  the 
discussion.)  If  a]^  and  a^+i  are  of  opposite  polarity  the  undis- 
torted waveform  x(t)  will  have  a zero  crossing  in  this  interval. 
Although  the  position  of  this  zero  crossing  will  depend  on  the 
data  sequence,  its  mean  value  (averaged  over  all  data  sequences) 
will  be  exactly  at  the  mid-point  of  the  interval.  The  distorted 
waveform  y(t)  will  also  have  a zero  crossing  in  this  interval  as 
long  as  the  "eye"  of  the  waveform  has  not  closed.  Since  eye 
closure  will  result  in  catastrophic  performance  degradation,  we 
will  not  seriously  restrict  the  analysis  by  assuming  the  eye  is 
open.  In  finding  the  mean-zero  crossing  two  preceding  and  two 
following  symbols  are  considered.  Of  these  256  cases  only  128 
will  have  a zero  crossing  in  view  of  the  requirement  on  the 
polarity  of  aj^  and  For  each  of  these  cases,  the  zero 

crossing  of  y(t)  [as  given  by  Eq . (3.33)]  is  found  and  the  results 
are  averaged. 
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3.1.3  IF  Modem 


3. 1.3.1  Principles  of  Operation 

The  modulation  technique  used  in  the  Harris  IF  modem  is 
phase-continuous  FSK.  The  only  unconventional  feature  of  the 
modulator  is  that  the  deviation  is  set  to  the  unusually  low 
value  of  1/8  to  achieve  a packing  density  of  2 bits/sec/Hz.  In 
the  receiver  portion  of  the  modem,  however,  a very  unconven- 
tional structure  has  been  implemented  to  achieve  high  performance. 
In  order  to  provide  background  necessary  to  understand  the  model 
used  for  performance  evaluation,  we  must  briefly  describe  the 
principle  of  operation  of  this  receiver.  For  a more  complete 
description,  as  well  as  for  details  of  implementation,  the  reader 
should  consult  reference  [3.2]. 

For  a modulation  index  of  1/8,  the  separation  between 
adjacent  tone  frequencies  is  1/8T,  where  T is  the  symbol  dura- 
tion. Thus,  the  set  of  four  frequencies  (referenced  to  the  low- 
est) are 


f . = ^ i = 0,1, 2, 3 (3.34) 

Considering  a single  baud,  the  transmitted  phase  waveform  assum- 
ing a transmitted  frequency  f^  will  be 

Oq  + i J (§)  0 < t < T (3.35) 


where  we  have  assumed  a starting  phase  oq  and  are  again  referenc- 
ing to  the  lowest  frequency  fQ.  It  is  clear  from  (3.35)  that  the 
phase  increment  over  a baud  will  be  a multiple  of  7r/4  and  hence 
the  allowable  phase  at  any  symbol  boundary  will  be  one  of  eight 
possible  values  spaced  at  intervals  of  ff/4  around  a circle. 

Assuming  the  signal  structure  described  above,  if  one  could 
measure  the  phase  at  each  symbol  boundary  with  a reference  at 
frequency  fQ  one  could  demodulate  the  transmitted  frequencies  by 
comparing  adjacent  phase  measurements.  In  essence,  this  is  the 
approach  used  in  the  demodulator  in  question.  However,  the  pre- 
sence of  a bandpass  filter  at  the  receiver  input  (Fig.  3.2) 
introduces  variations  in  the  received  phase  which  depend  on 
adjacent  frequencies.  This  filter  is  a fourth-order  Bessel  filter 
with  noise  bandwidth  equal  to  the  symbol  rate,  so  that  (in  the 
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Figure  3.2  Simplified  Modem  Block  Diagram 


absence  of  other  degradation)  the  phase  noise  variance  is  given  by 


(3.36) 


In  summary,  if  the  phase  at  t=0  is  zero  before  the  BPF,  the 
phase  at  its  output  is 


(k-j)eQ  k,  j=0,l,2,3 


(3.37) 


where  f.  is  the  frequency  during  the  interval  (-T,0),  fj^  the  fre- 
quency ^during  (0,T)  and  is  5°.  As  a result  of  this  phase 
distortion,  a decision  as  to  whether  phase  0°  or  one  of  its  near- 
est neighbors  at  ±45°  was  sent  becomes  more  difficult.  The 
correct  threshold  for  making  this  decision  depends  on  the  adjacent 
frequencies  fj  and  fj^.  which  are  unknown  to  the  receiver.  To 
circumvent  this  difficulty,  the  FSK  demodulator  uses  an  iterative 
procedure.  In  essence,  a series  of  tentative  phase  decisions  are 
made  using  compromise  thresholds.  These  tentative  decisions  are 
used  for  selecting  proper  thresholds,  and  final  phase  decisions 
are  made.  Both  simulations  and  tests  on  the  breadboard  receiver 
[3.2]  indicate  that  performance  of  this  ad-hoc  procedure  is 
essentially  equivalent  to  that  calculated  for  making  a noisy  phase 
measurement  with  correctly  chosen  thresholds. 

Other  aspects  of  the  FSK  receiver  which  are  of  significance 
for  performance  evaluation  are  the  extraction  of  references  for 
phase  measurement  and  symbol  timing.  As  shown  in  Fig.  3.2, 
an  eighth-order  nonlinear  operation  is  followed  by  PLL's  at  Sfg 
and  8f^.  The  former  supplies  the  reference  for  phase  measure- 
ment, while  the  difference  between  the  two  supplies  the  symbol 
timing.  The  effect  of  multipath  on  these  operations  will  be 
discussed  in  the  next  section. 

3. 1.3. 2 Effect  of  Multipath 

As  described  above,  the  fundamental  operation  in  the  FSK 
receiver  is  the  measurement  of  phases  at  the  symbol  boundaries. 
Thus,  we  will  analyze  the  phase  measurement  process.  In  doing  so, 
we  will  assume  that  the  receiver  has  correctly  detected  the  two 
adjacent  phases.  Thus  we  will  not  need  to  deal  with  the  com- 
plexity of  the  iterative  procedure  employed  by  the  actual  receiver. 
On  the  other  hand,  we  know  that  actual  receiver  operation  closely 
approaches  this  ideal  behavior. 
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Consider  the  following  situation:  the  receiver  has  correctly 
detected  the  transmitted  phases  C5(-T)  and  (p(T) . With  no  loss  of 
generality,  we  may  assume  that  the  true  phase  at  t=0  is  0.  Thus 
two  error  events  are  of  significance: 

Event  A:  Receiver  decides  o(0)  = ^ 

Event  B:  Receiver  decides  ip(0)  ~ ~ 

In  the  absence  of  multipath  the  actual  phase  y at  the  bandpass 
filter  output  depends  in  a simple  manner  on  the  preceding  and 
succeeding  frequencies  fj  and  f^^: 

y(0)  = (k-j)0Q  (3.38) 

For  o(0)  = the  preceding  and  succeeding  frequencies  would  have  to 
be  fj+i  and  fk_2»  giving  an  actual  phase 

y(0)  Ia  " 4 (3.39) 

So  that  the  correct  threshold  to  be 
o(0)  = 0 and  o(0)  = tt/4  is 


\ [r(0)  + y(0)  = I + (k  - j - 1)6^  (3.40) 

Thus  the  probability  of  event  A is 


PaCJ  »^)  = Pt  I y(0)  4-n>^  + (k-j  - i)0Q| 


= 1 


+ (k  - j - 1)9^  - y(0) 


(3.41) 


2 

where  the  noise  variance  a has  been  defined  above.  Similarly, 
event  B requires  that  the  adjacent  frequencies  be  fj-i  and  ^k+l ’ 
so  that  in  the  absence  of  distortion 
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(3.42) 


y(0)  Ig  = - f + (k  - j +2)0q 


Thus,  for  the  decision  between  o(0)  = 0 and  ^(0)  = -vik  the 
correct  threshold  is 


\ [r(0)  + y(0)  Ig]  = - f + (k  - j +1)0^  (3.43) 

and  the  probability  of  event  B is 

Pg(j  ,k)  = Pr  I y(0)  + n < - I + (k  - j - DcPq} 


= 1 


y(0)  + f ■ (k  - j + Doq 


(3.44) 


To  obtain  average  error  rates  these  expressions  must  be  averaged 
over  j and  k,  taking  into  account  the  fact  that  frequencies  out- 
side the  range  (0,3)  are  not  allowed.  Thus  we  have 


2 3 3 2 

Pr  (phase  error]  S L P^(j,k)  S E Pg(j,k) 

j=0  k-1  j=l  k=0 

(3.45) 

Since  each  phase  error  will  result  in  two  frequency  errors,  the 
symbol  error  rate  (and  the  bit  error  rate,  if  grey  coding  is 
used)  will  be  twice  this  quantity. 

In  the  presence  of  multipath,  the  simple  expression  for 
y(0)  given  above  [ Eq . (3.38)]  cannot  be  used  in  evaluating  P^(j ,k) 
and  Pg(j ,k) . First  of  all,  the  actual  sampling  time  will  not  be 
zero;  but  will  be  determined  by  the  operation  of  the  receiver 
reference  extraction  as  will  be  described  later.  Secondly,  the 
calculation  of  the  received  phase  y(tsamp)  raust  take  into  account 
the  multipath  distortion.  Referring  to  Fig.  3.3,  we  observe  that 
since  the  receiver  BPF  and  channel  model  are  linear  operations, 
we  may  interchange  them.  Although  we  have  previously  discussed 
0(t)  as  the  phase  relative  to  a carrier  at  Fq,  it  will  be  con- 
venient in  what  follows  to  refer  phases  to  a hypothetical  cent.i 
frequency.  With  this  reference,  the  transmitted  phases  ar. 
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